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ABSTRACT 


This  report  documents  work  carried  out  largely  over  the  fifth  and  final  year  of  the 
ONR  sponsored  University  Research  Initiative  (URI)  entitled  “Materials  for  Adaptive 
Structural  Acoustic  Control.”  This  program  has  continued  to  foster  the  successful 
development  of  new  electroceramic  single  crystal  and  composite  material  combinations  for 
both  sensing  and  actuation  functions  in  adaptive  structural  systems. 

For  the  classical  perovskite  relaxor,  dielectrics  typified  by  lead  magnesium  niobate, 
continuing  studies  of  properties  in  the  temperature  region  above  the  dielectric  maximum  T^ 
have  added  strong  additional  support  to  the  superparaelectric/spin  glass  model  for  the 
behavior  developed  earlier  in  the  IMRL.  The  most  exciting  and  important  discovery  of  the 
year  has  been  the  ultra  high  strain  capability  of  relaxor  ferroelectric  single  crystal  actuators. 
For  crystal  in  the  lead  zinc  niobate:lead  titanate  (PZN;PT)  solid  solution  system,  at 
compositions  in  the  rhombohedral  phase  close  to  the  morphotropic  phase  boundary  to  the 
tetragonal  ferroelectric  phase  at  9  mole  %  PT  in  PZN,  crystals  cut  and  poled  along  the  001 
cube  axis  exhibit  massive  field  induced  quasi  linear  anhysterestic  strains  up  to  0.6%.  For  this 
poling  d33  values  up  to  2,300  pC/N  and  coupling  coefficients  k33  of  94%  have  been  achieved 
and  it  was  the  original  hypothesis  that  these  extreme  numbers  must  be  largely  due  to  extrinsic 
domain  wall  motion.  Now  however  it  is  very  clear  that  the  exact  equivalence  of  the  effect  of 
an  001  oriented  E  field  on  the  lll,ill,lil,andni  rhombohedral  domains  precludes  this  field 
from  driving  domain  waU  motion  so  that  quite  contrary  to  our  earlier  expectation  the 
polarization  and  associated  strain  phenomena  are  purely  intrinsic.  At  higher  field  levels  there 
is  an  obvious  step  in  both  polarization  and  strain  into  an  induced  tetragonal  phase  which  gives 
total  reproducible  induced  strains  up  to  1.7%.  Clearly  the  PZN:PT  crystals  represent  a  major 
breakthrough  into  a  completely  new  regimen  for  piezoelectric  actuation  and  sensing. 

For  antiferroelectric:ferroelectric  switching  compositions  in  the  lead  lanthanum 
zirconated  titanate  stannate  family,  new  experimental  studies  have  proven  that  the  induced 
polarization  P3  and  the  strain  X33  onset  at  different  field  levels.  A  new  domain  re-Orientation 
model  has  been  invoked  to  explain  this  startlingly  unusual  behavior.  Both  barium  and 
strontium  additives  have  also  been  explored  to  control  hysteresis  between  forward  and 
backward  switching  with  good  success.  As  well  as  being  interesting  for  transduction  we 
believe  these  compositions  are  sure  to  be  important  for  energy  storage  dielectrics. 

In  composite  sensing  it  is  pleasing  to  report  that  the  moonie  flextensional  patent  has 
now  been  licensed  to  the  Input:Output  Corporation  who  have  successfully  fabricated  and  sold 
more  than  80,000  moonie  sensors.  Work  is  continuing  on  the  cymbal  type  modification  of 
the  moonie  with  focus  now  on  array  structures  for  large  area  panels.  This  topic  is 
transitioning  to  a  joint  study  between  the  IMRL  and  Penn  State’s  ARL,  on  a  new  MURI 
initiative.  For  the  very  small  hollow  PZT  spheres  produced  by  blowing,  the  emphasis  has  been 
upon  both  poling  and  driving  from  outer  surface  electrodes,  and  exploring  both  by 
experiment  and  by  finite  element  theoretical  methods,  the  resonant  mode  structures  which  can 
be  induced.  Studies  of  the  2:2  composite  structures  confirm  the  very  high  effective 
hydrostatic  sensitivity  and  are  permitting  closer  consonance  between  measurement  and 
theoretical  analysis. 

Actuation  studies  have  been  dominated  by  the  initial  exploration  of  the  fantastic  strain 
capability  of  the  relaxor  ferroelectric  MPB  single  crystals.  Obviously  the  induced  strains  are 
on  order  of  magnitude  larger  than  for  conventional  PZT  ceramics,  but  the  blocking  force  has 


not  yet  been  determined.  It  is  expected  that  djj  will  also  be  large  and  anhysteritic  in  these 
crystals,  as  spontaneous  strain  depends  on  Q44  which  is  a  pure  shear  constant.  The  dj5 
however  may  be  significantly  more  complex  as  an  Ej  field  will  certainly  drive  domain  walls  in 
these  E3  poled  crystals. 

Reliability  studies  of  conventional  actuators  are  continuing  with  emphasis  on  using 
acoustic  emission  to  explore  and  separate  domain  wall  motion  and  crack  propagation.  Most 
earlier  studies  were  indeterminate  and  difficult  to  interpret,  recently  for  these  strongly 
piezoelectric  samples  we  have  shown  that  electrical  noise  in  the  power  supply  induces  very 
strong  mechanical  noise  in  the  sample  giving  high  spurious  emission  counts.  New  studies 
using  a  long  time  constant  filter  in  the  supply  have  permitted  clear  and  effective  separation. 
Over  the  last  few  years  there  has  been  a  strong  re-awakening  of  interest  in  bimorph  type 
transducer  amplitiers  with  new  concepts  like  rainbow,  cerambow  and  thunder  appearing. 
Under  our  ONR  program  with  Virginia  Polytechnic  it  has  been  necessary  to  sort  out  the 
conflicting  claims  for  these  ‘morph’  types  and  these  data  are  included  for  completeness.  We 
have  also  begun  serious  study  of  the  large  electrostriction  in  the  soft  polyurethane  elastomers 
where  it  has  been  necessary  to  derive  new  techniques  to  measure  strain  with  ultra  low 
constraint  on  the  films. 

Processing  studies  now  involved  both  single  crystal  flux  growth  and  a  wide  range  of 
powder  and  ceramic  processing.  Current  needs  for  integrity  and  better  mechanical  properties 
are  driving  new  needs  for  fine  grained  PZT  piezoceramics  and  new  processing  is  permitting 
retention  of  excellent  properties  down  to  submicron  grain  sizes. 

From  the  wide  range  of  thin  ferroelectric  film  activities  in  the  laboratory,  only  those 
which  refer  to  the  thicker  films  being  produced  on  silicon  for  MEMS  devices  are  included. 
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Evaluation  of  Piezocomposites  for  Ultrasonic  Transducer  Applications 
Influence  of  the  Unit  Cell  Dimensions  and  the  Properties  of  Constituents  on  the 
Performance  of  2-2  Piezocomposites 

X.  Geng  and  Q.  M.  Zhang 

Materials  Research  Laboratory  and  Department  of  Electrical  Engineering 
Penn  State  University,  University  Park,  PA  16802 


Abstract: 

A  theoretical  model  on  piezocoamic  polymer  cmiqxjsites  with  laminar  periodic  stracture  is 
presented.  A  salient  feature  of  this  model  is  that  it  can  treat  explicitly  how  the  unit  cell 
dimensions  and  other  material  properties  influence  the  performance  of  an  ultrasonic  transducer 
made  of  2-2  piezocomposites.  The  model  predicts  that  there  exist  a  series  of  modes  associated 
with  the  periodic  structure  of  a  composite,  which  is  beyond  the  stop-band  edge  resonance 
prediction.  One  of  the  main  concerns  in  designing  a  con^xisite  transducer  is  how  the  surface 
vibration  profile  changes  with  fiequency  and  how  this  is  influenced  by  the  aspect  ratio  of  the 
ceramic  plate.  It  was  predicted  that  as  long  as  the  thickness  resonance  is  below  the  first  lateral  • 
mode  frequency,  there  is  always  a  fi'equency  f,  which  is  near  the  thickness  resonance  and  at 
which  the  polymer  and  ceramic  vibrate  in  unison.  The  effect  of  aspect  ratio  is  to  change  the 
position  of  f,  with  respect  to  the  thickness  resonance  frequency  and  the  bandwidth  in  which 
polymer  and  ceramic  have  nearly  the  same  vibration  amplitude  and  phase.  It  is  also  predicted 
that  when  operated  in  a  fluid  medium  such  as  water,  there  will  be  a  resonance  mode  which  has 
a  frequency  determined  by  the  velocity  of  the  fluid  medium  and  the  unit  cell  length  d  and  is 
associated  with  the  oscillation  of  the  fluid.  The  behavior  of  a  composite  plate  as  an  acoustic 
transmitter  and  receiver  and  the  influence  of  the  aspect  ratio  of  the  ceramic  plate  on  them  are 
also  investigated. 
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1.  Introduction 

To  produce  a  high  performance  ultrasonic  transducer  requires  the  transducer  material,  which 
performs  the  energy  conversion  between  the  mechanical  form  and  electrical  form,  to  have  a 
high  electromechanical  coupling  factor,  broad  operation  frequency  bandwidth,  and  adjustable 
acoustic  impedance  which  can  be  tuned  to  match  that  of  the  medium.  With  single  phase 
piezoelectric  materials,  it  is  difficult  to  simultaneously  meet  all  the  requirements.  Piezoceramic 
polymer  composite  materials,  which  combine  the  high  electromechanical  activity  of 
piezoceramics  and  the  low  acoustic  inqiedance  of  polymeric  materials,  have  provided  new 
opportunities  to  meet  these  requirements.*’^  Since  their  inception  in  the  seventies,  the 
piezoceramic  polymer  coiiqiosites  have  become  one  of  the  most  impratant  transducer  matmaic 
and  are  being  widely  used  in  many  areas  such  as  medical  imaging,  nondestructive  evaluation  of 
materials,  and  underwater  vision,  etc. 

Being  a  diphasic  material,  tihe  properties  of  a  piezo-conqx>site  can  be  tailored  ovw  a  wide 
range  by  adjusting  the  material  properties  and  geometric  shapes  of  constituent  phases.’-^  It  has 
also  been  observed  that  the  properties  of  a  composite  vary  with  frequency.’  The  diallenge  of 
understanding  the  seemingly  complex  relationship  between  the  performance  of  a  composite  and 
the  properties  of  its  constituents  and  the  great  opportunities  provided  by  these  Tnat«Tii^i<f  have 
stimulated,  in  the  past  two  decades,  extensive  investigations,  both  experimental  and  dieotetical, 
on  this  class  of  materials. 

The  classic  work  of  Newnham  al.,’  which  classified  piezoconqwsite  materials  according 
to  the  connectivity  of  tlw  constituent  phases,  has  greatly  facilitated  the  analysis  of  the 
composites  as  the  connectivity  is  one  of  the  key  parameters  in  determining  the  performance  of  a 
composite.  For  composites  with  1-3  and  2-2  connectivities,  both  analytical  and  finite  element 
modelings  have  been  carried  out  which  have  provided  useful  guidelines  in  the  design  of 
composite  transducers.’’”  The  isostrain  models  developed  by  Smith  et  al.’  and  Hashimoto  et 
al."*  linked  the  material  parameters  of  the  constituents  to  the  effective  piezoelectric  properties  of 
1-3  and  2-2  composites  respectively  and  predicted  that  the  thickness  coupling  factor  k,  of  a 
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composite  can  approach  the  longitudinal  coupling  factor  ife33of  the  piezoceramic  rod  and  of 
the  plate  respectively,  which  is  in  good  agreement  with  experiments  for  composites  with  a  high 
aspect  ratio  t/d,  where  t  is  the  thickness  and  d  is  the  periodicity  of  the  composites.  Auld  et  al., 
using  the  Floquet  theory,  investigated  wave  propagations  in  both  the  2-2  and  1-3  composites 
and  showed  that  due  to  the  periodic  structure  of  these  composites,  there  exist  pass  bands  and 
stop  bands,  similar  to  the  band  structure  in  a  crystal  solid,  and  that  there  are  piezoelectric 
resonances  associated  with  the  stop  band  edge  resonances.''"’  For  the  design  of  composite 
transducers,  the  recognition  of  the  existence  of  these  modes  and  the  precise  piedictitm  of  their 
frequencies  are  of  prime  importance  since  quite  often,  it  is  the  interference  of  these  modes  with 
the  thickness  mode  of  a  composite  that  deteriorates  the  performance,  especially  at  high 
fiequency  operations.  Craciun  et  al.  examined  the  coupling  between  these  latoal  modes  and 
thickness  mode  using  a  phenomenological  approach  and  the  results  provided  qualitative 
understanding  between  the  coupling  of  the  two  modes  and  the  material  properties."  The  results 
from  these  investigations  have  played  important  roles  in  the  development  of  ultrasonic 
composite  transducers.  However,  due  to  the  approximations  used  in  the  analysis,  there  are 
severe  limitations.  For  instance,  various  feature  related  to  the  dynamic  behavior  of  a 
conqmsite  transducer  were  not  treated  in  a  consistent  rhanner,  and  the  effect  of  finite  thickness 
of  a  composite  on  the  material  properties,  that  has  been  sho^^'n  to  be  crucial  in  determining  the 
performance  of  a  piezocomposite,  carmot  be  treated  in  these  analyses.  To  address  these  realistic 
issues  of  a  composite  material,  finite  element  analysis  has  been  employed  by  many  authors.'^*'’ 
For  example,  the  dispersion  curves  have  been  evaluated  for  various  modes  in  a  composite  and 
the  dependence  of  the  electromechanical  coupling  factor  on  the  ceramic  volume  content  and  the 
ceramic  rod  shape  was  investigated. 

For  a  piezoceramic  polymer  composite,  it  has  to  be  recognized  that  it  is  the  ceramic  phase 
which  performs  the  energy  conversion  between  the  electric  and  mechanical  forms  and  the 
polymer  phase  merely  acts  as  a  carrier  which  transfers  acoustic  energy  between  the 
piezoceramic  and  the  external  medium.  Hence,  if  the  elastic  coupling  between  the  two 
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constituents  is  not  veiy  effective,  even  if  the  material  exhibits  a  perfect  acoustic  impedance 
matching  with  the  medium  and  the  electromechanical  coupling  factor  is  large,  the 
electromechanical  performance  of  the  material  is  still  poor.  These  observations  clearly  indicate 
that  in  modeling  piezoceramic  polymer  composites,  one  cannot  singly  use  an  effective  medium 
s^oach  and  has  to  take  into  account  explicitly  this  internal  degree  of  freedom. 

Recently,  we  developed  an  analytical  nxxlel  on  the  dynamic  problem  of  a  piezocomposite 
material  with  the  2-2  connectivity.**"’*  In  the  model,  we  avoided  the  approximations  made  in 
the  earlier  works  and  hence,  can  address  the  dynamic  responses  of  2-2  piezocomposites,  such 
as  the  frequencies  of  various  modes,  the  mode  coupling,  the  electnmiechanical  coiqiling  factor, 
the  vibration  profiles  of  composites  under  different  external  driving  conditions,  etc.,  in  a 
realistic  and  consistent  manner.  For  exanqile,  one  of  the  misconcqitions  in  the  early  studies  of 
the  dynamic  behavior  of  piezocomposites  is  the  direct  linkage  betwera  the  non-uniftmn  surface 
viteation  profile  in  a  composite  and  the  aspect  ratio  of  the  ceramic  plate  (in  2-2  conqiosites)  or 
rod  (in  1-3  composites).  Here,  we  will  show,  which  was  verified  by  experiment,  that  for  a 
composite  plate,  as  Itmg  as  the  thickness  lescmance  frequency  is  below  that  of  the  lateral  mode, 
there  is  always  a  frequency  near  the  thickness  resonance  where  the  vibration  profile  of  the 
composite  is  uniform.  The  influence  of  the  aspect  ratio  is  on  the  frequency  bandwidth  in  which 
the  polymer  and  ceramic  vibrate  in  unison. 

The  approach  taken  here  to  solve  the  vibration  problem  in  a  finite  thickness  composite  plate. 
is  based  on  the  method  of  partial  wave  expansion  where  the  various  elastic  and  electric  fields 
are  expended  in  terms  of  the  eigenmodes  of  the  structure  as  shown  in  figure  1(a)  and  the 
coefficients  for  each  eigen-mode  cuu  determined  by  the  boundary  conditions  at  X3  =  t  /2  (figure 
1(b)).  Some  of  the  results  have  been  presented  in  early  publications.  In  this  paper,  we  would 
like  to  summarize  briefly  these  early  results  and  then,  discuss  many  issues  pertinent  to  the 
design  of  piezocomposite  ultrasonic  transducers  based  on  the  results  from  the  model  analysis. 
The  paper  is  organized  as  the  following;  First,  the  details  of  derivation  of  the  eigenfunctions 
for  a  2-2  piezoceramic  polymer  cotrq)osite  will  be  presented.  Based  on  the  results,  the  wave 
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propagation  in  an  unbounded  2-2  composite  is  analyzed.  For  the  finite  thickness  composites, 
various  features  of  a  piezocomposite  plate  under  an  external  driving  electric  field  in  both  air  and 
water  media  are  treated.  The  response  behavior  of  the  composite  plate  under  external  pressure 
is  also  treated.  Experiments  were  conducted  to  provide  cotrqKirisoD  with  the  theoretical  results. 
From  the  r^ults,  we  will  show  explicitly  how  the  aspect  ratio  t/d  influences  the  transduction 
performance  of  a  composite  plate. 

n.  General  Solutions  for  Wave  Propagation  in  a  2-2  Piezoconposite 

Shown  in  fiigure  1(a)  is  a  schonatic  drawing  of  a  2-2  conposite  (unbounded  in  the  Xj- 
direction),  where  plates  of  piezocoamic  and  polymer  form  a  paraUel  array.  The  coordinate 
system  is  chosen  such  that  the  Xj-axis  is  along  the  ceramic  poling  direction,  the  x,-axis  is 
perpendicular  to  the  ceramic  polymer  interface,  and  Xj-axis  is  in  the  plane  of  the  plates.  For  a 
typical  2-2  composite,  die  dimensions  in  the  x,-  and  x^-directions  are  much  larger  than  die 
period  d  and  thickness  t  In  the  treatment  here,  they  can  be  taken  as  infinite  without  much  error 
in  the  results.  Under  these  conditions,  the  composite  is  clamped  in  the  Xj-direction  so  that 
82=0,  where  Sj  is  the  strain  in  the  Xj-direction,  and  the  problem  becomes  a  two  dimensional 
one  with  no  dependence  on  the  Xj-coordinate. 

The  governing  equations  for  the  dynamics  of  a  2-2  composite  are”'^’ 


Hi 

•t- 

Hi 

+Hi- 

dDi 

II 

=  0 

dx\ 

*3 

The  symbols  adopted  in  this  paper  are  summarized  here:  T,  and  are  the  stress  and  strain 
tensor  components,  where  the  Voigt  notation  is  used,  u,  is  the  elastic  displacement  vector,  p  is 
the  density,  and  D;  is  the  electric  displacement  vector,  E,  is  the  electric  field.  The  relevant 
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material  coefficients  are:  e^^  is  the  piezoelectric  coefficient,  c^^  is  the  elastic  stiffness,  and  e,  is 

the  dielectric  permittivity.  Equation  (1)  holds  for  both  the  polymer  and  ceramic  phases. 

The  constitutive  equations,  relating  stress  T,  strain  S,  electric  displacement  D,  and 
electric  field  E,  are 

[r]  =  [c^][5]-[e,][£]  (2  a) 

[D]  =  [cIS1+[£^][£]  (2b) 

For  the  polymer  phase,  e^  in  the  equation  (2)  are  zero.  The  superscripts  E  and  S  indicate  that 
the  coefficients  are  under  the  constant  E  field  and  constant  strain  conditions,  respectively.  Here 
[e,]  is  the  transposed  [e]  array.  Under  the  quasi-electrostatic  tqiproximation,  die  electric  field  E 
can  be  expressed  as 

E  =  -V4>  (  3  ) 

where  O  is  the  electrical  potential. 

Combining  eqs;  (1),  (2),  and  (3)  yields  three  second  order  differential  equations, 
governing  the  elastic  displacement  u„  Uj,  and  electric  potential  4>  in  the  ceramic  plate, 
respectively: 

4"3.U  +(4  +4)«1.I3  +4“3.33  +(%<I‘.33  +«l5*ll)  = 

4“l.ll  +(4  +4Ki3  +4«1.33  +(«3I  +<?15)‘I».I3  =  P«I  (  4  ) 

1  (^15  ^31  )“l.l3  ^33^3.33  “  (4^,33  4^.1 1 )  —  0 

For  the  polymer  phase,  ejj  in  eq.  (4)  should  be  taken  to  zero. 

For  an  unbounded  composite,  the  solutions  to  eq.  (4)  have  the  form 
U3  =  A  exp(J(hx,  +  -  ax)) 

M,  =  £exp(y(fcc,  +  -ax))  ( 5  ) 

d>  =  Ce\p(j(hx^  +  Px^  -OX)) 
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where  A,  B,  and  C  arc  three  constants,  to  is  the  angular  frcquency,  h  and  P  arc  the  wave 

vector  components  in  the  x,-  and  Xj-  directions,  respectively. 

Substituting  eq.  (5)  into  (4)  yields  three  homogeneous  equations  with  the  undetermined 
constants  A,  B,  and  C, 

=  0  (6) 

where  [Aj]=[A,  B,  and 


w- 

-  pm* 

(<^i  +  Cu)hp 

(«l5+<3l)*^ 

(Clj  +  ^3, 

The  condition  for  a  nontrivial  solutions  is  such  that  the  detemnnant  of  the  coefficient  matrix 
vanishes,  i.e.,  IMjjM  (8) 

Equation  (8)  is  a  cubic  equation  of  hi  For  a  given  ca  and  P ,  eq.  (8).  in  general,  has  three  roots 

of  h^,  denoted  as  h|,  and  b^,  corresponding  to  the  quasi~electromaghetic,  quasi- 

longitudinal,  and  quasi-shear  waves  in  the  piezoelectric  plate,  respectively.  For  each  h‘,  the 
ratio  among  A,  B,  and  C  can  be  determined  ftom  eq.  (6).  Since  we  are  concerned  only  with  the 
waves  in  the  X3  direction  which  correspond  to  piezo-active  modes  in  a  finite  thickness 
composite  plate,  the  general  solutions  become 

“3  =  '^R^f^cos{h‘x^)siniPxJ) 

i 

<  =  sin(A;x,)cos(^3)  ( 9 ) 

r 

cos(h|x^)sinipXy) 

I 

where  i  runs  from  1  to  3.  f^,  gj  and  tj  are  the  cofactors  of  A^.d).  A^jCi),  and  A^3(i)  of  the 
determinant  (8)  (where  h  is  replaced  by  h'  for  i=l,  2,  and  3.  respectively).  The  exp(-ja)t)  term 
in  eq.  (9)  is  omitted. 
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Following  the  similar  procedure,  the  solutions  for  the  polymer  phase  can  be  obtained 
(the  center  of  the  polymer  plate  is  at  x,  s  d/2) 

>4  =  f/'cosih!'  (x^ -^))sinifixj) 

i  ^ 

“T  =  -  7))cos(^)  ( 10 ) 

I  ^ 

coshiP(x^  - ^))sin{Px^) 


where  i  =1, 2.  and  gf  are  the  cofactors  of  Ak,(i),  Aa(i)  of  the  determinant  (8)  with  aD  the 
material  parameters  replaced  by  those  of  the  polymer  phase,  and 

(hP)2  =  (kf)2 .  and  (kf)2  -  p^ 

where  ~  —p ,  kj.  ~  and  are  die  Itxigitudinai  and  shear  wave  velocities  of  the 

polymer  phase,  respectively. 

The  expressions  of  the  stresses  and  the  electric  displacement  in  the  ceramic  plate  can  be 
obtained  by  substituting  eq.  (9)  into.eq.  (2);  :  .  .  . 

=  X  cos(/»f X,  )cos(j3xj) 

i 

Tf  =  cos(hFx,  )cos(fixj)  (11) 

t 

"^5^  =  X^(0/^‘^sin(hfx,)sm(jSbr3) 

I 

where 


T,U0  =  cf,hfgf+(4f+e,,tf)p 

7^3(0  =  C, 3^,  gj  +(^33/  +^33^1  (12) 


where  i=l,  2, 3  for  ceramic  phase.  The  electric  displacements  are: 
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i 

-^D^(i)Rfcosihfx^)cos{fixJ) 


(13a) 

(13b) 


where  ^  + hff)+  ef^hftf  and Dfii)  =  -  4)5  tf.  SimUar 

expressions  can  be  obtained  for  the  polymer  plate.  For  instance,  the  electric  displacements  are: 


D,'*  =  -4^'’sinh()5(jc,  --))sin(^3) 
Dj  = —E^^pc  cosh(^(X|  —  — ))cos()Sxj) 


(14c) 

(14d) 


Making  use  of  the  boundary  conditions  at  the  ceramic  polymo:  intei£K%  (x,=  vd/2,  where  v 
is  the  volume  fraction  of  the  ceramic  in  a  connposite),  which  are 

4=4,  4  =  4,  21^  =  7;'“  and  2f  =  (15a) 

and  and  ( 15^ ) 

and  C**'  (i  =1. 2, 3  and  j=l,  2)  can  be  determined.  Expanding  Eqs.  (15)  in  tenhs  of 
eqs.  (9),  (10),  (11),  (13)  and  (14)  yields  the  following  six  homogeneous  linear  equations: 


i^.3 

^.4 

0- 

^2. 

i^^23 

^25 

0 

^32 

^33 

^35 

0 

^4, 

^42 

^^43 

^44 

0 

0 

0 

Ks, 

^54 

Kss 

0 

0 

i^64 

^65 

(16) 


where  the  matrix  elements  Kjj  are  functions  of  p  ,CD,  d,  v,  and  the  material  parameters  of  both 
the  polymer  and  piezoceramic.  The  condition  that  the  determinant 

Ik,I=0  (17) 
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yields  the  relationship  between  ©  and  p,  the  dispersion  relations  in  the  composite.  For  each 
pair  of  ©  and  p ,  the  relationships  among  Rf,  and  C?*  can  be  obtained.  Hence,  the  various 
stress,  strain,  electric  field  distributions  in  the  composite  can  be  determined. 

in.  Dispersion  Curves,  Modes  and  Modes  Coupling  of  a  2-2  Piezocomposite 

Eq.  (17)  allows  us  to  determine  the  relationship  between  P  and  ©,  the  dispersion  curves,  for 

a  composite  if  the  materials  parameters  of  piezoceramic  and  polymer  and  the  geometric 

parameters,  such  as  d  and  v,  are  known.  Eq.  (17)  is  a  transcendental  function  which  cannot  be 

solved  analjrtically.  A  computer  code  was  developed  and  the  dispersion  curves  were  evafaiated 
numericaUy. 

\ 

Shown  in  figure  2  (the  solid  curves)  are  the  two  lowest  branches  of  tfie  dispersion  carves 
for  a  2-2  coiiqwsite  made  of  PZT-5H  piezoceramk  and  Spurr  epoxy  with  the  ceramic  vdume 
fi'action  of  44%. ”  For  all  the  composites  discussed  in  this  paper,  excqx  otherwise  specified, 
PZr-5H  piezoceramic  and  Spurr  epoxy  are  used  as  the  constituents.  The  parameters  of  PCT- 
5H  and  Spurr  epoxy  are  presented  in  Table  I.  The  general  trend  of  the  dispersion  curves  in 
figure  2  resembles  that  of  the  S3mimetric  Lamb  waves  in  a  plate.*’  As  will  be  shown  later,  at 

small  Pd  limit,  the  first  branch  corresponds  to  the  longitudinal  wave  pr(^)agation  along  the  Xj- 

direction,  that  is,  Uj  is  the  dominant  displacement  which  is  more  or  less  uniform  and  in  phase 
in  the  x, -direction,  and  its  phase  velocity  is  the  effective  longitudinal  wave  velocity  of  the 
composite.  The  second  branch  corresponds  to  the  lateral  resonance  which  arises  from  the 
periodic  structure  of  the  composite,  that  is,  it  is  a  stationary  shear  wave  along  the  x, -direction, 
and  it  can  be  shown  that  the  displacement  u,  of  the  polymer  phase  is  much  larger  than  that  of 

the  ceramic  phase  and  the  phase  difference  between  them  is  180°.  Hence,  it  is  the  so  called 

stop-band  edge  resonance  as  predicted  by  Auld  et  al.^'® 
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Under  the  assumption  that  the  wave  length  X  in  the  Xj-direction  (P  =2n/  X)  is  equal  to  two 

times  the  composite  thickness  t  (the  condition  for  the  thickness  resonance),  the  theoretical 
dispersion  curves  can  be  conqnred  with  the  experimental  results,  obtained  from  composites 

with  different  thickness.  At  low  Pd  on  the  dispersion  curves,  where  d  is  the  period  of  die 

composite,  the  frequency  of  the  fundamental  thickness  mode  from  the  experiment  falls  on  the 
first  branch  as  marked  by  the  Mack  dots,  and  the  frequency  of  the  first  lateral  mode  falls  on  the 
second  branch  as  marked  by  the  open  circles.  Afro-  the  crossover  region  B,  the  modes 
interchange  the  positions  on  die  dispersion  curves  where  the  thickness  mode  is  in  the  seomd 
branch  while  there  is  a  weak  icstMiance  at  a  frequency  below  the  thickness  mode  and  it  falls  on 
the  first  branch.  Here,  the  thidmess  mode  is  defined  as  the  one  widi  higher  electromechanical 

coupling  factor.  The  mode  on  the  first  branch  gradually  diminishes  at  high  pd  values.  The 

result  is  presented  in  figure  2.  Qearly,  there  is  an  excellent  accord  between  the  theoretically 
derived  resonance  frequencies  and  experimentally  observed  ones. 

Shown  in  figures  3(aL)"and  figure  3(b)  are'  the  3-dimensional  dispersion  curves  for  2-2 
coinposites  with  15%  and  44%  ceramic  volunK  fraction,  respectively.  Obviously,  there  arc 
imaginary  and  complex  branches  of  the  dispersion  curves.  The  modes  on  these  branches  are 
non-propagating  modes  (imaginary  branches)  and  attenuated  modes  (complex  branches), 
respectively.  These  modes  do  not  exist  in  an  unbounded  composite,  however,  they  are 
important  in  the  vibration  problems  of  finite  thickness  plates  as  well  as  in  semi-infinite 
mediums,  where  they  correspond  to  the  evanescent  waves  at  the  surface  which  will  be 
discussed  later  in  the  paper. 

Figure  4  presents  the  effect  of  the  stiffness  of  the  polymer  phase  on  the  dispersion  curves  of 
a  2-2  composite  with  30%  ceramic  content  where  PZT-5H  is  used  as  the  piezoceramic  phase.  It 
is  apparent  that  the  frequency  position  of  the  second  branch  is  very  sensitive  to  properties  of 
the  polymer  phase.  This  is  quite  understandable  since  as  will  be  shown  later,  the  frequency 
position  of  this  branch  is  directly  related  to  the  shear  velocity  of  the  polymer  matrix  for 
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composites  at  this  ceramic  content.  The  lower  the  shear  wave  velocity  of  the  polymer*phase  is, 
the  lower  the  frequency  of  the  second  branch  will  be.  Therefore,  for  a  transducer  operated  at 
high  frequencies,  in  order  to  avoid  the  interference  from  the  lateral  modes,  a  polymer  matrix 
with  a  high  shear  velocity  should  be  utilized  even  though  the  thickness  coupling  factor  k,  of  this 
composite  may  be  reduced  as  a  result  of  the  stiffer  polymer  matrix.^-^ 

In  order  to  further  elucidate  the  origins  of  these  resonant  modes,  it  is  instructive  to  maift*  a 
comparistNi  between  the  dispersion  curves  of  a  composite  and  those  of  single  piezoceramic  and 
polymer  plates  with  appropriate  boundary  conditions.®  In  general,  in  piezoceiamic  polymer 
conq)osites,  die  elastic  stifbiess  of  the  piezoceramic  is  mote  than  one  ord^  of  magnitude,  higlytr 
than  diat  of  die  polymer.  It  is  reasonable  to  assume  that  the  piezoceramic  plates  in  a  2-2 
composite  are  stress  free  at  the  ceramic-polymer  int^ace,  that  is,  the  stress  components  T,  and 
Tj  ate  zero  at  the  interface.  Similarly,  the  polymo*  plates  can  be  tqiproxitnated  as  under  the 
fixed  boundary  conditions,  that  is,  u,  and  u,  are  zero  at  the  interface.  Under  these  assumptions, 
the  dispersion  curves  of  both  single  ceramic  and  single  polymer  plates  are  calculated  for  the 
plates  with  different  widths  corresponding  to  2-2  composites  with  different  coamic  volume 
fractions.  The  results  are  presented  in  figures  5(a),  5(b),  and  5(c)  corresponding  to  2-2 
composites  with  15%,  44%  and  .80%  ceramic  content,  respectively.  In  these  figures,  the  solid 
lines  are  the  dispersion  curves  of  the  2-2  composites,  the  dashed  lines  are  those  of  the  ceramic 
plate,  and  doted  lines  are  the  dispersion  curves  of  the  polymer  plate. 

The  results  from  these  figures  reveal  that  there  are  many  resemblance  between  the 
dispersion  curves  of  2-2  composites  and  the  dispersion  curves  of  the  single  ceramic  and  single 

polymer  plates  with  appropriate  boundary  conditions.  For  instance,  at  small  Pd,  the  first 

branch  in  the  dispersion  curves  of  the  2-2  composites  with  44%  and  80%  ceramic  volume 
content  is  vei>'  close  to  the  first  branch  of  the  ceramic  plate.  On  the  other  hand,  for  2-2 
composites  with  low  ceramic  volume  content,  a  large  difference  between  these  two  is  found  for 
the  first  branch.  These  results  are  consistent  with  those  from  Smith  et  al.  and  from  Hashimoto 
et  al.  based  on  an  effective  medium  model.^'^  For  the  second  branch  which  corresponds  to  die 
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lateral  resonant  mode  in  a  composite  plate,  at  low  and  medium  ceramic  volume  content,  it  is 
close  to  the  first  branch  of  the  ptriymer  plate  which  frequency  at  small  Pd  is  equal  to  V/  lldp 

where  and  dp  are  shear  wave  velocity  and  width  of  the  polymer  plate,  respectively. 
However,  for  composites  with  high  ceramic  volume  content  such  as  the  one  shown  hoe 
(80%),  the  second  branch  is  related  to  the  longitudinal  resonance  of  the  ceramic  plate  along  the 
width  (or  x,-)  direction  while  the  shear  resonance  associated  with  the  polymer  plate  lies  on  the 
third  branch  of  the  dispersion  curves. 

Obviously,  the  coupling  between  the  two  phases  through  the  interface  boundary  conditions 
will  influence  the  dispersion  curves  of  the  waves  in  tfie  two  phases.  It  is  well  known  diat  die 
dispersitm  curves  for  the  uncoupled  waves  are  split  at  their  crossover  points  when  coupling  is 
introduced."  Far  from  the  crossover  region,  the  coupled  wave  dispersion  curves  should  nearly 
coincide  with  those  of  the  uncoiqrled  waves.  If  the  coupling  is  very  strong  as  for  the  15%  and 
44%  piezocon^sites,  the  coupled  waves  exhibit  large  departure  from  the  uncoupled  curves  in 
the  crossover  region,  which  is  dearfy  shown  in  figure  5(a)  and  figure  5(b).  While  for  the 
composite  with  80%  ceramic  content,  the  coupling  between  the  first  and  the  second  branches  is 
through  the  coupling  of  P  wave  aiKl  SV  wave  in  the  ceramic  plate  where  the  interface  does  not 
have  a  significant  effect  on  it. 

IV.  Vibration  of  a  Finite  Thickness  (Composite  Plate  under  an  Electric  Field  in  Air 

In  the  previous  section,  the  properties  of  guided  wave  propagation  in  laminatpH  2-2 
piezocomposites  have  been  analyzed.  In  spite  of  the  fact  that  many  effective  parameters  of  the 
material  can  be  derived  by  this  simple  method,  one  key  issue  in  the  design  of  a  composite 
transducer,  i.e.,  the  influence  of  die  aspect  ratio  of  the  ceramic  plates  (or  the  unit  cell)  in  a 
composite  on  the  performance  of  the  transducer  cannot  be  addressed.  In  this  section,  we  will 
treat  the  vibration  problem  of  a  finite  thickness  composite  plate  under  an  external  driving 
electric  field  and  situated  in  air.  ftom  the  analysis,  one  can  obtain  detailed  information  on  how 
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the  surface  vibration  profile  changes  with  frequency  and  its  dependence  on  the  aspect  ratio  of 
the  ceramic  plate  (or  the  unit  cell  dimension),  the  possible  resonant  modes  in  a  composite 
transducer,  and  the  dependence  of  the  electromechanical  coupling  factor  on  the  aspect  ratio  of 
the  ceramic  plate,  etc.  We  will  also  show  that  as  long  as  the  thickness  resonance  ftequency  is 
below  the  lateral  mode  frequency,  the  aspect  ratio  will  not  have  a  direct  effect  on  the  vibration 
uniformity  of  a  composite  near  the  thickness  resonance.  The  influence  is  on  the  bandwidth  in 
which  the  ceramic  and  polymer  vibrate  with  nearly  same  amplitude  and  phase. 

A  composite  plate  with  a  thickness  t  is  drawn  schematically  in  figure  1(b).  For  the  problem 
treated  here,  the  two  free  surfaces  of  the  composite  are  electroded  with  conducting 
material  and  an  AC  electric  field  of  a  frequency  f  is  tq)plied  between  the  two  electrodes. 

It  is  well  known  that  there  exist  no  simple  solutions  for  the  vibration  problem  of  a  finite 
thickness  plate  such  as  the  piezocon^site  treated  here.  One  of  die  most  frequently  used 
approximation  mediods  is  to  expand  the  elastic  and  electric  fields  in  a  material  in  terms  of  the 
eigenfunctions  in  an  unbounded  one.**  Different  techniques  such  as  the  variational  technique 
and  the  method  of  least  squares  can  be  used  to  detennine'the  expansion  coefficients.  For  the 
conqiosite  plate  treated  here,  we  found  that  the  variational  technique  is  more  qipropriale  in 
treating  the  boundary  problem  than  the  method  of  least  squares  because  of  die  large  difference 
in  the  vibration  amplitudes  between  the  two  phases  in  the  composite. 

For  the  ceramic  phase  in  the  composite  plate,  the  elastic  displacements  Uj  and  u,,  and  the 

electric  potential  are  expended  in  terms  of  the  eigenfunctions  derived: 

11  =  1  1  =  1 

wf  =  +  CR^sm{h^,x,)  ( 18  ) 

#1=1  #=l 

cos(/i„^x,)sin(^„X3)/l,  +  Oc, 

n=I  i=l 

Similarly,  for  the  polymer  phase  in  the  composite  plate: 
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(19) 


cos(/»^(x,  -^))sin(Aar3)M, 

#1*1  •*!  ^ 

■«,'  =  -r))cos(Aj,H.  +  C/lf  sinCAlCi,  -|)) 

#1*1 1*1  ^  2 


‘J»^=  Xqrcosh()3„(x,-^))sin()5„X3)^  +  Car3 

/i=l  2 

where  /^‘j,  /<f,  and  are  determined  in  eq.  (16).  p„,  and /i'  are  the 

wave  vector  conqx>nents  of  the  nth  mode  in  the  X3-  and  x, -directions  where  the  superscripts  c 
and  p  stand  for  the  ceramic  and  polymer,  respectively.  ^  and  are  related  to  toe  mode  in 


which  P  is  equal  to  zero:  and  *  1^©.  R§  and  Rg  are  determined  by  eq. 

\  V  ^11 


-  P 


(16)  in  which  P  is  set  to  zero.  The  mode  with  P=0  is  generated  due  to  the  fact  that  the  velocity 

of  toe  electromagnetic  wave  is  much  faster  than  that  of  the  elastic  waves  and  the  composite 
plate  is  clamped  in  the  Xj-direction  (S3  =  0).  A,  and  C  are  the  coefficients  which  will  be 

determined  by  toe  boundary  conditions  which  are  traction  fiee  and  d>=±V/2  (here  exp(-j(0t)  is 
omitted)  at  X3=d  t/2. 

The  stresses  T3,  Tj  and  electric  displacement  D3  are  expressed  as: 

m  3 

^  =  SS^<"’»)cos(/i^^,)cos(^,X3)4,  +  (C33  +  7^  cos(h^^x^))C  ( 20a ) 

#1=1  f*l 

Tf  =  1  iri3(/i,i)sin(/iSx,)sin()3„X3)A„  (20b) 

n=li=l 

m  3 

=  ZZ  (/i,i)cos(ft^x, )cos()3„X3)A,  +  (-4  + cos(4.x, ))C  ( 20c ) 

#1=1  f=l 

T!  =  XZ^33(”>0cos(4(x,  -^))cos(P„X3)4  +  cos(/jo';(x,  -4))C 

n=l  1=1  ^  2 


(20d) 


(20e) 


Ts  =  1  lTi^(n,i)smih^ixi-^))sin{P„X2)A„ 

n=li~l  2 

^3'’  =  XS ^:r («.t)cos(A^(x,  -  ^))cos0„Xj)\  - £^^c  ( 20f ) 

R«1  is|  2 

Where  7^3 (n,{),  7j3(n,i),  Df{n,i),  T^{n,i),  T^j(n^i)  and  Dj(n,i)  arc  the  sanw*  as 

^33(0,  73?(0.  A^(0.  233(0.  2'3^(0.  and  D3'’(i)  in  equations  (11)  to  (14)  if  the  subscript!  there 
is  replaced  by  n  and  i  for  i-th  partial  wave  of  n-th  mode  here.  And 

2?^  =  cf3A,X.  =  «3,^W- 

The  number  of  die  eigenfunctions,  m,  required  in  die  expansion  is  determined  by  the 
accuracy  needed  for  die  solution.  For  the  problem  treated  here,  we  found  that  it  is  aH/./piatA  to 

use  eight  eigen-modes  in  the  expansion.  In  the  frequency  range  smdied  (  fd<2  MHz*mm), 
there  are  two  branches  having  real  P  and  the  other  branches  having  either  imaginaiy  or  complex 

p,  which  corresponds  to  the  modes  confmed  at  the  surface  of  Xy=±X/2. 

For  the  problem  treated  here,  all  the  stress  components  in  air  are  zero  and  for  the  galr^  of 
simplicity,  D  in  air  is  also  assumed  to  be  zero  since  the  dielectric  permittivity  of  the  con^site 
is  much  higher  than  air.  Under  these  conditions,  the  variational  formula  takes  the  following 
form: 

where  the  integration  over  time  has  been  performed  to  take  into  account  of  the  complex 
notations  for  the  quantities  in  the  integrands.*’-^'  The  *  represents  the  complex  conjugate  of  the 

corresponding  quantity.  The  first  integral  is  over  the  volume  of  the  2-2  composite  plate  and  it 
can  be  shown  that  it  is  equal  to  zero  since  all  the  quantities  in  the  integrand  satisfy  eq.  (1).  The 

second  integral  is  over  the  surfaces  of  the  composite  plate  at  x,=±t  H  where  O  =±V/2  at  x^=±xJ2 

where  V  is  the  applied  voltage.  The  periodic  condition  of  the  composite  in  the  x, -direction  and 
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the  symmetric  condition  of  the  solutions  and  the  boundary  conditions  about  the  plane  of  XjS  0 
allow  the  second  integration  to  be  performed  over  one  unit  cell  at  Xj  =  1 12. 

Substituting  Tj,  T,,  UpUj,  <I>  and  D,  into  equation  (21)  yields  the  following  linear  algebraic 

equation: 

<“) 

where  (M^)  is  a  9x9  matrix,  (Aj)=(A,,A2,Aj,A4,A5,Aj,A7,A,,C)^,  (Vj)  is  a  9x1  matrix  whose 

elements  depend  only  on  V,  the  {q>plied  voltage.  For  a  given  frequency  f,  (mjc  can  solve  eq. 
(22)  to  obtain  Aj  and  C.  From  Aj  and  C,  eqs.  (18),  (19),  and  (20)  yield  all  the  diaracteristic 
properties  related  to  the  vibration  of  a  2-2  conq>osite  plate  with  different  thickness  and  different 
ceramic  volume  fraction  such  as  the  electrical  impedance,  surface  displacement  distribution, 
resonant  modes,  electromechanical  coupling  coefficient,  etc. 

For  the  composite  plate  treated  here,  die  electrical  impedance  for  a  single  repeating  unit  can 
be  found  from: 

V 

Z  =  y  (23) 

where  I  is  the  current  which  is  equal  to  I  =  —  = -jcoQ  = -jlcob  \  D^dx^,  where  Qand  b  are 

dt  i 

the  electric  charge  and  the  length  in  the  Xj-diiection  of  the  plate  respectively.  To  cnmpam  with 
the  experimental  result,  the  current  I  should  be  multiplied  by  N,  the  number  of  rq)eating  unit  in 
a  composite  plate. 

Shown  in  figure  6(a)  is  the  electric  impedance  curve  calculated  from  eqs.  (22)  and  (23)  for  a 
composite  plate  with  t/d=4.5  and  the  ceramic  volume  fraction  v=44%.  The  electric  impedance 

measured  experimentally  from  the  same  composite  plate  is  shown  in  figure  6(b)  and  clearly  the 
theoretical  impedance  curve  reproduces  the  experimental  data  quite  well.  The  difference  in  the 
sharpness  of  the  resonant  peaks  between  the  experimental  data  and  the  theoretical  curve  is  due 
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to  the  fact  that  in  the  theoretical  analysis,  the  electrical  and  mechanical  losses  of  the  ceramic  and 
pol)rmer  phases  were  not  included. 

In  figure  6(c),  the  peak  positions  from  the  experimental  data  are  compared  with  the 
dispersion  curves  for  this  composite  which  shows  excellent  agreement  between  the  two.  In 
order  to  elucidate  the  nature  of  these  modes,  the  spatial  distribution  of  u,  at  each  mode  at  the 
composite  surface  is  jnesented  in  figure  7  which  are  evaluated  based  on  eqs.  (18)  and  (19). 
Apparently,  4,  is  the  fundamental  thickness  resonance  and  f„  is  the  first  lateral  mode  as 

revealed  by  the  fact  that  the  ceramic  and  polymer  vibrate  180°  out  of  phase  at  this  mode,  as 

predicted  in  the  earlier  theoretical  work.^**  The  frequency  position  and  the  distribution  of  u, 
along  the  Xj  axis  indicate  that  fjj  is  the  third  harmonic  of  the  thickness  mode.  However,  die 
appearance  of  f^  is  not  expected  from  the  earlier  theoretical  works  in  which  the  rt^ramir  and 

polymer  vibrate  in  phase.  By  examining  the  equations  of  the  boundary  conditions  at  X3=±t/2,  it 
can  be  deduced  that  a  resorumce  will  occur  whenever  P  =  (l+2n)jc/t,  i.e.,  cos  (Pt/2)  =  0. 
Hence,  the  dispersion  curvw-of  real  p,  as  shown  in  figure  6(c),  reveal  that  the  fimHamAntai 
thickness  resonance  and  the  first  lateral  resonance  occur  at  P  =  7c/t  (4,  and  ^,>.  Similarly,  when 

P  =  3  jc^,  the  third  harmonic  of  the  thickness  mode  will  occur  at  f^,.  In  addition,  a  mode  f^  will 

also  show  up  at  the  branch  1  which  is  at  a  frequency  near  and  above  f„.  By  the  same  argument, 
it  would  be  expected  that  fj^,  f,j,  etc.  may  also  be  observed,  depending  on  the 
electromechanical  coupling  factors  of  these  modes.  As  shown  later,  the  effective  coupling 

factor  for  the  modes  in  the  first  branch  will  decrease  with  increasing  Pd,  i.e.,  reducing 
thickness,  and  on  the  other  hand,  the  effective  coupling  factor  for  the  modes  in  the  second 
branch  will  increase  with  p.  As  the  ratio  of  t/d  decreases,  the  frequency  of  the  thickness  mode 

will  gradually  move  towards  the  first  lateral  mode  which  will  become  stronger  (coupling  factor 
increases),  and  the  second  lateral  mode  (fjj)  becomes  weaker  and  finally  it  will  disappear. 
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experimental  observations  on  how  various  modes  change  with  temperature  (which  causes 
reduction  of  the  shear  velocity  of  the  polymer  phase)  and  the  composite  thickness,^ 

Both  the  experimental  results  and  the  theoretical  data  indicate  that  the  ceramic  and  polymer 
vibrate  in  phase  for  all  the  modes  on  the  first  branch  and  out  of  phase  for  modes  on  the  second 
branch. 

The  electromechanical  coupling  factor  for  the  thickness  resonance  can  be  evaluated  based  on 
the  definition  of  IEEE:” 

(24) 

24  2  4 

where  k,  is  the  thickness  mode  coupling  factor,  f,  and  fp  are  the  series  and  parallel  resonance 
frequencies,  respectively.  Eq.  (24)  is  used  here  to  calculate  the  coupling  factor  for  the  modes  in 
both  the  first  and  second  branches.  Shown  in  figure  9  is  the  results  for  a  2-2  composite  plate 
with  44%  ceranaic  volume  content  for  different  d/t,  where  both  theoretical  and  experimental 
results  are  presented.  As  the  ratio  of  d/t  increases,  the  coupling  factor  of  the  mode  (thickness 
mode)  in  the  first  branch  gradually  decreases,  while  the  coupling  factor  of  the  mode  in  die 
second  branch  gradually  increases  due  to  the  modes  coupling.  As  d/t  is  further  increased,  die 
thicknftS'g  mode  will  jump  to  the  second  branch  when  the  coupling  factor  in  the  second  branch 
surpasses  that  of  the  first  branch.  Although  in  this  region  the  coupling  factor  for  the  thickness 
mode  can  still  be  quite  high,  the  distribution  of  the  Uj  is  not  uniform  on  the  composite  surface 

the  ^ftramir.  anH  polymer  vibrate  180®  out  of  phase  which  is  not  desirable  since  the  polymer 

phase  will  not  be  able  to  perform  properly  the  function  of  transfering  the  acoustic  energy 
between  the  ceramic  plates  and  the  external  medium. 

For  a  composite  plate  to  work  effectively  as  an  electromechanical  transduction  material,  it  is 
required  that  the  ceramic  and  the  polymer  plates  in  the  composite  vibrate  in  phase  with  nearly 
the  same  amplitude  in  the  Xj-direction.  The  evolution  of  the  vibration  pattern  in  the  two  phases 
with  frequency  and  the  effect  of  the  aspect  ratio  t/d  of  a  composite  plate  on  this  distribution  are 


19 


investigated.  Shown  in  figure  10(a)  is  the  variation  of  the  ratio  u^lu^  at  the  surface  of  the 

composite  plate,  where  u,'’ and  ii|  are  Ujat  the  centers  of  the  polymer  (x,ad/2)  and  the  ceramic 
plates  (x,=0),  respectively,  with  frequency  for  the  composite  plate  of  t/d=4  where  both  the 
experimental  results  and  theoretically  calculated  curve  are  presented.  The  experimental  data 
were  acquired  using  a  laser  dilatometer."  Hence,  at  frequencies  far  below  any  resonant  mode, 

is  always  less  than  one,  which  is  true  as  long  as  the  composite  is  driven  electrically. 
For  a  given  frequency,  as  t/id  increases,  this  ratio  increases  and  rqrproaches  one.  These  are 
consistent  with  the  results  of  the  earlier  theoretical  work  on  the  static  properties  of 

conqx)sites.“  As  frequency  increases,  the  ratio  increases  towards  one.  At  a  fiequency 
f,  which  is  near  f,  of  the  thickness  mode,  itf  1.  This  ratio  will  surpass  one  as  the 
fiequency  is  further  increased.  This  is  true  as  long  as  fL,<f„.  In  figure  10(b),  the  change  of 

f,/f,  vs.  the  ratio  of  d/t  is  presented.  Qearly,  f,/f,  is  close  to  but  larger  than  one  except  for 
con^ite  plates  with  a  small  aspect  ratio.  Hence,  Ae  aspect  ratio  t/d  does  not  have  a 
significant  effect  on  the  ratio  of  lu^  at  fiequencies  very  near  f,  of  the  thickness  mode,  where 
/of  is  always,  near  one  if  the  thickness  mode  is  below  the  first  iateral  mode  fiequency. 
However,  it  wUl  affect  the  bandwidth  in  which  /u^  is  near  one.  For  example,  the  bandwidth 
Afif,,  where  Af  is  defined  here  as  the  fiequency  width  in  which  0.9  <af /uf  <1.1,  increases  as 

the  ratio  of  d/t  decreases,  which  is  shown  in  figure  10(b).  The  experimental  data  points  are  also 
presented  in  figure  10(b)  and  the  agreement  between  the  two  is  quite  good. 

V.  Forced  Vibration  of  a  Piezocomposite  Plate  in  a  Fluid  Medium 

Following  a  similar  procedure  as  outlined  in  the  preceding  section,  the  vibration  problem  of 
a  composite  plate  in  a  fluid  medium  can  also  be  analyzed  which  are  more  relevant  to  the 
practical  design  and  application  of  a  composite  transducer.  The  fluid  medium  chosen  for  the 
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study  is  water  and  two  situations  will  be  investigated:  the  con^site  plate  as  a  transmitter,  i.e., 
under  a  harmonic  electric  field,  and  the  composite  as  a  receiver,  i.e.,  under  a  hanoonic 
acoustic  pressure  in  water. 

(5.1)  Forced  Vibration  of  a  Composite  Plate  in  Water  under  an  AC  Electric  Field 

For  the  sake  of  convenience,  we  will  treat  the  system  as  a  symmetric  one  in  which  the 
composite  plate  is  loaded  by  water  on  both  surfaces,  hence,  the  boundary  conditions  at  the 
composite-water  interface  are: 

7^=0,  and  =  atX3=±-  (25) 

where  superscripts  Cm  and  W  denote  the  quantities  in  the  conqx)site  and  water,  and  the  factor 
e~^“‘  is  omitted  in  electric  potential  4>.  The  variational  formula  for  this  problem  can  be  derived 

)dS  =  o  (26) 

where  S*  indicates  that  the  integration  is  over  the  two  composite-water  interfaces  (at  Xj*  ±  t/2), 

and  the  integration  over  the  volume  has  been  omitted  since  it  is  equal  to  zero. 

Here,  the  expressions  for  the  elastic  displacement  vectors  and  the  electrical  potential  in  the 

composite  are  those  in  eq.  (18)  and  (19).  Because  of  the  periodic  nature  of  the  composite  plgt** 
in  the  x, -direction,  the  solutions  in  water  have  the  form: 

=;■  sin(/j^  Xi')e\p(±jP^  x2)R^ 

n=0 

cosih^ xi  )exp(±y^^ X3  )R^  ( 27  ) 

n=0 

Icos(/^ xi)exp(,±jp^ X2)R^ 

11=0 
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where  s  ^  _ / ii^  \2  *pw  f  R^  ^  i/*^  ^  i_ 

wnerc  «.  -  ^  ,  p,  -c  (^0  )  .  Po  =-^.  ^  ^  J~w  “  *® 


longitudinal  wave  velocity  of  water,  c*'  is  the  bulk  modulus  of  water.  ±  correspond  to  Ae 


solutions  in  X3  >0  region  and  Xj  <  0  region,  respectively,  h  and  P  are  the  wave  vector 

components  in  the  x,-  and  Xj-directions  in  water.  is  determined  by  the  boundary  c<mditions 
through  the  variational  principle  formula  eq.  (26). 

Substituting  equations  (18),  (19)  and  (27)  into  equation  (26)  yields  a  set  of  linear  algetaaic 
equadcHis: 

(Mij)(Aj)=(Vi)  (28) 

where  (M^p  is  a  matrix  in  which  its  elements  are  related  to  the  parameters  of  water,  ceramic  and 
polymer  and  the  geometrical  parameters  v  and  d,  as  well  as  m,  P  and  h.  (Aj)=(A,,  ....  A„,  C, 
w  w 

Rq  , ...,  Rj  ),  where  m  and  J  are  the  numbers  of  the  eigenfunctions  used  in  the  expansions 
for  the  quantities  in  the  compoate  and  in  water,  respectively,  and  how  many  eigenfunctions 
should  be  used  in  the  expansion  depends  on  the  accuracy  desired.  In  this  calculation,  m=8  and 

J=6.  Hence,  (Vp  is  a  16x1  matrix  which  elements  depend  on  the  applied  voltage  V.  Equation 
(28)  is  solved  numerically. 

Shown  in  figure  1 1(a)  is  a  theoretical  electrical  impedance  curve  for  a  2-2  composite  plate  of 
44%  ceramic  volume  content  and  the  aspect  ratio  t/d  =  3.8  loaded  with  water.  For  the 
comparison,  the  experimental  curve  for  the  same  composite  plate  is  shoun  in  figure  1 1(b).  The 
agreement  between  the  experimental  result  and  the  theoretical  one  is  quite  good.  The  relatively 
sharp  resonant  peaks  in  the  theoretical  curve  compared  with  the  experimental  one  are  due  to  the 
fact  that  in  the  theoretical  calculation,  the  elastic  and  dielectric  losses  of  the  composite  plate  ate 
not  included.  Comparison  betueen  figure  1 1  with  figure  6,  which  is  the  electrical  impedance 
for  a  similar  composite  plate  in  air,  reveals  that  the  resonance  is  severely  Hamp^H  in  water  as 
shown  by  the  marked  broadening  of  the  resonant  peaks  in  the  impedance  curve.  For  the  lateral 
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modes  (modes  1  and  2  in  figure  II)  which  is  mainly  determined  by  the  shear  resonance  of  the 
polymer  phase  in  the  2-2  composite  investigated,  because  the  acoustic  impedance  between  the 
polymer  and  water  is  very  close,  the  change  in  the  amplitude  is  quite  significant.  One 
interesting  feature  of  the  influence  of  the  water  loading  on  a  composite  plate  is  the  appearance 
of  a  mode  which  is  labeled  3  in  the  figure.  It  can  be  shown  that  the  resonant  frequency  for  this 
mode  does  not  change  very  much  as  the  thickness  of  the  composite  plate  changes.  From  the 
stress  distribution  pattern  in  the  water,  it  is  not  difficult  to  show  that  this  mode  is  related  to  the 
coupling  between  the  composite  and  water  and  its  fiequency  is  determined  by  the  periodicity  d 

and  the  acoustic  wave  velocity  of  water  (f  *  V/d).  And  it  corresponds  to  the  local  oscillatioD 
of  water  within  one  unit  cell. 

In  analogous  to  the  situation  in  air,  it  is  also  found  that  in  water,  at  frequencies  far  below 
any  resonant  mode,  the  ratio  of  /u^  is  always  less  than  one  and  wiU  approach  one  as  the 
aspect  ratio  t/d  increases.  As  the  frequency  increases  towards  the  thickness  resonant 
frequency,  u^/i^  increases  towards  one.  At  a  fiequency  f,  near  the  thickness  resonant 
fiequency  f,,  this  ratio  becomes  orre  ^d  abpve  that,  ^  ratio  is  larger  than  one  (figure  12(a)). 
Figure  12(b)  presents  f,/f,  and  ^f^  ys.  the  ratio  of  dA  foj  a  coiD^site  wi^^  cerarnk 

content.  Clearly,  the  effect  of  the  aspect  ratio  t/d  on  the  surface  uniformity  of  a  conqxtsite  piatp 
is  to  change  the  frequency  position  f,  with  respect  to  the  thickness  resonance  and  the  frequency 

width  Af  in  which  /u^  is  near  one  which  are  very  similar  to  that  found  in  air.  However,  the 

bandwidth  Af/f,  in  water  is  larger  than  that  in  air,  indicating  that  the  surface  distribution  of  the 

displacement  of  the  two  phases  is  much  flatter  in  water  than  in  air  which  is  quite 
understandable.  Because  of  the  water  loading,  the  vibration  amplitude  of  the  polymer  phase  is 
significantiy  reduced.  Figure  12(a)  shows  how  the  ratio  of  Uj^/Uj'  at  the  composite  surface 
varies  with  frequency  (t/d=4).  One  noticeable  change  between  the  surface  profiles  in  air  and  in 
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water  is  that  in  water,  even  at  the  lateral  mode  frequency,  this  ratio  does  not  become  very  large 
(not  shown  in  the  figure  since  the  lateral  mode  frequency  is  at  about  1.5  MHz). 

(5.2)  Forced  Vibration  of  a  Composite  Plate  in  Water  under  Harmonic  Acoustic  Pressure 
Now  we  turn  to  investigate  die  vibration  behavior  of  a  composite  plate  in  water  under 
harmonic  acoustic  pressure,  that  is,  to  study  issues  related  to  the  receiving  sensitivity  of  a 
sm^le  2-2  composite  transducer.  For  simplicity,  we  shall  restrict  the  treatment  to  die 
symmetric  system.  In  this  case,  as  sketched  in  figure  13,  the  incident  acoustic  waves  impinge 

normally  on  the  two  surfaces  of  the  conqwsite  plate  from  opposite  directions  and  the  boundary 
conditions  for  this  problem  are; 

=  T^.  Jf "  =  0,  U3°”  =  oj""  md  Df”  =  0  at  X3  =  ±1  (29) 

The  appropriate  variational  formula  is 

-  jf”Ki(u3"')*  +«(a3^”)')+(«3C"  -«3"')(i(r3"')*  +5(r3°»)*)jB 

=  0  '  (30) 

The  expressions  of  the  displacements  and  electrical  potential  in  the  composite  plate  are  those  in 

equations  (18)  and  (19)  except  now  the  mode  of  P=0  should  not  be  included  since  there  is  no 

externally  applied  electric  field.  To  account  fw  the  incident  wave  in  water,  equation  (27)  is 
modified  and  the  solutions  in  the  Xj  >  0  region  are: 

«r  =7  iCsin(h^Xi),xpijP^x,)lC 

n=0 

=  -Pq  exp(-j^0^3 )+  cos(h^  Xi)txp{jp^  (31) 

n=0 

^3*^  =yTo''[exp(-;^0'^3)+  Icos(/i^x,)exp(y^^j:3)/?,J''  ] 

n-0 


24 


where  the  term  of  -Pq  ^xpi—jPQX^)  corresponds  to  the  incident  plane  wave.  By  substituting 
equations  (18),  (19)  and  (31)  into  (30),  we  obtain  a  set  of  linear  algebraic  equations: 

(Mij)(Aj)=(Vi)  (32) 

where  (Mjj)  is  a  15x15  matrix,  (Aj)=(A,,...,  A„,  ...,  r'^)  with  m=8  and  J=6,  and  (V)  is 

a  15x1  matrix  related  to  the  incident  wave.  By  solving  equation  (32)  one  can  obtain  aU  of  the 

properties  related  to  the  behaviors  of  a  composite  plate  under  a  harmonic  acoustic  (oessure. 

Shown  in  figure  14  is  tfte  open  circuit  voltage  receiving  sensitivity  V/p,  where  V  is  the 
voltage  ou^ut  and  p  is  the  pressure  of  the  incident  wave,  fw  a  44%  2-2  conqx>site  plate  at  die 
ttiickness  resonance  vs.  d/t  for  the  thickness  resonance  mode  where  the  cross  over  region 
(point  B  in  figure  2)  is  at  d/t  =  0.65.“  That  is,  at  dA  <  0.65,  the  thickness  mode  is  in  the  first 
branch  of  the  dispersion  curves  and  above  that,  the  thickness  resonance  is  at  the  second  branch 
of  the  dispersion  curves  where  the  polymer  and  ceramic  vibrate  out  of  phase.  Since  die 
problems  treated  in  this  paper  are  related  only  to  the  piezo-materials,  the  open-circuit  receiving 
sensitivity  is  used  here  which,  as  pointed  cot  by  Kojima,  is  a  system  independeiit  parameter.^ 
The  bandwidth  is  defined  ^  die  3dB', width  about  tte  peak'^fie^hcy  tfp  for  the  receiving 
mode).  For  a  single  ph^  mktefiki;  if  (iah  be  d^nv^  'fifom^  KLM  model  that  the  sensitivity 
here  should  be  proportional  to  t,  the  thickness  of  the  transducer.  Hence,  in  figure  14,  the 
vs.  d/t  curve  for  a  single  phase  mat^al  should  fall  off  as  t  (solid  line  in  the  figure)  while  the 
bandwidth  should  stay  constant.  The  results  show  that  the  sensitivity  of  the  thickness  mode 
for  a  2-2  composite  decreases  slowly  as  d/t  becomes  larger  than  0.4,  but  the  bandwidth 
increases  graduaUy.  The  increase  in  the  bandwidth  is  due  to  the  merger  of  the  two  resonant 
modes  in  water.  After  that,  there  are  anomalous  changes  in  both  the  bandwidth  and  the 
sensitivity  in  the  cross-over  region.  At  higher  values  of  dA  (thin  composite  plates),  both  the 
sensitivity  and  the  bandwidth  fall  much  below  the  values  of  single  phase  material  (solid  line). 
If  we  define  the  figure  of  merit  (FOM)  here  as  the  product  of  sensitivity  and  bandwidth,  as 
shown  in  figure  14(b),  at  dA  <  0.5,  the  FOM  falls  off  with  t.  At  dA  >  0.5,  die  FOM  drops  to 
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much  smaller  value.  Therefore,  in  order  to  gain  a  high  receiving  sensitivity  and  a  broad 
bandwidth  of  a  2-2  composite  transducer,  it  is  desirable  to  have  d/t  less  than  0.5.  The  results 
here  can  be  compared  with  what  shown  in  figures  2,  8,  9,  and  12.  At  dA  above  0.5,  the 
coupling  factor  shows  a  precipitous  drop  (figure  9)  and  the  thickness  mode  frequency  also 
shows  an  apparent  deviation  from  the  extrapolated  value  (the  dashed  line  in  figure  2).  Figure  8 
also  shows  that  there  is  a  significant  coupling  between  the  two  modes.  Although  the  results 
presented  are  for  a  composite  with  44%  ceramic  content,  it  is  approximately  jnie  for 
composites  with  other  volume  fractions. 

VI.  Summary 

The  details  of  a  theoretical  model  on  pi^oceramic  polymer  composites  with  laminar 
poiodic  structure  are  presented.  The  result  shows  tfiat  the  various  resonant  modes  in  a 
composite  structure  can  be  traced  back  to  the  modes  in  either  an  isolated  ceramic  plate  or 
polymer  plate  with  ^[ipropriate  boundary  conditions  (stress  free  for  the  ceramic  plam  and  strain 
fiec  for  the  polyiner  pl^).  It  also  shows  that  there  east  a  series  of  modes  associated  with  the 
periodic  structure  of  a  composite,  which  is  beyoi^  the  stop-band  edge  resonance  prediction.  ‘ 
One  of  the  main  concerns  in  designing  a  con^site  transducer  is  how  the  surface  vibration 
profile  changes  with  fiequency  and  how  this  is  influenced  by  the  aspect  ratio  t/d.  It  was 
predicted  and  verified  by  experiment  that  as  long  as  the  thickness  resonance  is  below  the  first 
lateral  mode  fiequency,  there  is  always  a  fiequency  f,  which  is  near  the  thickness  resonance 
and  at  which  the  polymer  and  ceramic  vibrate  in  unison.  The  effect  of  t/d  is  to  change  the 
position  of  f,  with  respect  to  the  thickness  resonance  fiequency  and  the  bandwidth  in  which 
polymer  and  ceramic  have  nearly  the  same  vibration  amplitude  and  phase.  It  is  also  predicted 
that  when  operated  in  a  fluid  medium  such  as  water,  there  wiU  be  a  resonance  mode  whose 
frequency  is  determined  by  the  velocitj-  of  the  fluid  medium  and  the  unit  cell  length  d  and  is 
associated  with  the  oscillation  of  the  fluid,  causing  the  polymer  and  ceramic  to  vibrate  180°  out 
of  phase.  The  difference  in  the  surface  vibration  profiles  between  in  air  and  in  water  indicates 
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the  need  to  characterize  the  vibration  pattern  of  a  composite  in  a  fluid  medium  since  it  is  much 
closer  to  the  real  {plication  environment 

In  general,  the  maximum  transmitting  voltage  sensitivity  of  a  transducer  is  at  a  hequency 
near  f,  and  the  maximum  open  circuit  receiving  sensitivity  is  near  f„  of  the  thickness  mode  (f  > 
f,),  and  hence  it  is  expected  that  the  influence  of  the  lateral  naode  (hence,  the  aspect  ratio  t/d) 
will  be  more  severe  on  the  receiving  sensitivity.  From  the  data  analysis,  it  is  shown  that  when 
an  FOM  which  is  the  product  of  the  sensitivity  and  the  bandwidth  is  introduced  as  the  criterion 
of  the  performance  of  a  composite  transducer  as  a  receiver,  the  performance  deteri<»ates  when 
d/t  >  0.5  for  the  ronqxasite  discussed  here.  Experimental  results  confirm  this  where  as 
shown  in  figure  15,  the  insertion  loss  of  a  2-2  piezocomposite  shows  large  increases  at  d/t  near 
and  above  0.5. 
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Table  I.  The  material  properties  of  PZT-5H  and  Spurr  epoxy  used  for  the  2-2  composites  in 
the  investigation. 


PZT-5H:  ejj  =  23.09  am\  e,,  =  -  6.603  am^  e,,  *  17.0  am^  c„  =  12.72*10'®  N/m^ 

=  2.3*10'®  N/m^  Cjj  =  11.74*10'®  N/m^  c,,  =  8.47*10'®  N/m^  K„  =  1700,  K33  =  1470,  p  = 
7500  kg/m\ 

Spurr  epoxy:  c„  =  5.4*10*  N/m\  C44  *  1.3*10*  N/m^  p  =  1 100  kg/m’. 
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Figure  captions: 

Figure  1.  (a)  Schematic  drawing  of  a  2-2  piezoceramic  polymer  composite  which  is 
unbounded  in  all  three  orthogonal  directions.  The  period  of  the  composite  is  d.  (b)  .^rhpinaHr 
drawing  of  a  2-2  piezoceramic  polymer  composite  plate  with  a  thickness  t  situated  in  air.  An 
external  voltage  is  applied  to  the  composite. 

Figure  2.  The  first  and  second  branches  of  the  dispersion  curves  (solid  lines)  and  the 
experimentally  measured  thickness  mode  (solid  circles)  and  lateral  mode  (open  circles)  for  a  2- 
2  composite  made  of  PZT-5H  and  Spun  epoxy  with  44%  ceramic  content.  After  the  crossover 
region,  the  thickness  mode  jumps  to  the  second  branch.  The  criterion  for  the  thickness  mode  is 
that  it  is  the  one  with  the  larger  effective  coupling  factor. 

Figure  3.  The  dispersion  curves  for  a  2-2  conqwsite  (figure  l(a))made  of  PZT-5H  and 
Spurr  epoxy  for  (a)  15%  ceramic  volume  content  and  (b)  44%  ceramic  volume  content 

Figure  4.  The  effect  of  the  properties  of  the  polymer  matrix  on  the  dispersion  curves  of  2-2 
composites  with  PZT-5H  ceramic,  derived  firom  the  model.  The  properties  of  the  polymer 

phase  are:  curve  1:  c,,*  7.72*10’  N/mVc^=  1.588*10’  N/m^  p  =  1 160  kg/m’;  curve  2:“ 
c„=  1.3634*10’°  N/m’,  944=  3.432* l()’  N/m’,  p  =  1610  curve  3:  c,,=  3.173*10’ 
N/m’,  C44=  0.696*10’  N/m’,  p  =  1060  kg/m’;  curve  4:  c„=  1.622*10’  N/m’,  C44=  1.646*10* 
N/m’,  p  =  890  kg/m’. 

Figure  5.  Comparison  of  the  dispersion  curves  in  2-2  composites  with  those  in  a  single 
piezoceramic  plate  and  a  polymer  plate  for  composites  with  (a)  15%  ceramic  content;  (b)  44% 
ceramic;  and  (c)  80%  ceramic  content.  PZT-5H  is  used  as  the  piezoceramic  and  Spurr  epoxy 
as  the  polymer  phase. 

Figure  6.  The  electric  impedance  cur\'e  for  a  2-2  composite  plate  with  PZT-5H  and  Spurr 
epoxy  (44%  ceramic)  measured  in  air  (a)  the  theoretical  curve  and  (b)  the  experimental  curve. 
The  thickness  t  of  the  composite  is  t/d  =  4.5  and  d=0.635  mm.  fL,  and  fLj  are  the  fundamental 
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and  third  harmonic  of  the  thickness  mode  and  tl  and  t2  are  the  modes  arose  from  the 
periodicity  of  the  conqwsite.  (c)  The  dispersion  curves  elucidating  die  origin  of  the  resonant 

modes  f„  and  f^.  In  general,  any  modes  on  the  dispersion  curves  will  show  up  whenever  ^d  = 
(2n+l)/2,  n=0,l,2, ...  is  satisfied,  where  P=2jt/Jl  and  X  =2t. 

Figure  7.  The  distribution  of  the  elastic  displacement  u,  for  f^,,  f  f  „,  and  f  ^  of  a  2-2 
composite  with  t/da4,  where  (a)  and  (b)  are  the  distributions  along  the  x  3-direction  when  x  ,=0 
(at  the  center  line  of  the  ceramic  plate,  figure  1(b)),  Xj^O  is  at  the  center  and  X3=2  is  at  the 
ceramic  surface;  (c)  and  (d)  arc  surface  vibration  profiles  at  X3=2,  x,=0  corresponds  to  the 
center  of  the  ceramic  plate  and  x  ,=0.5  is  at  the  center  of  the  polymer  plate. 

Figure  8.  The  evolution  of  various  modes  in  a  2-2  composite  made  of  PZT-5H  (44% 
ceramic)  and  Spurr  epoxy  with  the  thickness  t  of  composite  plate:  (a)  t/d  =  4,  (b)  t/d  =  2,  and 
(c)  t/d  =  1.  There  is  already  substantial  coupling  between  the  thickness  and  lateral  modes  at  t/d 
—  2  for  this  composite.  At  low  t/d  values  (thin  samples),  the  lateral  mode  will  disappear. 

Bgure  9.  'Die  eyolution  <rf  the  coupling  factor  for  the  modes  on  the  first  and  second 
branches  of  the  dispersion  curves  with  the  composite  thickness  d/t  for  a  composite  maA.  of 
PZr-5H  and  Spurr  epoxy  with  44%  ceramic  content 

Figure  10.  (a)  The  ratio  of  vs.  frequency  for  a  composite  made  of  PZT-SH  and 
Spurr  epoxy  (44%  ceramic  content)  at  a  thickness  of  t/d  =  4  (d=0.635  mm)  measured  in  air. 
U3P  and  U3*  are  the  surface  displacements  at  the  centers  of  the  polymer  plate  and  ceramic  plate, 
respectively.  The  black  dots  are  the  experimental  data  measured  using  a  laser  dilatometer  and 
the  solid  line  is  derived  from  the  model. 

(b)  The  ratio  f,/fj  as  a  function  of  the  ratio  d/t  of  the  composite  plate,  where  f,  is  the 
frequency  at  which  U3'’/U3'=l  and  fj  is  the  series  resonant  frequency  for  the  thickness  mode. 
The  black  dots  are  experimental  data  and  the  solid  lines  are  derived  from  the  model  calculation. 

(c)  The  bandwidth  as  a  function  of  the  ratio  d/t  of  the  composite  plate.  The  black  dots  are 
experimental  data  and  the  solid  lines  are  derived  from  the  model  calculation. 
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Figure  11.  The  electrical  impedance  magnitude  for  a  2-2  conqwsite  of  PZT-5H  and  Spurr 
epoxy  with  44%  ceramic  content  measured  in  water.  The  thickness  of  the  composite  is  t/d=3.8 
and  d=0.635  mm.  (a)  is  derived  from  the  model  and  (b)  is  the  experimental  results.  The 
modes  1  and  2  are  the  lateral  modes  (f„  and  in  figure  6(c))  and  mode  3  arises  fium  the 
coupling  of  water  to  the  periodic  structure  of  the  compnsitp.  surface. 

Figure  12.  (a)  The  ratio  of  Uj^Uj*  vs.  frequency  for  a  composite  made  of  PZT-5H  and 
Spun-  epoxy  (44%  ceramic  content)  at  a  thickness  of  t/d  =  4  (d=0.635  mm)  measured  in  water. 
^3^  srtd  Uj  are  the  surface  displacements  at  the  centers  of  the  pol3nDer  plate  and  ceramic  plate, 
respectively.  The  results  are  derived  from  the  model 

(b)  The  ratio  f,/f,  as  a  function  of  d/t  of  the  composite  plate,  where  f,  is  the  frequency  at 
which  U3'’/U3'=l  and  f,  is  the  series  resonant  frequency  for  the  thickness  mode.  The  results  are 
derived  from  the  model. 

(c)  The  bandwidth  as  a  function  of  d/t  of  the  composite  plate.  The  results  here  can  be 
compared  with  those  in  figure  10  and  apparently,  and  the  water  loading  improves  die 
uniformity  of  the  vibration  profile  at  the  conqxisite  surface.  The  results  are  derived  from  the 
model. 

Figure  13.  Schematic  drawing  of  a  2-2  composite  plate  under  normal  incident  of  a  pressure 
wave  from  the  water. 

Figure  14.  (a)  The  influence  of  the  unit  cell  dimension  ratio  of  d/t  on  the  open  circuit 
voltage  sensitivity  of  a  2-2  composite  with  44%  PZT-SH  and  Spurr  epoxy  matrix.  Solid  dots 
are  V/p  at  the  peak  frequency  (fp  and  the  solid  line  is  the  sensitivity  of  a  single  phase  material 
with  the  effective  properties  of  the  composite.  The  results  are  derived  from  the  model. 

(b)  The  3  dB  bandwidth  as  a  function  of  d/t  for  a  composite  transducer  derived  from  the 
model. 

(c)  FOM  (solid  dots,  the  product  of  the  sensitivity  and  the  bandwidth)  for  the  2-2 
composite  in  (a)  as  a  function  of  d/t.  The  solid  line  is  FOM  for  a  single  phase  material  with  the 
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effective  properties  of  the  composite.  At  d/t  >  0.5,  FOM  drops  much  below  the  single  phase 
material  value.  The  results  are  derived  from  the  model. 

Figure  15.  The  experimental  data  on  the  insertion  loss  (open  circles)  and  6  dB  bandwidth 
(solid  dots)  of  2-2  composites  with  44%  volume  content  of  PZT-5H  and  Spurr  epoxy  polymer 
matrix  measured  by  the  pulse  echo  method  for  different  ratio  of  d/t  (d=0.635  mm).  The  solid 
lines  are  drawn  to  guide  eyes. 
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These  features  are  summarized  in  figure  8  which  provide  understanding  on  the  earlier 
experimental  observations  on  how  various  modes  change  with  tenqjeratuie  (which  causes 
reduction  of  the  shear  velocity  of  the  polymer  phase)  and  the  composite  thickness.® 

Both  the  experimental  results  and  the  theoretical  data  indicate  that  the  ceramic  and  polymer 
vibrate  in  phase  for  all  the  modes  on  the  first  branch  and  out  of  phase  for  modes  on  the  second 
branch. 

The  electromechanical  coupling  factor  for  the  thickness  resonance  can  be  evaluated  based  on 
the  definition  of  IEEE:” 

tf  (24) 

^  jp  ^  Jp 

where  k,  is  the  thickness  mode  coupling  factor,  f,  and  f^  are  the  series  and  parallel  resonance 
frequencies,  respectively.  Eq.  (24)  is  used  here  to  calculate  the  coupling  factor  for  the  modes  in 
both  the  first  and  second  branches.  Shown  in  figure  9  is  the  results  for  a  2-2  composite  plate 
with  44%  ceramic  volume  content  for  different  d/t,  where  both  theoretical  and  experimental 
results  are  presented.'  As  the  ratio  of  d/t  increases,  the  coupling  factor  of  the  mode  (thickness 
mode)  in  the  first  branch  gradually  decreases,  while  the  coupling  factor  of  the  mode  in  the 
second  branch  gradually  increases  due  to  the  modes  coupling.  As  d/t  increases  further,  the 
thickness  mode  will  jump  to  the  second  branch  when  the  coupling  factor  in  the  second  branch 
surpasses  that  of  the  first  branch.  Although  in  this  region  the  coupling  factor  for  the  thickness 
mode  can  still  be  quite  high,  the  distribution  of  Uj  is  not  uniform  on  the  composite  surface  and 

the  ceramic  and  polymer  vibrate  180°  out  of  phase  which  is  not  desirable  since  the  polymer 

phase  will  not  be  able  to  perform  properly  the  function  of  transferring  the  acoustic  energy 
between  the  ceramic  plates  and  the  external  medium. 

For  a  composite  plate  to  work  effectively  as  an  electromechanical  transduction  material,  it  is 
required  that  the  ceramic  and  the  polymer  plates  in  the  composite  vibrate  in  phase  with  nearly 
the  same  amplitude  in  the  Xj-direction.  The  evolution  of  the  vibration  pattern  in  the  two  phases 
with  frequency  and  the  effect  of  the  aspect  ratio  t/d  of  a  composite  plate  on  this  distribution  are 
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Abstract: 

The  acoustic  properties  of  a  medium-composite  interface  and  the  field  distributions  in  a  2-2 
composite  are  analyzed  based  on  an  approach  developed  recently.  It  is  shown  that  the 
reflection  coefBcient  from  the  interface  is  a  complex  number  and  the  reflected  wave  suffers  a 
more  than  180®  phase  change.  In  addition  to  that  the  effective  input  acoustic  impedance  Zj„  of 
the  composite  is  a  complex  number,  it  is  found  that  Z„,  also  shows  dependence  on  the  shear 
stiffness  constant  of  the  medium  in  contact,  which  is  related  to  the  non-uniform  vibration 
distribution  at  the  composite  surface.  Since  for  a  piezoceramic  polymer  composite,  it  is  the 
ceramic  phase  which  performs  the  energy  conversion  between  the  acoustic  and  electric  forms, 
the  amount  of  acoustic  energy  which  can  enter  the  ceramic  region  is  one  of  the  most  important 
parameters  in  a  composite  transducer  design.  In  the  paper,  we  show  that  even  though  the 
effective  transmission  coefficient  increases  as  the  frequency  is  increased,  the  amount  of 
acoustic  energy  entering  the  ceramic  region  actually  decreases.  From  the  fact  that  there  is  more 
than  180°  phase  change  in  the  reflection  wave  from  the  medium-composite  interface,  it  is 
shown  that  the  matching  layer  thickness  is  no  longer  equal  to  the  quarter  wavelength,  but 
smaller.  These  predictions  are  confirmed  experimentally. 

PACS  No:  43.35.  -c,  68.35.  Gy,  85.50.Ly 
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I.  Introduction: 


Transmission  and  reflection  of  a  wave  at  an  interface  is  of  great  importance  in  many  areas  of 
modem  technology  such  as  the  optical,  electromagnetic,  acoustic,  and  materials  nondestructive 
evaluations.  For  homogeneous  materials,  a  quantitative  understanding  of  how  the  material 
properties  at  the  interface  influence  the  transmission  and  reflection  of  various  wave  components 
has  been  well  established.’’^  However,  for  materials  with  heterogeneity  such  as  composites, 
much  work  remains  to  be  done  in  order  to  provide  both  qualitative  and  quantitative 
understandings  on  how  a  composite  stmcture  affect  flows  of  waves  at  an  interface  between  a 
composite  and  a  homogeneous  fluid  or  solid  medium. 

For  the  piezoceramic  polymer  composites  used  for  ultrasonic  transducer  applications,^*  in 
addition  to  the  wave  transmission  and  reflection  at  the  interface,  one  has  also  to  deal  with  the 
acoustic  energy  transfer  between  the  two  constituents,  i.e.,  the  piezoceramic  and  polymer 
phases.  In  a  piezoceramic  polymer  composite,  because  the  energy  conversion  can  only  be 
carried  out  in  the  ceramic  phase  and  the  polymer  phase  acts  as  a  carrier  to  transfer  the  acoustic 
energy  between  the  piezoceramic  and  the  external  medium  to  which  the  transducer  is 
interacting.  Clearly,  to  establish  a  stracture  performance  relationship  in  this  type  of  matprial, 
both  the  propagation  of  the  acoustic  energy  at  the  composite-medium  interface  and  the 
effectiveness  of  the  energy  transfer  between  the  polymer  and  ceramic  inside  a  composite  have 
to  be  considered. 

The  objective  of  this  paper  is  to  examine  quantitatively  how  the  material  properties  at  the 
two  sides  of  the  interface,  i.e.,  the  properties  of  a  piezoceramic  polymer  composite  and  a 
uniform  medium,  either  a  liquid  or  a  solid,  influence  the  wave  transmission  and  reflection  at  the 
interface  and  the  acoustic  energy  transfer  inside  the  composite,  and  how  these  affect  the 
parameters  of  the  matching  layer  at  the  interface.  The  approach  taken  here  is  based  on  the 
eigen-mode  expansion  method  and  variational  technique.^’  ® 

The  organization  of  the  paper  is  as  follows:  In  section  11,  the  general  formula  used  to  treat 
the  wave  propagation  at  the  interface  is  briefly  presented  and  the  results  for  the  fluid-composite 
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interface  are  derived.  The  definition  of  the  input  acoustic  impedance  of  the  composite  at  the 
interface  is  also  discussed.  The  acoustic  energy  transfer  between  the  ceramic  and  polymer  is 
analyzed.  In  section  III,  the  wave  propagation  at  a  interface  between  a  solid  and  a 
piezocomposite  is  analyzed.  In  section  IV,  the  issues  related  to  the  anti-reflection  quarter  wave 
matching  layer  at  a  composite-medium  interface  are  examined.  And  in  section  V,  a  brief 
summary  and  discussion  are  presented. 

In  the  paper,  the  piezoceramic  polymer  composites  used  in  the  model  calculations  and  in 
experiments  are  those  made  of  PZT-5H  piezoceramic  and  Spurr  epoxy  polymer  matrix.’  The 
composites  have  a  2-2  connectivity  (a  terminology  introduced  by  Newnham),^  i.e.,  the  ceramic 
and  polymer  form  a  periodic  laminated  structure  as  shown  schematically  in  figure  1(a). 

n.  Plane  acoustic  wave  propagation  at  a  fluid-2-2  composite  interface 

(2.1)  Eigen-modes  in  a  2-2  piezocomposite  and  fluid  medium  and  variational  formula  at 

their  interface 

The  schematic  of  the  interface  system  is  shown  in  figure  1(b)  where  a  plane  acoustic  wave 
is  normally  incident  upon  the  boundary  between  a  fluid  and  a  2-2  piezocomposite.  The 
coordinate  system  is  chosen  such  that  the  Xj-axis  is  perpendicular  to  the  interface,  the  x,-axis  is 
parallel  to  the  interface  as  shown  in  the  figure,  and  the  Xj-axis  is  perpendicular  to  both  the  x, 
and  X3  axes  (pointing  out  of  the  paper).  The  fluid  occupies  the  space  of  Xj  <  0  and  the  2-2 
composite  is  in  Xj  >  0.  Since  the  dimension  in  the  Xj-direction  is  much  larger  than  the  period  d 
of  the  2-2  composite  and  the  wave  fields  are  uniform  in  that  direction,  it  can  be  taken  as  infinite 
without  much  error  in  the  results.  Hence,  for  the  dynamic  problem  to  be  considered,  the  2-2 
composite  can  be  regarded  as  elastically  clamped  in  the  Xj-direction  so  that  $2=  0,  where  is 
the  strain  in  that  direction,  and  we  have  a  two  dimensional  problem  with  no  dependence  on  the 
Xj-variable. 

The  governing  equations  for  this  problem  are 
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Where  i  and  j=l  and  3.  The  symbols  adopted  in  this  paper  are  summarized  here:  Ty  and  Sy  are 
the  stress  and  strain  tensor  components,  where  the  matrix  notation  is  used,  Uj  is  the  elastic 

displacement  vector,  p  is  the  mass  density,  is  the  electric  displacement  vector,  and  Ej  is  the 
electric  field.  The  relevant  material  coefficients  are:  Cy  is  the  piezoelectric  coefficient,  Cy  is  the 
elastic  stiffness,  and  e,  is  the  dielectric  permittivity. 

The  constitutive  equations  in  the  composite  are: 

tr>l=W[S]+[£*][E] 


The  boundary  conditions  at  the  fluid-2-2  composite  interface  for  the  problem  are 


^33  ~  ^33  *  ^13  ~  ^>“3  ~  “3 

^Cm  _  j^m  _  ^ 


atX3=0  (3) 


The  superscripts  Cm  and  W  represent  the  2-2  composite  and  the  fluid.  And  4)  is  the  electric 
potential. 

In  the  2-2  composite,  as  has  been  shown  in  earlier  publications,*  '®  the  displacements  and 
electric  potential  in  the  ceramic  phase  can  be  expressed  as 

u]  =SRfffcos(hfx,)exp(jpX3) 

i 

Ui  =  jl  Ri  gf  Sin(hfx,  )expG  pXj)  (4) 

i 

=IRftfcos(hfx,)exp(jPx3) 


where  i  mns  from  1  to  3.  In  eq.  (4),  hj  and  P  are  the  wave  vector  components  in  the  x,- 


and  Xj-directions,  respectively.  For  each  P,  there  are  three  h,  correspond  to  the  quasi¬ 
longitudinal,  quasi-shear,  and  quasi-electric  waves,  respectively.  For  the  polymer  phase, 
we  have  the  similar  expressions. 


(5) 


•4  =  <^os(hP{xi  -■^))exp(;j3.r3) 

4  =  sin(;i,^(x,  -^))exp(jPx2) 

i  ^ 

=  C^coshipixi  -^))e\pijpxi) 

where  i=l  and  2.  In  equations  (4)  and  (5),  ff ,  g®  ,  t?,  fj*’ ,  and  gf  are  the  coefficients  and 
superscripts  c  and  p  refer  to  the  ceramic  and  polymer  regions.  In  these  equations,  the  time 
dependent  term,  i.e.,  exp(-jtDt),  is  omitted  for  the  sake  of  simplicity. 

By  imposing  the  boundary  conditions  at  the  ceramic-polymer  interface,  one  can  obtain  the 
dispersion  relations  in  a  2-2  composite  for  waves  in  the  Xj-direction.  As  an  example,  the 
dispersion  curves  for  a  2-2  composite  made  of  PZT-5H  piezoceramic  and  Spurr  epoxy  with 
44%  ceramic  volume  content  are  shown  in  figure  2.  In  addition,  the  relations  among 

Rf  ,R2,R3,Rf  ,R2  and  are  determined  and  for  a  given  mode  (a  given  P),  there  is  only  one 

independent  variable  among  these  6  coefficients.  Therefore,  they  can  be  expressed  as 

(Rf,R^,R§,Rf,R5,CP)‘A,,where  A  is  an  undetermined  constant.  Hence,  eqs.  (4)  and  (5) 

constitute  the  eigenfunctions  in  a  2-2  composite  for  a  wave  propagating  in  the  Xj-direction. 

In  the  fluid  medium,  the  constitutive  relation  has  the  form 

COTj  C«3 

where  is  the  bulk  modulus  of  the  fluid.  The  periodic  condition  in  the  x, -direction  due  to  the 

composite  structure  also  imposes  the  constraints  on  the  solutions  to  the  equation  (1)  in  the  fluid 
region: 

ur=jh:sin(h^'x,)exp(jP^X3)R:f^ 

U3'^=p^'cos(h^^x,)exp(jp^3)R„^ 
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where 


andp:=^ 


ft.'M  \2 


(n=0, 1,2,3...) 


are  the  wave  vector 


components  in  the  x,-  and  Xj-directions  in  the  fluid  region, 


is  the  acoustic 


wave  velocity  of  fluid.  In  eq,  (7),  is  an  undetermined  constant. 

In  order  to  satisfy  the  boundary  conditions  at  the  composite-fluid  interface,  the  elastic  and 
electric  fields  in  the  two  regions  are  expanded  in  terms  of  the  eigenfunctions  obtained.  The 
coefficient  of  each  mode  is  determined  by  the  variational  method.'® 

+  -  u^)(Srff  -I- ST^  f  +  (Z3f  -  (8) 

+  +  SD^'”)*]dS  =  0 


where  the  surface  integral  is  performed  at  XjsO  over  one  unit  cell  ^  xj  <  — ),  W  and  Cm 

2  2 


refer  to  the  fluid  and  composite,  respectively,  and  *  denotes  the  complex  conjugate. 

For  the  situation  considered  here,  in  the  fluid  medium,  there  are  incident  wave  and  reflected 
wave,  therefore,  in  terms  of  the  eigenfunctions  in  the  fluid  medium,  the  elastic  displacement 
and  electric  potential  in  the  fluid  can  be  expressed  as 


«r  =7  I  h^smih^xOcxp(-jP^x^)R^ 

n=0 

^3  =Po  cxpUpQXz)-  s  P^ cosih^ Xi)exp(_-jp^ Xi)R^ 

«=0 


=  I  cos(/i„%)exp(^rx3)i?„^ 

n=0 


where  Pj"  =  yW  Po  txpUPoX‘i)corTesponds  to  the  incident  wave  normal  to  the 
interface.  are  the  coefficients  determined  by  equation  (8). 
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In  the  composite  region,  only  the  transmitted  wave  exists.  In  the  ceramic  phase,  the  elastic 
displacements  and  the  electric  potential  are 
ni  3 

M?=  Z  Z^/S/nfcos(/i5x,)exp(y)3„.r3M„ 

«=li=l 
fti  3 
/j=li=l 

m  3  (10) 

I  lR^itniCOS(h^iXi)e\pUPnX3)An 

n=l/=l 

and  in  the  polymer  phase,  they  are 

ffl  2 

^^3  ~  Z  Z^/ii.^/ ~~))cxp(y^„X3)A^ 
n=l/=l  2 

«r  =  ;  Z  IRnigniSinih^iixi  - ^))expUPnX3)An  (1 1) 

n=l/=l  2 

^  p  d 

<I>  =  ZC„  cosh(^„(xi--))exp(y^„X3)A„ 

«=1  2 

where  A,  are  coefficients  to  be  determined  from  eq.  (8).  In  the  expansions  (9),  (10)  and  (1 1), 
the  number  of  modes  used  are  determined  by  the  accuracy  needed.  Making  use  of  the 
constitutive  equations  in  the  ceramic  and  polymer  phases,  the  stress  and  electric  displacement  in 
the  composite  region  can  be  obtained.  Substituting  these  quantities  into  eq.  (8)  yields  the 
foUowing  algebraic  equation: 

(M,)(Rj)  =  (Y)  (12) 

where  (Rj)=(R(,,  R,,...Rq,  A,,...A„)^,  (Mjj)  is  a  (G+m+l)x(G+m+l)  matrix,  and  (Vp  is  a 

(G+m+l)xl  matrix.  In  the  problem  treated,  G=7  and  m=9  are  used.  From  eq.  (12),  the 
coefficients  Rj  and  A^  are  determined. 

(2.2)  Reflection  coefficient  and  input  acoustic  impedance  at  the  interface  and  field 
distributions  in  the  composite 
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In  this  section,  we  will  first  derive  the  reflection  coefficient  for  a  composite  in  contact  with 
water  based  on  the  results  derived  in  the  preceding  sections.  As  we  have  pointed  out,  in  the 

expression  (9),  Pq  txpijpQ  A:3)is  the  incident  wave  and  other  terms  are  from  the  reflected 
wave.  The  zeroth  order  reflected  wave,  i.e.,  n=0  and  ho^sO,  represents  the  normal  reflection. 
For  n^O,  in  the  frequency  of  interest,  i.e.,  below  the  first  lateral  mode  frequency,  it  can  be 

shown  that  is  nearly  equal  to  j  (2n7c/d)  which  implies  that  the  amplitude  of  these  waves  will 

decay  as  exp(-2nKXj/d)  as  the  waves  depart  fi'om  the  interface  and  hence  they  are  evanescent 

waves  which  are  important  only  in  the  boundary  region.  In  this  sense,  the  reflection  coefficient 
which  can  be  measured  experimentally  is  equal  to  -R,,. 

Figure  3(a)  shows  the  reflection  coefficient  derived  for  a  plane  acoustic  wave  on  the 
boundary  between  water  and  a  2-2  composite  in  which  the  volume  fraction  of  ceramic  phase  is 
44%.  At  lower  frequencies  where  the  periodicity  d  of  the  2-2  composite  in  the  x, -direction  is 
much  smaller  than  the  wavelength  of  the  incident  wave,  the  2-2  composite  behaves  like  an 
effective  homogeneous  medium,  and  the  reflection  coefficient  is  real  and  can  be  fully 
characterized  by  the  acoustic  impedance  of  the  two  materials,  i.e.,  the  effective  acoustic 
impedance  of  the  composite  derived  based  on  effective  medium  theory  and  the  characteristic 
impedance  of  water.  "  As  the  frequency  increases,  however,  the  vibration  displacement  at  the 
surface  of  the  2-2  composite  becomes  nonuniform,  which  implies  that  it  could  no  longer  be 
regarded  as  a  homogeneous  medium.  Hence,  the  reflection  coefficient  at  the  water  and 
composite  interface  will  be  dispersive  and  become  a  complex  number  as  shown  in  figure  3(a) 
where  the  magnitude  of  the  reflection  coefficient  decreases  gradually  and  meanwhile,  the  phase 

of  the  reflection  coefficient  becomes  larger  than  180°  and  increases  with  frequency. 

To  verify  the  results,  an  experiment  was  conducted  to  measure  the  reflection  coefficient 
from  the  interface  between  water  and  a  2-2  composite  of  44%  ceramic.  The  schematic  of  the 
experimental  set-up  is  shown  in  figure  4,  where  the  reflected  signals,  both  the  amplitude  and 
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time  delay,  from  a  2-2  composite  of  44%  ceramic  volume  and  a  standard  reflector  made  of  a 
stainless  steel  rod,  were  compared.  Since  the  acoustic  impedance  of  the  stainless  steel  standard 
is  know,  the  reflection  coefficient  from  the  water-composite  interface  can  be  determined.  To 
approximate  the  plane  wave  condition,  the  distance  between  the  transducer  and  the  water- 
composite  interface  is  smaller  than  the  radius  of  the  transducer.  Because  the  main  purpose  of 
this  experiment  is  to  determine  the  phase  angle  of  the  reflection  coefficient  from  the  water- 
composite  interface,  it  is  crucial  to  ensure  the  separation  between  the  probing  transducer  and 
the  composite  the  same  as  that  between  the  transducer  and  stainless  steel  standard. 

The  result  thus  obtained  is  shown  in  figure  3(a).  At  fiequencies  below  the  lateral  resonant 
mode  (fd  <  1  MHz  mm),  the  agreement  between  the  model  result  and  experimental  result  is 
quite  good,  especially  the  phase.  The  fluctuation  in  the  experimental  data  in  the  amplitude  of  R 
around  fd  —  0.65  (MHz  mm)  is  probably  due  to  the  existence  of  other  weak  modes  at  the 
composite  surface.  At  fipequencies  near  the  lateral  mode  (fd  —  1  MHz  nun),  even  though  the 
model  prediction  does  not  reproduce  the  experimental  data,  it  still  follows  the  trend  of  the  data. 
Near  a  resonance,  any  small  error  in  the  input  parameters  in  the  model,  especially  the  elastic 
loss  in  the  polymer  phase  which  is  not  included  in  the  model,  can  cause  large  effect.  This  could 
be  the  reason  for  the  discrepancy  between  the  model  results  and  experimental  data. 

In  figure  3(b),  the  reflection  coefficient  for  a  2-2  composite  with  15%  ceramic  content  is 
presented.  Qearly,  the  trend  is  quite  similar  to  that  in  figure  3(a),  although  the  magnitude  and 
the  phase  of  the  reflection  coefficient  are  more  sensitive  to  the  frequency  of  incident  wave 
which  is  apparently  due  to  the  small  volume  content  of  the  ceramic  phase  in  the  composite  and 
lower  lateral  mode  frequency  (fd  —  0.8  MHz*mm). 

One  important  parameter  related  closely  to  the  reflection  and  transmission  of  an  acoustic 
wave  at  an  interface  is  the  acoustic  impedance  of  the  media  at  the  two  sides  of  the  boundary 
which  also  depends  on  the  nature  of  the  wave  such  as  plane  wave  or  spherical  wave.^  ’^  For 
homogeneous  isotropic  materials,  the  acoustic  impedance  for  a  plane  wave  is  the  characteristic 

acoustic  impedance  which  is  simply  the  product  of  the  mass  density  p  and  the  wave 
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(longitudinal  wave  or  shear  wave)  velocity  V.'-^  ’’  For  nonuniform  materials  (such  as 
piezocomposite  materials),  on  the  other  hand,  no  characteristic  impedance  can  be  defined  in 
such  a  simple  manner  due  to  the  dispersive  nature  of  the  properties  as  has  been  shown. 
However,  the  acoustic  impedance  of  a  material  can  still  be  found,  for  example,  from  the 
reflection  coefficient  or  other  methods.  It  should  be  pointed  out  that  although  there  exist 
several  definitions  for  determining  the  acoustic  impedance  at  an  interface  (input  acoustic 
impedance),  for  a  heterogeneous  material,  the  results  obtained  by  using  different  definitions 
may  not  be  the  same.'^  '*  However,  since  the  input  acoustic  impedance  is  not  a  direct  physical 
quantity  but  rather  a  parameter  introduced  for  the  convenience  of  the  analysis,  a  situation  in 
analogy  to  the  complex  notation  introduced  in  many  engineering  fields,  which  definition  is 
more  appropriate  really  depends  on  the  situation  where  it  is  used.  In  dealing  with  the  energy 
transfer  across  an  interface,  one  might  have  to  use  the  formula  derived  under  energy  flow 
consideration.  On  the  other  hand,  in  designing  matching  layers  and  dealing  with  the  reflection 
and  transmission  of  acoustic  waves  in  a  multilayer  medium,  the  amplitude  and  phase  of  the 
wave  components  are  crucial  which  may  not  be  included  in  the  coefficients  related  to  the 
transmission  and  reflection  of  the  acoustic  energy.  In  the  discussion  here,  we  will  adopt  the 
definition  related  to  the  wave  reflection  and  transmission  problems. 

From  the  reflection  coefficient  R,  the  input  acoustic  impedance  of  the  composite  at  the 
interface  can  be  found  as 

Z,=1^ZW  (.3) 

where  Z„  is  the  effective  input  impedance  of  the  2-2  composite,  is  the  characteristic 
impedance  of  water.  Using  the  results  in  figure  3,  the  effective  input  acoustic  impedance  of  2-2 
composites  with  44%  and  15%  ceramic  volume  content  is  evaluated  and  presented  in  figures 
5(a)  and  5(b).  Apparently,  there  is  a  large  change  of  the  magnimde  and  phase  of  the  input 
acoustic  impedance  as  the  frequency  increases. 
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In  analogy  to  the  electrical  impedance,  the  specific  impedance  at  the  interface  can  also  be 
found  from  the  ratio  of  the  stress  (in  analogy  to  the  voltage)  to  the  displacement  velocity  (in 
analogy  to  the  current): 


Z(x„0)  = 


T/^(Xi.O) 

vf’(x„  0) 


(14) 


where  v^^  is  the  particle  displacement  velocity  in  the  2-2  composite.  Eq.  (14)  is  used  widely 

in  the  equivalent  circuit  model  of  transducers.'*’'*  Apparently,  for  a  2-2  composite  considered, 
Z  from  eq.  (14)  is  a  function  of  x,  due  to  the  variation  of  T3  and  Vj  in  that  direction.  To 
eliminate  this  variation,  the  approach  taken  by  Miller  and  Pursey  is  adopted  here.'*  That  is,  the 
averaged  T3  and  averaged  V3  in  the  x,-direction  are  used: 

,  _  JT°"(x„0)<ic, 

-  Jvr(x„0)dx, 

where  the  integration  is  taken  in  one  unit  cell.  This  expression  seems  physically  meaningful 
since  it  reflects  an  averaged  mechanical  impedance  of  the  composite  at  the  interface. 
calculated  from  eq.  (15)  for  composites  with  44%  and  15%  ceramic  contents  is  shown  in  figure 
6  which  is  quite  close  to  those  in  figure  5.  In  the  following  discussion,  eq.  (15)  is  used  to 
calculate  Z„  in  the  composite. 

One  interesting  question  is  whether  the  acoustic  impedance  of  the  fluid  phase  will  affect  the 
effective  input  impedance  of  the  composite  at  the  interface.  To  answer  this  question,  the 
effective  input  impedance  of  the  composite  at  the  interface  is  evaluated  for  fluid  phases  with 
different  acoustic  impedance  (by  varying  the  bulk  modulus)  and  the  results  are  presented  in 
figure  7.  Clearly,  the  effective  acoustic  impedance  of  a  2-2  composite  is  independent  of  the 
acoustic  impedance  of  the  fluid  at  the  interface. 

In  order  to  shed  light  on  the  large  change  of  the  input  acoustic  impedance  with  frequency, 
the  surface  displacement  U3  is  evaluated  at  the  center  of  the  ceramic  plate  (x,=0)  and  polymer 
plate  (x,=d/2)  and  the  results  are  presented  in  figure  8(a).  At  low  frequencies,  the  surface 
displacements  in  the  polymer  and  ceramic  regions  are  in  unison  which  indicate  that  the  isostrain 
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model  used  in  many  earlier  modelings  on  ultrasonic  composite  transducers  is  valid  in  this 
frequency  region."  '^  As  the  frequency  increases,  even  at  frequencies  still  far  below  the  first 
lateral  mode  (at  fd  near  1),  the  surface  vibration  amplitudes  in  the  two  regions  become  quite 
different.  At  the  frequency  near  the  first  lateral  mode,  the  vibrations  in  the  polymer  and 
ceramic  regions  are  180  out  of  phase  and  in  this  frequency  region,  the  effective  input  acoustic 
impedance  of  the  composite  is  very  near  that  of  water  and  the  reflection  coefficient  exhibits  a 
minimum  as  shown  in  figure  3. 

As  has  been  pointed  out  in  the  introduction,  for  an  ultrasonic  piezoceramic  polymer 
composite,  one  of  the  most  important  factors  in  determining  the  performance  is  the  effective 
acoustic  energy  exchange  between  the  ceramic  and  the  external  medium.  As  seen  in  figure  3 
and  5  where  at  frequencies  near  the  lateral  mode,  the  input  acoustic  impedance  of  the  composite 
is  close  to  that  of  water  and  the  reflection  coefficient  reaches  a  minimum,  the  question  is  how 
much  of  that  energy  enters  into  the  ceramic  plates.  Figure  8(b)  shows  how  the  total  acoustic 
powers  entering  into  the  ceramic  region  and  polymer  region  vary  with  frequency.  Clearly,  at 
the  high  frequency  region,  in  spite  of  the  fact  that  the  effective  transmission  coefficient  of  the 
composite  increases,  the  amount  entering  into  the  ceramic  plate  actually  decreases. 

It  is  also  interesting  to  examine  how  this  acoustic  energy  distribution  in  the  polymer  and 
ceramic  regions  changes  along  the  Xj-direction  since  the  interaction  between  the  ceramic  and 
polymer  through  the  joint  region  will  cause  the  acoustic  energy  transfer  between  the  two.  In 
figure  9,  the  acoustic  energy  distribution  along  the  Xj-direction  at  three  frequencies  ( f  «  f,,  f  - 
-  f,  12,  and  f  —  f,,  where  f,  is  the  first  lateral  mode  frequency)  is  shown.  In  all  the  cases,  there 
is  a  redistribution  of  the  acoustic  energy  along  the  Xj-axis  and  the  acoustic  energy  in  the 
polymer  is  gradually  transferred  to  the  ceramic  plate.  However,  at  low  frequencies,  the 
transition  region  is  much  shorter  than  that  at  high  frequencies  (with  respect  to  the  wave  length 

at  that  frequency).  For  instance,  in  figure  9(a),  the  reduced  length  Xj/X  is  about  0.01  for  the 

acoustic  power  in  the  ceramic  to  reach  90%  of  the  final  value,  while  in  figure  9(b),  this  region 
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increases  to  about  0. 1 ,  and  in  figure  9(c),  it  becomes  near  0.2.  It  should  be  reminded  that  in  a 
thickness  mode  transducer,  the  thickness  of  the  piezocomposite  is  X/2. 

In  figure  10,  the  change  of  the  displacement  Uj  in  the  polymer  and  ceramic  regions  along  the 
Xj-direction  at  the  three  frequencies  is  presented  and  it  shows  that  at  high  frequencies,  the 
displacement  amplitude  in  the  polymer  and  ceramic  regions  is  no  longer  the  same  even  deep 
inside  the  composite. 

The  change  of  the  input  acoustic  impedance  from  eq.  (15)  along  the  Xj-direction  in  the 
composite  is  also  evaluated.  Figure  1 1  are  the  results  of  the  2-2  composite  with  44%  ceramic 
content  at  different  frequencies.  Apparently,  at  high  frequencies  there  is  a  large  change  of  the 
effective  acoustic  impedance  from  the  surface  into  the  interior  of  the  composites  and  even  in  the 
interior  of  the  composite,  the  effective  impedance  is  not  the  same  as  that  derived  from  the 
effective  medium  theory."  '* 

in.  W ave  reflection  and  transmission  at  solid-2-2  composite  interfaces 

In  general,  a  piezoelectric  ultrasonic  traiisducer  has  a  multilayer  structure  in  which  there  are 
one  or  two  quarter  wavelength  impedance  matching  layers  in  front  of  the  piezoelectric  element 
and  a  backing  material  in  the  back.  The  piezoelectric  element  in  the  transducer  does  not  contact 
direcdy  with  fluid  medium.  Therefore,  the  investigation  of  the  interaction  between  a  2-2 
composite  and  solid  medium  has  practical  importance  for  composite  transducers. 

The  issues  investigated  in  this  section  are  quite  similar  to  those  studied  in  the  preceding 
section.  The  difference  is  that  a  solid  medium  can  support  shear  waves  while  a  fluid  cannot. 
Due  to  this  difference,  the  boundary  conditions  (3)  must  be  modified  to: 

Tff  =  Tji  Tg"  = 

=  4.^  Cf”  =  Di.  ‘  ' 

where  the  superscript  S  refers  to  the  variables  in  the  solid  medium.  Hence,  the  variational 
formula  for  the  problem  becomes 


(17) 


-  J35”’)(&3  +&?')• +(r,!  -  r,f'")(5,f'”  +&f)* 

+(U3^" -uixsifi" +sri,f +(uf^"‘ -ufmS” +sr^\)‘ 

+(£^  -Zf”)(M>^  +«>'^”)'  +(0^  -$'^'’■)(^Df  +SD^'")*]ilS.  =  0 
where  the  surface  integral  is  performed  at  X3=0  over  one  repeat  unit  in  the  x, -direction  of  the 

composite.  The  elastic  displacement  wf,  1/3  ,  and  the  electrical  potential  <I>^in  the  solid  medium 
are  expanded  in  terms  of  the  eigenfuntions: 


2 

wf  =J1  I^mSin()i^xj)exp(-;jS^X3)/?^ 

n  /=I 

2 

«3  =  Po  exp(7^^X3)-7l  X/„^  cos()i^xi)exp(-yi8^X3)/?^ 

n  /=1 


(18) 


=  Icos()i;Jxi)exp(;i^x3)i?^ 

n 


where  is  the  incident  wave  and  other  terms  are  from  the  reflected  waves  in 

which  both  the  longitudinal  and  shear  components  exist  except  the  one  with  n=0  where  the 
reflected  wave  is  a  pure  longitudind  one.  Substituting  the  elastic  and  electrical  variables  in  the  2- 
2  composite  and  eq.  (18)  into  equation  (17),  the  coefficients  of  the  expansions  can  be  obtained. 

Figure  12  presents  the  dispersion  curve  of  the  reflection  coefficient  at  the  interface  for  2-2 
piezocomposites  with  44%  ceramic  volume  content  and  the  acoustic  impedance  calculated  from 
eq.  (13)  where  the  characteristic  impedance  of  the  solid  medium  is  Z=2.4  MRayls. 

By  comparing  the  results  in  figure  12  and  figure  5,  it  can  be  found  that  the  input  acoustic 
impedance  of  2-2  composites  calculated  here  is  different  from  that  in  the  fluid  case.  In 
addition,  it  is  also  found  that  at  a  solid  medium-2-2  composite  interface,  the  input  acoustic 
impedance  also  depends  on  the  characteristic  impedance  of  the  solid  medium  at  the  interface, 
which  seems  to  be  quite  different  from  the  fluid  medium  interface. 

To  understand  this  phenomenon,  we  notice  that  the  difference  between  a  solid  medium  and 
a  fluid  one  is  that  a  solid  medium  has  a  non-zero  shear  stiffness  constant.  Therefore,  Z,„  for  a 
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2-2  composite  (44%  ceramic  content)  in  contact  with  solid  medium  is  evaluated  where  p  (mass 

density)  and  c,,  are  kept  as  constant  and  C44  is  varied  and  the  results  are  presented  in  figure 
13(a).  Similarly,  Zj„  for  a  2-2  composite  in  contact  with  solid  medium  is  also  evaluated  where 

p  and  C44  are  kept  constant  and  c,,  is  varied  and  the  results  are  shown  in  figure  13(b).  The 

results  demonstrate  that  for  a  plane  incident  wave,  even  if  the  characteristic  longitudinal 
impedance  of  the  solid  medium  is  kept  constant,  the  effective  input  impedance  of  a  composite 
and,  hence,  the  reflection  coefficient  from  the  interface  will  change  if  the  shear  stiffness 
coefficient  of  the  solid  medium  changes.  On  the  other  hand,  if  the  shear  stiffness  constant  is 
kept  constant  and  the  characteristic  longitudinal  impedance  is  varied  in  the  solid  medium,  the 
effective  input  acoustic  impedance  of  a  composite  will  not  have  much  change.  This 
phenomenon  is  directly  related  to  the  surface  uniformity  of  the  composite  and  as  has  been 
shown  earlier,  the  enhanced  stress  transfer  between  the  ceramic  and  polymer  regions  is  through 
the  shear  action.'® 

To  illustrate  this,  the  surface  vibration  distribution  in  the  ceramic  and  polymer  regions  is 
also  evaluated  for  the  2-2  composite  in  contact  with  solid  media  of  different  c,,  and  C44  and  the 
results  are  presented  in  figure  14.  There  is  very  little  difference  in  the  surface  vibration  profile 
for  solid  media  with  the  same  C44  and  different  c,,.  In  contrast,  the  surface  vibration  profile 
will  change  as  C44  is  changed  (figures  14(a)  and  14(b)). 

Figure  15  illustrates  the  input  acoustic  impedance  of  a  composite  as  a  function  of  X3  for 
different  solid  media  at  the  interface.  It  shows  the  length  of  the  transition  region  in  the 
composite  where  the  evanescent  waves  are  important  decreases  as  the  acoustic  impedance  of 
the  medium,  especially  the  shear  modulus,  increases. 

rv.  Reflection  and  transmission  from  the  boundary  between  fluid  and  acoustic  impedance 
matching  layer 
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We  now  proceed  to  investigate  how  the  various  effects  observed  in  the  preceding  sections 
affect  the  selection  of  anti-reflection  matching  layers  at  the  interface.  Figure  16  is  a  schematic 
of  an  acoustic  system  which  consists  of  a  fluid  medium,  a  matching  layer  and  a  2-2  composite. 
In  this  case,  there  are  two  interfaces  with  different  boundary  conditions  to  be  treated.  One  of 
the  interfaces  is  between  the  2-2  composite  and  matching  layer  where  the  boundary  conditions 
are  those  of  equation  (16).  Another  one  is  the  fluid  and  matching  layer  interface  where  the 
boundary  conditions  are: 


T-S  rrW  ’pS  t\  s  w  r,S  T-,W 

^3  =  ^33 . 7i3  =  0,  M3  =  M3  ,0  =  d>  ,  at  XjS-t 


(19) 


where  t  is  the  thickness  of  the  matching  layer  and  the  other  notations  are  the  same  as  before. 
The  variational  formula  for  this  problem  can  be  derived: 


Js,  ^(^3  -  J33”  X&l  +  &3^" )’  +  (3ii  -  X&f" 


■Kuf"  -ufxffjf” +ar3^)’  Huf"  -afxfflif"' +aif  )• 

+  (D^-  +  53)^'”)*  +  (0)^  -  <I>^'”)(^D^  -t-  a3f'”)*]</5 

+Js.  ^(533  -  Jf3Xa-f + SuY)' + 2r,i(&f  )• + (ul  -  u^KSif, + srS)' 

+  (D^  -  +  (3>'^  -  4>^XSDi  +  SC^)*]dS  =  0 

where  the  surface  integrals  S,  and  Sj  are  performed  at  x,=0  and  Xj^-t  planes,  respectively, 
from  x,=0  to  x,=d/2. 

In  the  matching  layer  region,  there  are  waves  which  propagate  in  both  the  -hXj  and  -X3 
directions  due  to  the  two  interfaces  at  XjSsO  and  XjSt.  Hence,  the  expressions  for  the  elastic  and 
electrical  variables  become; 


2 

uf  =  JI  Ig«/(i3fi)sin(/iJxi)[-exp(-7)S^X3)/?i  -hexp(j)3^jX3)/?2] 
n  f=l 
2 

«3  =  I  Ifni(^ni)cos(h^xi)[-exp(-Jfi^iX3)/ii  +  expUPmX2)J^]  (21) 
n  £=1 

=  Icos(/j^xi)[-exp(-/^j:3)/^  -Hexp(;^X3)/?^] 
n 
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The  solutions  in  the  composite  and  in  the  fluid  medium  are  those  of  eqs.  (9),  (10)  and  (11). 
Substituting  those  expressions  into  equation  (20),  the  undetermined  coefficients  in  these 
equations  can  be  found  and  various  quantities  in  different  regions  can  be  evaluated. 

In  general,  for  homogeneous  materials,  the  reflection  coefficient  from  the  boundary  between 
a  fluid  and  a  matching  layer  will  be  zero  when  the  thickness  of  the  matching  layer  is  quarter 
wavelength  and  its  acoustic  impedance  Z“  is 

Z^=Vz^  (22) 

where  7^  and  Z®  are  the  acoustic  impedance  of  the  media  on  the  two  sides  of  the  matching 
layer,  respectively.*'^ 

For  a  composite  material,  from  the  results  obtained  in  the  preceding  sections,  it  is  expected 
that  the  parameters  of  the  anti-reflection  matching  layer  will  be  modified.  Figure  17  illustrates 
the  derived  reflection  coefficient  from  the  boundary  between  the  fluid  and  matching  layer  for 
different  matching  layer  materials  where  the  thickness  of  the  matching  layer  is  chose  so  that  the 
ftequency  of  the  reflection  coefficient  minimum  is  at  fd=0.335  MHz*mm.  From  results,  it  is 
found  that  the  acoustic  impedance  of  the  matching  layer  from  which  the  reflection  coefficient 
becomes  zero  satisfies  approximately  equation  (22)  if  Z®  is  replaced  by  the  effective  input 
acoustic  impedance  of  the  2-2  composite  at  the  low  frequency.  However,  the  thickness  of  the 

matching  layer  is  no  longer  quarter  wavelength  and  in  fact,  it  is  shorter  than  X/4.  These  results 

are  summarized  in  figure  18.  The  data  in  figure  18(a)  is  the  acoustic  impedance  of  the 
matching  layer  in  order  for  the  reflection  coefficient  to  be  zero  as  a  function  of  frequency.  In 
figure  18(b),  the  change  of  the  matching  layer  thickness  as  a  function  of  frequency  is  presented 

where  t=l  corresponds  to  the  thickness  of  X/4  at  that  frequency.  Therefore,  as  the  frequency 

increases,  both  the  acoustic  impedance  and  the  thickness  of  the  matching  layer  decrease.  For 
example,  when  the  frequency  changes  from  0.05  MHz*mm  to  0.5  MHz*mm,  the  acoustic 
impedance  of  the  matching  layer  changes  from  4.70  MRayls  to  4.45  MRayls  and  the  thickness 

t  changes  from  about  X/4  to  0.86*X/4.  The  reduction  in  the  matching  layer  thickness  is  a 
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direct  result  of  the  fact  that  the  incident  plane  wave  suffers  more  than  180°  phase  loss  upon  the 
reflection  from  the  interface. 

From  the  fact  that  in  the  frequency  range  of  interesting,  the  evanescent  waves  at  the 
interface  will  decay  approximately  as  exp(-2nTCX3/d),  the  effect  of  the  interface  of  Xjst, 


where  t  is  approximately  X/4,  to  the  surface  vibration  distribution  at  the  composite¬ 


matching  layer  interface  will  be  quite  small.  Therefore,  the  reflection  coefficient  introduced 
in  the  section  in  between  the  composite  and  a  solid  medium  can  be  used  to  approximate  the 
reflection  at  the  matching  layer-composite  interface.  Hence,  the  reflection  coefficient  R 
from  the  system  in  figure  16  can  be  approximated  as: 


R  +  R 

p  _  *^12^  ^23^ 


2jP|r 


(23) 


1+  R12R23® 

where  R,2  and  R23  are  the  reflection  coefficients  at  the  fluid-matching  layer  interface  and 
matching  layer-2-2  composite  interface,  respectively.  In  eq.  (23),  R,2  is  the  reflection 
coefficient  from  the  interface  of  two  semi-infinite  media  of  the  fluid  and  solid,  and  R23  is  that 
derived  in  the  section  HI  (the  reflection  coefficient  from  the  solid  medium-composite  interface). 
From  equation  (23),  the  requirements  to  the  matching  layer  to  achieve  the  total  transmission 
(R=0)  can  be  derived.  The  acoustic  impedance  of  the  matching  layer  is 


rM 


f'yW  1^7  Cm| 

ZCm 

-z'^cos(e) 

r  1 

ZCm 

cos(0)-z'^ 

(24) 


where  0  is  the  phase  angle  of  the  acoustic  impedance  of  the  2-2  composite.  The  thickness  of 


the  matching  layer  is 


t 


1  2Z“  Z°"|sm(«) 

^arctan(  .  '  7 

n  Z^l  -(Z'^)^ 


(25) 


In  the  frequency  range  of  interest,  the  results  from  eqs.  (24)  and  (25)  are  almost  identical  to 
those  derived  from  the  numerical  calculation  presented  in  figure  18. 
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In  order  to  verify  the  theoretical  predictions,  the  reflection  coefficient  from  a  fluid¬ 
matching  layer  interface  was  measured  as  a  function  of  frequency  utilizing  the  experimental 
set-up  shown  in  figure  4  where  a  stainless  steel  standard  was  used  as  the  reference.  For  a 
2-2  composite  with  44%  PZT-5H  ceramic  in  a  Spurr  epoxy  matrix,  the  required  acoustic 
impedance  for  the  matching  layer  in  frequencies  near  f  d  =  0.5  (MHz  mm)  is  about  4.45 
MRayls.  A  series  of  polymer  were  examined  and  the  acoustic  velocity  was  determined 
from  the  time  of  flight  technique.  Combined  with  the  data  on  the  density,  the  acoustic 
impedance  of  the  material  can  be  obtained.  A  Hysol  epoxy  (EE1068)  was  chosen  because 

the  acoustic  impedance  Z  =  4.68  Mrayls  (p  =  1610  kg/m2,  c„  =  1.3634x10'®  N/m^  and 

c^  =  3.432x10®  N/m-)  is  quite  close  to  the  required  value.^®  For  a  fixed  thickness  t  of  the 
matching  layer,  the  frequency  position  of  the  minimum  reflection  can  be  determined  and  the 

ratio  of  AxTk  is  obtained.  By  varying  the  thickness  of  the  matching  layer,  the  dependence  of 

4tA,on  the  frequency  can  be  obtained.  The  data  is  presented  in  Fig.  18(b)  (open  circles) 

and  apparently,  there  is  an  excellent  accord  between  the  theoretical  prediction  and 
experimental  data. 

V.  Summary  and  Acknowledgment: 

The  reflection  and  transmission  of  a  plane  acoustic  wave  at  a  medium-composite  interface 
and  the  issues  related  to  the  design  of  matching  layer  for  a  composite  are  analyzed  based  on  an 
approach  developed  recently,  which  can  address  explicitly  the  non-uniform  field  distributions 
due  to  the  heterogeneity  structure  of  a  piezoceramic  polymer  composite.  It  has  been  shown  that 
the  reflection  coefficient  from  the  interface  is  a  complex  and  the  reflected  wave  suffers  a  more 
than  180“  phase  change.  The  effective  input  acoustic  impedance  Zj^  of  the  composite  at  the 
interface  was  evaluated  and  both  the  amplitude  and  phase  show  a  strong  frequency 
dependence.  For  a  fluid  medium,  it  was  found  that  Z;,  does  not  change  with  the  acoustic 
impedance  of  the  medium.  However,  for  a  solid  medium,  Z,„  will  change  if  the  shear  stiffness 
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constant  of  the  medium  changes.  It  was  demonstrated  that  this  difference  originates  from  the 
non-uniformity  of  the  surface  vibration  distribution  of  the  composite  at  the  interface  which 
depends  crucially  on  the  shear  stiffness  constant  of  the  medium. 

Since  for  a  piezoceramic  polymer  composite,  it  is  the  ceramic  phase  which  performs  the 
energy  conversion  between  the  acoustic  and  electric  forms,  how  much  acoustic  energy  can 
enter  the  ceramic  region  is  one  of  the  most  important  parameters  in  a  composite  transducer 
design.  In  the  paper,  we  show  that  even  though  the  effective  transmission  coefficient  increases 
as  the  frequency  is  increased,  the  amount  of  acoustic  energy  entering  the  ceramic  region 
,  actually  decreases.  Therefore,  there  may  be  a  trade-off  between  the  bandwidth,  which  is 
related  to  the  transmission  coefficient,  and  the  sensitivity  in  the  composite  transducer  design. 

From  the  fact  that  there  is  more  than  180°  phase  change  in  the  reflection  from  the  medium- 
composite  interface,  it  is  shown  that  the  matching  layer  thickness  is  no  longer  equal  to  the 
quarter  wavelength  but  smaller  than  that.  In  addition,  the  acoustic  impedance  of  the  matching 
layer  will  also  be  affected  by  the  phase  of  7.^  of  the  composite,  although  the  effect  is  not  as 
sigmficant  as  that  in  the  thickness  of  the  matching  layer. 

Reflection  experiment  was  conducted  from  both  the  fluid-composite  interface  and  Ae  fluid¬ 
matching  layer  interface  and  the  agreement  between  the  experimental  results  and  theoretical 
prediction  is  quite  good. 

The  authors  wish  to  thank  Dr.  J.  Yuan  of  ATL/Echo  Ultrasound  for  the  assistance  during 
this  investigation.  This  work  was  supported  by  the  Office  of  Naval  Research  under  the  Grant 
No:  N00014-96-0357. 
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Figure  captions: 

Figure  1,  (a)  Schematic  of  a  2-2  piezoceramic  polymer  composite  where  the  dimensions  in 
both  the  X,  and  Xj  directions  are  much  larger  than  the  period  d.  (b)  Schematic  of  an  interface 
between  a  fluid  and  a  2-2  piezocomposite  where  a  plane  acoustic  wave  normally  incidents  to 
the  interface  from  the  fluid  medium.  The  composite  occupies  the  semi-infinite  space  Xj  >  0. 

Figure  2.  The  dispersion  curves  for  the  2-2  composite  of  PZT-5H  ceramic  and  Spurrs 
epoxy  where  the  ceramic  content  is  44%. 

Figure  3.  The  reflection  coefficient  (both  the  amplitude  and  phase)  as  a  function  of 
frequency  for  2-2  composites  with  (a)  44%  ceramic  and  (b)  15%  ceramic  content.  The 
anomalous  changes  at  fd  =  1  in  (a)  and  fd  =  0.8  in  (b)  are  related  to  the  lateral  resonant  mode  in 
these  composites,  respectively.  The  experimental  data  are  shown  in  (a)  where  the  open  circles 
are  the  anqilitude  and  solid  circles  are  the  phase  of  the  measured  reflection  coefficient. 

Figure  4.  Schematic  of  the  set-up  used  to  measure  the  reflection  coefficient  of  2-2 
composites.  In  order  to  obtain  phase  information  of  the  reflection  coefficient,  the  distance 
between  the  probe  transducer  and  the  2-2  composite  is  maintained  the  same  as  that  between  the 
probe  and  the  stainless  steel  standard. 

Figure  5.  (a)  The  effective  input  acoustic  impedance  Z;,  of  the  2-2  composite  with  44% 
ceramic  content  calculated  from  the  data  in  figure  3  (the  reflection  coefficient),  (b)  for  the  2- 
2  composite  with  15%  ceramic  volume  content. 

Figure  6.  The  effective  input  impedance  Z,„  calculated  from  eq.  (16)  for  (a)  the  composite 
of  44%  ceramic  volume  content  and  (b)  the  composite  of  15%  ceramic  content.  The  results 
here  are  nearly  the  same  as  those  in  figure  5  except  at  frequencies  above  the  first  lateral  mode. 

Figure  7.  Zi„  of  the  composite  of  44%  volume  content  derived  from  the  reflection 
coefficient  for  fluid  media  with  different  acoustic  impedance:  Z,  =  1.48  Mrayls  (water),  and 
Z2=  2.09  Mrayls.  Apparently,  Zj„  of  the  composite  for  the  two  cases  overlaps  with  each  other 
indicating  that  it  does  not  change  with  the  acoustic  impedance  of  the  fluid  medium. 
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Figure  8.  (a)  The  surface  vibration  distribution  as  a  function  of  frequency  for  the  2-2 

composite  with  44%  ceramic  content  at  the  interface  with  water  and  under  an  acoustic  field.  Uj" 
and  u/  are  the  displacement  at  the  center  of  the  polymer  plate  (x,=d/2)  and  ceramic  plate 
(x,ssO).  (b)  The  acoustic  energy  (in  reduced  units)  entering  the  ceramic  and  polymer  regions 

as  a  function  of  frequency  evaluated  at  the  interface  (x3=0)  where  is  the  total  energy  in  the 
incident  wave. 

Figure  9.  The  energy  redistribution  inside  the  composite  along  the  Xj-direction  where 
X3=0  is  the  interface.  There  is  an  energy  transfer  from  the  polymer  to  the  ceramic  inside  the 
^  composite.  As  shown  in  the  figures,  the  width  of  the  transition  region  increases  as  the 
frequency  increases  (in  reduced  length  unit). 

Figure  10.  The  vibration  profile  as  a  function  of  the  distance  from  the  interface  with  water 
(XjsO)  where  Uj**  and  Uj'  are  the  displacement  at  the  center  of  the  polymer  plate  (x,=d/2)  and 
ceramic  plate  (x,=0).  At  low  frequencies,  the  polymer  and  ceramic  vibrate  with  the  same 
amplitude  and  phase  while  at  high  frequencies,  even  far  away  from  the  interface,  the  vibration 
amplitude  of  the  two  is  still  not  the  same  while  the  phase  becomes  the  same. 

Figure  1 1.  The  effective  input  impedance  from  eq.  (16)  as  a  function  of  the  distance  from 
the  interface  with  water  (X3=0)  at  three  typical  frequencies:  far  below  the  lateral  mode,  at  about 
the  half  of  the  lateral  mode  frequency,  and  near  the  lateral  mode  frequency. 

Figure  12.  (a)  The  reflection  coefficient  at  a  interface  between  a  solid  medium  (Z=2.4 
MRayls)  and  a  2-2  composite  of  44%  ceramic  volume  content  as  a  function  of  frequency,  (b) 
The  effective  input  acoustic  imepdance  of  a  2-2  composite  of  44%  ceramic  volume  content  as  a 
function  of  frequency.  Z,„  is  calculated  from  the  reflection  coefficient. 

Figure  13.  Data  illustrate  the  dependence  of  the  input  acoustic  impedance  of  a  composite 
on  the  elastic  properties  of  the  solid  medium  at  the  interface,  (a)  The  solid  medium  for  the 

curve  1  is  p=1.61  g/cm\  c„  =  1.36*10'°  N/m%  and  c«  =  3.43*10’  N/m',  and  for  the  curve  2 
and  curve  3,  both  p  and  c„  are  kept  constant  and  c.^,  =  6.86*10’  N/m*  and  c^^  =  1.72*10’ 
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N/m^  (b)  Here  both  p  and  C44  are  kept  constant,  and  c,,=  0.68  *10'®,  c,,=  1.36  *10'“,  and 

c,,=  2.72  *10'“  N/m^  The  amplitude  of  does  not  change  much  with  c,,  except  a  small 
change  in  the  phase. 

Figure  14.  The  figures  illustrate  the  effect  of  c,,  and  c^  of  the  solid  medium  on  the  surface 
vibration  profile  of  the  composite,  (a)  The  parameters  of  the  solid  medium  for  the  data  of 

dashed  curve  is  p=1.61  g/cm^  c„  =  1.36*10'“  N/m%  and  C44  =  3.43*10“  N/m^  and  for  the 
solid  curve,  both  p,  C44  are  the  same  except  c,, =2.72*  10'“  N/ml  (b)  The  parameters  for  the 

solid  medium  is  p=  1.61  g/cm^  c,,  =  1.36*10'“  N/m^  and  c^  =  6.86*10“  N/m^  Notice  the 

effect  of  C44  on  the  surface  vibration  profile  of  the  composite. 

Figure  15.  The  figures  illustrate  the  effect  of  the  elastic  properties  of  the  solid  medium  at 
the  interface  on  the  effective  input  impedance  distribution  along  the  Xj-axis  at  a  fi^quency  f*d  = 

0.446  Mhz*mm  where  X3=0  is  the  interface,  (a)  The  dashed  line:  p=1.16  g/cm^  c„= 

7.72*10“  N/m^  and  c^^  =1.588*10“  N/ml  The  soUd  line:  p=l.l  g/cm^  c„  =  5.41*10“  N/m^ 
and  C44  =  1.307*10“  N/m^  (Spurrs  epoxy),  (b)  The  parameters  for  the  solid  medium  is: 
p=1.61  g/cm^  c,,  =  1.36*10'“  N/m^  and  =  3.43*10“  N/ml 

Figure  16.  Schematic  of  an  interface  system  consisting  of  a  fluid  medium  (medium  1),  a 
matching  layer,  and  a  2-2  composite.  The  thickness  of  the  matching  layer  is  t  and  the  period  of 
the  composite  is  d. 

Figure  17.  The  reflection  coefficient  from  the  matching  layer  system  where  the  total 
transmission  occurs  at  f*d  =  0.334  MHz*mm  for  the  matching  layer  with  t  =  2.02*d  = 

0.91*(X/4)  (at  f*d  =  0.334  MHz*mm)  and  Z'^=4.68  Mrayls  (p=1.61  g/cm^  c,,  =  1.36*10'“ 

N/m^,  and  C44  =  3.43*10“  N/m^)  (the  curve  1).  For  the  comparison,  the  reflection  coefficient 
from  the  matching  layers  with  other  parameters  is  also  shown  here:  for  the  curve  2,  Z“  =  3.0 
Mrayls  and  t  =  1.85*d,  and  for  the  curve  3,  =  2.44  Mrayls  and  t  =  1.58*d. 
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Figure  1 8.  The  required  parameters  for  (a)  and  (b)  t  of  the  matching  layer  to  have  total 
transmission  as  a  function  of  frequency  for  a  water-2-2  composite  (44%  ceramic  volume 
content)  interface  (solid  lines  are  the  theoretical  predictions).  The  experimental  data  was 
presented  in  (b)  as  open  circles. 
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ACTUATOR  STUDIES 


APPENDIX  51 


ULTRAHIGH  STRAIN  AND  PIEZOELECTRIC  BEHAVIOR  IN 
RELAXOR  BASED  FERROELECTRIC  SINGLE  CRYSTALS 


Seung-Eek  Park  and  Thomas  R.  Shrout 

Materials  Research  Laboratory,  The  Pennsylvania  State  University,  University  Park,  PA  16802 

Abstract 

The  piezoelectric  properties  of  relaxor  based  ferroelectric  single  crystals,  such  as 
Pb(Zn,/3Nb2/3)03  -  PbTi03  (PZN-PT)  and  Pb(Mgi/3Nbw)03  -  PbTi03  (PMN-PT)  were 
investigated  for  electromechanical  actuators.  In  contrast  to  polyciystalline  materials  such  as 
Pb(Zr,Ti)03  (PZT’s),  morphotropic  phase  boundary  (MPB)  compositions  were  not  essential  for 
high  piezoelectric  strain.  Piezoelectric  coefficients  (d33’s)  >  2500  pC/N  and  subsequent  strain 
levels  up  to  >0.6%  with  minimal  hysteresis  were  observed.  Crystallographically,  high  strains  are 
achieved  for  <001>  oriented  rhombohedral  crystals,  though  <11 1>  is  the  polar  direction. 
Ultrahigh  strain  levels  up  to  1.7%,  an  order  of  magnitude  larger  than  those  available  from 
conventional  piezoelectric  and  electrostrictive  ceramics  could  be  achieved,  being  related  to  an  E- 
field  induced  phase  transformation.  High  electromechanical  coupling  (k33)  >  90%  and  low 
dielectric  loss  <1%,  along  with  large  strain  make  these  crystals  promising  candidates  for  high 
performance  solid  state  actuators. 

Keywords  :  Relaxor  ferroelectrics,  single  crystal,  actuator,  piezoelectrics,  phase 
transformations,  hysteresis,  domains. 
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1. Introduction 


Electromechanical  actuators  directly  transform  input  electrical  energy  into  mechanical  energy. 
Of  the  many  types  of  actuator  materials  including,  magnetostrictive,  photostrictive,  and  shape 
memory  alloys,  piezoelectric  and  electrostrictive  ceramics  are  widely  used  in  applications 
requiring  high  generative  force,  high  frequency  operation,  accurate  displacement,  qmck  response 
time,  or  small  device  size  [1].  Generally,  among  the  material  properties  determining  actuator 
performance,  the  electric  field  (E-field)  induced  strain  is  the  most  important  parameter  for 
actuators.  This  is  demonstrated  by  strain  energy  density  which  is  a  measure  of  the  energy  per  unit 
mass  an  actuator  can  deliver, 

CiHox  ~  l/p  '  1/4  ‘  (1/2  •  E(SmaK)  )  ^  ^ 

where  is  the  strain  energy  density,  £  is  the  actuator’s  elastic  modulus,  Sm«  is  the 
maximum  field  induced  strain,  and  p  is  the  actuator’s  density.  [1]  1/4  is  a  factor  of  appropnate 
for  an  actuator  with  its  impedance  matching  to  that  of  its  surroundings.  In  designing  an  actuator, 
the  maximum  strain  energy  density  should  be  as  high  as  possible.  In  electroactive  ceramics, 
density  and  elastic  modulus  vary  little  from  material  to  material,  therefore,  the  level  of  strain  and 
maximum  strain  achievable  with  a  reasonable  electric  field  (<50kV/cm)  dominates  the  energy 
density.  The  piezoelectric  coefficient  (dy),  determining  the  level  of  induced  strain  at  a  given 
electric  field,  is  the  most  widely  used  parameter  describing  actuator  performance. 

It  is  the  objective  of  this  paper  to  report  ultrahigh  piezoelectric  coefficients  (das)  and  ultrahigh 
strain  levels  with  low  hysteresis  observed  for  single  crystals  of  relaxor  perovskite 
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Pb(Zni/3Nb2/3)03  -  PbTi03  (PZN-PT)  and  Pb(Mgi/3Nb2/3)03  -  PbTi03  (PMN-PT).  Strain 
behavior  as  a  function  of  E-field  will  be  discussed  with  respect  to  crystal  structure,  orientation, 
domain  configuration  and  anticipated  actuator  performance. 

2.  Background 

2.1  Morphotropic  phase  boundary  and  electrostrictive  ceramics 

Pb(Zri.x,Tix)03  (PZT)  ceramics  have  been  the  mainstay  for  high  performance  actuator 
applications.  Compositionally,  PZT  ceramics  lie  near  the  morphotropic  phase  boundary  (MPB) 
between  the  tetragonal  and  rhombohedral  phases  as  shown  in  figure  1.  MPB  compositions 
exhibit  anomalously  high  dielectric  and  piezoelectric  properties  as  a  result  of  enhanced 
polarizability  arising  from  the  coupling  between  two  equivalent  energy  states,  i.e.  the  tetragonal 
and  rhombohedral  phases,  allowing  optimum  domain  reorientation  during  the  poling  process. 
Alternate  MPB  systems  can  be  found  in  Relaxor-PbTi03 ,  also  as  shown  in  figure  1.  Lead  based 
reiaxor  matftrials  are  complex  perovskites  with  the  general  formula  Pb(B, 62)03,  (Bi— Mg  ,  Zn  , 
Ni^"^,  Sc^"^  ...,  B2=Nb®'^,  Ta*^  ...).  Characteristic  of  relaxors  is  a  broad  and  frequency 

dispersive  dielectric  maxima  [2]. 

To  achieve  a  high  piezoelectric  coefficient,  MPB-based  ceramics  are  further  engineered  by 
compositionaly  adjusting  the  Curie  temperature  (Tc)  downward  relative  to  room  temperature. 
The  effect  of  transition  temperature  (Tc)  on  the  piezoelectric  properties  is  clearly  evident  in 
figure  2.  As  shown,  the  room  temperature  values  of  d33  are  plotted  as  a  function  of  Tc  for  a 
variety  of  modified  PZT  ceramics,  including  Relaxor-PT  systems.  Enhanced  piezoelectric 
activity  of  MPB-based  ceramics,  achieved  by  compositionaly  adjusting  Tc  downward  relative  to 
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room  temperature,  results  in  ‘soft’  piezoelectric  ceramics.  This  enhanced  piezoelectric  effect, 
therefore,  comes  with  the  expense  of  more  temperature  dependent  properties,  and  less 
polarization  stability,  i.e.  aging  and  loss  of  piezoelectric  activity.  Further  details  on  the 
relationship  between  dielectric/piezoelectric  properties  and  Curie  temperature  (Tc)  of 
piezoelectric  ceramics  can  be  found  in  the  article  by  S.-E.  Park  et  al.  [3] 

Most  importantly,  a  consequence  of  increased  piezoelectric  activity  for  these  “soft”  ceramics 
is  large  hysteresis  in  the  strain  vs.  E-field  behavior  as  a  result  of  domain  motion.  Strain  vs.  E- 
field  behavior  for  PZT-5H  (Navy  type  -VI)  is  shown  in  figure  3  as  an  example.  Though  the 
piezoelectric  coefficient  (das)  of  PZT-5H  ceramics  is  in  the  range  of  -600  to  700  pC/N  [4] 
(implying  -0.06  to  0.07%  strain  at  10  kV/cm),  strain  as  high  as  0.1%  can  be  observed  at 
lOkV/cm.  This  enhanced  nonlinear  strain  is  the  result  of  domain  motion,  and  therefore, 
accompanied  by  significant  hysteresis,  resulting  in  poor  positioning  accuracy.  The  area  within  the 
strain  vs.  E-field  curve,  or  dielectric  loss,  also  results  in  significant  heat  generation  during 
operation.  Heat  generation  combined  with  a  decreased  temperature  usage  range,  results  in  poor 
temperature  stability  and  limits  these  ceramics  to  low  frequency  applications. 

Strain  vs.  E-field  hysteresis  can  be  minimized  with  the  use  of  the  hard  piezoelectric 
ceramics.  Hard  piezoelectric  ceramics  such  as  PZT-8  (Navy  type  HI)  offer  very  low  hysteresis  as 
shown  in  figure  3.  However,  the  reduction  in  hysteresis  and  loss  comes  at  the  expense  of  das  and 
subsequent  strain  level.  Typically  daa  values,  for  ‘hard’  PZTs  range  from  -200  to  300  pC/N  [4]. 

Another  category  of  ceramic  materials  used  in  commercial  actuators  are  electrostrictors. 
Electrostrictive  strain  is  proportional  to  the  square  of  polarization.  A  few  matenals  such  as  PMN 
and  its  solid  solution  with  PT  exhibit  significant  electrostrictive  strain  (>  0.1%)  with  virtually  no 
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hysteresis  as  shown  in  figure  3.  Effective  das ’s  >  -800  pC/N  calculated  directly  from  the  strain 
vs.  E-field  curve  can  be  achieved,  but  only  over  a  very  narrow  range  of  E-field  and  temperature. 
For  hard  piezoelectric  and  electrostrictive  ceramics,  strain  level  with  low  hysteresis  does  not 
exceed  0.15%.  This  limitation  originates  from  the  material’s  dielectric  breakdown  strength  and 
polarization  saturation. 

In  summary,  piezoelectric  and  electrostrictive  ceramics  offer  strain  levels  up  to  ~  0.15%.  Soft 
PZT’s,  exhibiting  piezoelectric  coefficients  (das)  as  high  as  750  pC/N,  are  inherently  limited  due 
to  hysteresis  caused  by  domain  motion.  Hysteresis  can  be  minimized  with  the  use  of  hard 
piezoelectric  ceramics,  but  daa  values  of  only  -200  to  300  pC/N  are  available.  Even  though 
electrostrictive  ceramics  offer  effective  daa ’s-  800  pC/N,  maximum  strain  level  is  limited  by  its 
dielectric  breakdown  strength  and  polarization  saturation.  To  achieve  E-field  induced  strain 
levels  >  0.15%,  electroactive  materials  should  possess  high  piezoelectric  coefficients  (das 
>1000pC/N)  and  high  dielectric  breakdown  strength. 

2.2  Single  crystal  piezoelectrics 

Single  crystal  piezoelectrics  such  as  quartz  (SiOa),  lithium  niobate  (LiNbOa),  and  the  analogue 
lithium  tantalate  (LiTaOa)  are  widely  employed  in  specific  applications  that  include  oscillators, 
surface  acoustic  wave  (SAW)  devices,  and  in  optics.  In  contrast  to  PZT  ceramics,  however, 
these  single  crystals  offer  inferior  piezoelectric  properties,  with  daa  ‘s<  50  pC/N. 

Attempts  to  grow  single  crystals  of  MPB  PZTs  have  been  made  by  numerous  researchers, 
resulting  in  crystallites  too  small  to  allow  adequate  property  measurements  [5, 6, 7,8,9].  In 
contrast  to  PZT  crystal  growth,  relaxor-PT  materials  can  be  readily  grown  in  single  crystal  form. 
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This  key  distinction  was  first  realized  by  Nomura  and  co-workers  for  the  PZN  and  PZN-PT 
systems  [10,11]  and  later  by  Shrout  for  the  PMN-PT  [12].  In  general,  most  Pb(B,,B2)03-PT 
crystals  can  be  grown  by  high  temperature  solution  growth  using  Pb-based  fluxes.[13,14] 

Piezoelectric  coefficients  as  high  as  ~  1500  pC/N  have  been  reported  [11,12]  for  MPB 
Relaxor-PT  crystals.  However,  it  should  be  noted  that  piezoelectric  coefficients  are  generally 
determined  using  low-field  (<  O.lkV/cm)  techniques  such  as  the  resonance  method  (IEEE 
standard)  [15].  Therefore,  direct  observation  of  the  strain  vs.  E-field  behavior  is  essential  m  order 
to  investigate  hysteresis  and  maximum  levels  of  strain,  key  experiments  to  directly  confirm 
actuator  performance.  In  relation  to  actuators,  several  questions  arise,  1)  How  do  the  high  das 
values  determined  using  low  field  techniques  correlate  to  direct  measurements?  As  for 
piezoelectric  and  electrostrictive  ceramics,  will  the  strain  level  saturate  with  increased  E-field?  2) 
How  much  hysteresis  accompanies  the  strain?,  3)  Is  morphotrophy  essential  for  enhan 
piezoelectric  properties?,  4)  Are  there  optimum  crystallographic  cuts  as  in  the  case  of  the  other 
piezoelectric  crystals?  In  the  following  sections,  we  will  attempt  to  answer  these  questions, 
reporting  piezoelectric  properties  and  direct  observation  of  strain  behavior  as  a  function  of 
crystallographic  orientation  and  electric  field  for  Relaxor-PT  single  crystals. 

Commonalties  inherent  to  Relaxor-PT  systems  have  been  discussed  in  reviews  by  Shrout  [16] 
and  Randall  [17].  Based  on  these  commonalties,  our  research  was  limited  to  two  representative 
systems,  PZN-PT  and  PMN-PT.  Though  PMN-PT  MPB  crystals  exhibit  piezoelectric  properties 
comparable  with  PZN-PT,  more  focus  was  given  to  the  PZN  system  owing  to  its  relatively  lower 
PT  content  for  MPB,  allowing  more  uniform  crystal  growth  of  these  solid  solution  materials. 
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3.  Experimental  Procedure 
3.1  Crystal  growth 

Crystals  of  PZN,  PMN  and  their  solid  solutions  with  PT  were  grown  using  the  high 
temperature  flux  technique.  High  purity  (>99.9%)  powders  of  Pb304  (Aldrich,  WI),  ZnO 
(Johnson  Matthey,  MA),  MgCOa  (Johnson  Matthey,  MA),  NbaOs  (Aldrich,  WI),  and  T1O2 
(Aldrich,  WI)  were  used.  Raw  powders  were  weighed  with  desired  molar  ratio  with  excess  Pb304 
as  a  flux.  The  powders  were  dry  mixed  for  a  desired  period  of  time  using  a  tumbling  mill.  The 
mixed  powders  were  loaded  into  a  Platinum  crucible,  which  was  placed  in  an  alumina  crucible 
sealed  with  an  alumina  lid  and  alumina  cement  to  minimize  PbO  volatilization.  The  crucible  and 
powder  were  placed  in  a  tube  furnace  and  held  at  soak  temperatures  (1100  to  1200“C),  followed 
by  slow  cooling  (1  to  5“C/hr).  The  crucible  was  then  furnace-cooled  to  room  temperature.  Hot 
HNO3  is  used  to  separate  the  crystals  out  of  the  rest  of  the  melt.  Typically  crystal  size  ranged 
from  3  to  20  mm.  Further  details  on  the  flux  growth  technique  of  these  crystals  can  be  found  in 

ref.  [13,14]. 

2.2  X-ray  diffraction  (XRD)  analysis 

Single  crystals  were  ground  into  a  fine  powder  for  XRD  analysis  ranging  from  20  -  80  20 
with  a  step  size  of  0.01  and  a  counting  time  of  3  sec,  in  order  to  confirm  phase  pure  perovskite 
and  to  calculate  lattice  parameter.  Individual  crystals  were  oriented  along  various 
crystallographic  directions  such  as  their  pseudocubic  <001>  and  the  <11 1>  axis  using  a  Laue 
back  reflection  camera. 
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2.3  Electrical  characterization 


Dielectric  and  piezoelectric  properties  were  measured  using  direct  observation  of  strain  as  a 
function  of  electric  field  as  well  as  low-field  property  measurements  using  the  IEEE  resonance 
technique.  [15]  For  electrical  characterization,  samples  were  prepared  by  polishing  with  silicon 
carbide  and  alumina  polishing  powders  to  achieve  flat  and  parallel  surfaces  onto  which  gold 
electrodes  were  sputtered.  High-field  measurements  included  polarization  and  strain  hysteresis 
a  modified  Sawyer-Tower  circuit  and  linear  variable  differential  transducer  (LVDT)  driven 
by  a  lock-in  amplifier  (Stanford  Research  Systems,  Model  SR830).  Plate  shape  samples  with 
thickness  ranging  fi-om  0.2  mm  to  0.5  mm  were  used.  Electric  fields  as  high  as  -140  kV/cm  were 
applied  using  an  amplified  unipolar  wave  form  at  0.2  Hz,  using  a  Trek  609C-6  high  voltage  DC 
amplifier.  During  testing  the  samples  were  submerged  in  Fluorinert  (FC-40,  3M,  St.  Paul,  MN). 
an  insulating  liquid,  to  prevent  arcing.  For  piezoelectric  coefFicient(d33)  determination,  bar  shape 
samples  with  lengths  ranging  from  3  mm  to  5  mm  were  tested.  Samples  were  poled  either  by 
field  cooling  (lOkV/cm)  from  temperatures  above  the  dielectric  maximum  temperature  (Tmax)  or 
by  applying  40kV/cm  at  room  temperature.  Figure  4  presents,  from  left  to  right,  a  representative 
of  an  as  grown  crystal,  an  aligned  sample  using  the  Laue  camera,  a  cut  and  polished  sample,  an 
electroded  sample  for  direct  strain  observation  as  well  as  for  low  field  planar  mode  resonant 
sample,  and  a  low  field  longitudinal  mode  resonant  sample,  for  PZN-9.5%PT,  respectively. 

4.  Results  and  Discussion 
4.1  Low  field  measurements 
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Piezoelectric  coefficients  as  a  function  of  composition  and  crystal  orientation  for  PZN-PT, 
calculated  based  on  IEEE  standards,  are  presented  in  figure  5.  As  shown,  large  piezoelectric 
coefficients  (dja  ~  1600  pC/N)  were  found  for  PZN-PT  with  MPB  compositions  (9.5%  PT),  as 
previously  reported  by  Kuwata  [11].  PMN-PT  crystals  with  MPB  compositions  (PMN-35%  PT) 
also  exhibited  large  piezoelectric  coefficients  (dsa  ~  1500  pC/N).[3]  It  should  be  noted,  however, 
that  all  rhombohedral  crystals  oriented  along  their  pseudocubic  <001>  direction  exhibited  large 
piezoelectric  coefficients.  As  shown  in  figure  5,  das  increased  with  increased  amount  of  PbTiOs 
for  <001>  oriented  rhombohedral  crystal.  Maximum  das  values  of -2500  pC/N  were  determined 
with  domain  engineered  rhombohedral  PZN-8%PT  crystals  .  In  contrast  to  PZT  s,  daa 
dramatically  decreases  at  MPB  to  levels  -  500pC/N  for  tetragonal  composition.  Though  <1 1 1>  is 
the  polar  direction  for  rhombohedral  crystals,  such  cuts  exhibited  low  piezoelectric  coefficient  as 
shown  in  figure  5.  In  the  following  section  the  origin  of  apparent  anisotropy  will  be  discussed  by 
direct  strain  observation  as  a  function  of  electric  field. 

4.2  Domain  (in)stability  and  associated  anisotropy 

Table  1  presents  dielectric  and  piezoelectric  properties  for  two  rhombohedral  crystals,  pure 
PZN  and  PZN-8%PT,  as  a  function  of  crystallographic  orientation.  Dielectric  loss  values  <  1% 
and  electromechanical  coupling  constant  values  (kss)  >  90%,  attractive  for  transducer 
applications  and  essential  for  high  performance  solid  state  actuators,  could  be  achieved  tvith 
<001>  orirated  crystals,  with  further  details  found  in  ref.  [3].  Oriented,  <001>  rhombohedral 


*  Rhombohedral  crystals  oriented  and  poled  along  pseudocubic  <001>  direction.  Ciystallographically,  polar 
directio  of  ihombohedral  crystal  is  pseudocubic  <1 1 1>  duection. 
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crystals  exhibited  ulttahigh  piezoelectric  coefficients  (ds,)  of-1  lOOpC/N  (PZN)  up  to  2500  pON 
(PZN-8%PT),  and  djs  values  increased  with  increased  PbTiOs  content  as  presented  in  figure  5. 
However,  eleetromeehanical  couplings  (k„)  and  piezoelectric  coefficients  (ds,)  of  <11 1> 
oriented  rhombohedral  crystals  are  only  -35%  and  ~80pC/N,  respectively,  regardless  of 
......roririnn  This  composition-independence  of  piezoelectric  properties  must  have  the  origins 

Other  than  the  inherent  crystallographic  anisotropy. 

Polarization  and  strain  as  a  function  of  E-field  (bipolar)  for  <1  M>  oriented  PZN  crystal  are 
presented  in  figures  6a  and  6b,  respectively.  Subsequem  unipolar  strain  behavior  after  poling 
(field  cooling  under  20kV/cm  from  200°C  to  room  temperature)  is  shown  in  figure  7a.  As  <1 1 1> 
is  polar  direction,  compiete  poling  results  in  a  single  domain  state  and  piezoelectric  strain 
behavior  is  expeeted  to  be  hysteresis  free  for  <1 1 1>  poled  PZN  crystals.  However,  as  shown  in 
figure  7a,  significant  hysteresis  together  with  abnormally  high  strain  values  (0.1%  at  lOkV/cm) 
indicates  domain  reorientation  under  bias,  as  diseussed  in  section  2.1.  Domain  motion  along  with 
increased  strain  must  include  109*  or  71”  domain  reorientation,  starting  at  E  -  2kV/cm  m  this 
case  (figure  7a),  followed  by  strain  saturation.  It  should  be  noted  that  remnant  strain  values  of 
-0.04%  after  E-field  is  removed  and.  original  sample  dimension  recovered  its  zero  strain  point  a 
few  seconds  after  E-fieid  is  removed,  as  presented  in  figure  7a.  This  behavior  obviously  refleets 
domain  motion,  including  depoling  as  well  as  reorientation.  As  in  case  of  “soft”  ceramics,  this 
a.'m.in  motion  is  the  origin  of  abnormally  high  strain  values  and  hysteresis.  Piezoeleetnc 
coeffieient  (drs)  values  -  SOpCTN  determined  by  low-field  resonance  technique  could  be 
explained  by  the  slope  of  strain  vs.  E-field  at  low  field  region,  as  given  in  figure  7a.  Various  d„*s 
calculated  by  the  slop  of  strain  vs.  E-field  are  given,  too.  Therefore,  the  inferior  piezoelectrie 
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properties  of  <111>  oriented  rhombohedral  crystals  detected  by  low-field  resonance  method  are 
strongly  related  to  domain  instability.  Complete  domain  orientation  and  single  domam 
configuration  under  bias  may  cause  elastic  energy  in  crystals  to  increase,  however,  resulting  in 
depoling  after  removing  E-field. 

In  contrast,  the  domain  configuration  of  <001>  oriented  rhombohedral  crystals  was  found  to 
be  stable.  Polarization  and  strain  vs.  E-field  (bipolar)  curves  of  <001>  oriented  PZN  crystals  are 
presented  in  figure  6c  and  6d,  respectively.  Domain  switching  by  bipolar  field  abruptly  occurs  at 
Ec  (~3kV/cm),  similar  behavior  found  with  single  domain  ferroelectric  crystals.  Schematics  of 
domain  configurations  for  <00l>  poled  rhombohedral  crystals  are  shown  in  figure  8.  Since 
<1 1 1>  is  the  polar  direction,  <001>  poled  crystals  must  have  the  configuration  that  each  domain 
has  one  of  four  possible  polar  directions  -  <11 1>,  <-11 1>,  <1-1 1>, and  <-1-1 1>.  Remnant 
polarization  (Pr  ~  25  nC/cm^)  of  <001>  oriented  rhombohedral  PZN  crystal  must  be  iWS  of 
<1 1 1>  oriented  crystal’s  according  to  figure  8,  which  explains  ~  25  pC/cm  for  Pr.<ooi>  out  of 
Pr,<in>  ~  43  pC/cm^  .  Therefore,  abrupt  domain  switching  at  Ec  (~3kV/cm)  is  expected  to  be  71 
domain  switching,  as  shown  in  figure  9.  In  figure  7b,  <001>  poled  PZN  crystals  exhibit  almost 
hysteresis-free  strain  behavior,  a  consequence  of  engineered-domain  stability.  From  the  slope  of 
strain  vs.  E-field  curve  revealed  dss  value  of  — 1  lOOpC/N,  as  determined  by  low-field  resonance 
technique.  This  behavior  for  <00 1>  oriented  rhombohedral  crystals  will  be  discussed  in  the 

following  section. 

The  observed  domain  (in)stability  as  a  fimction  of  crystallographic  orientation  holds  for  all 
rhombohedral  crystals  of  PZN-PT  and  PMN-PT.  Owing  to  low  hysteresis  as  well  as  ultrahigh 
piezoelectric  coefficient  (dss),  <001>  oriented  rhombohedral  crystal  may  be  attractive  candidates 
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for  actuators.  Although  domain  (in)stability  is  clearly  a  fimction  of  crystallographic  orientation, 
in-situ  domain  observation  using  optical  microscopy  is  necessary  and  vdll  be  reported  elsewhere. 

[19] 

4.3  Optimum  crystallographic  orientation 

As  discussed  in  section  2.2,  optimum  crystallographic  cuts  was  one  of  the  merits  umquely 
utilized  with  single  crystals.  For  conventional  piezoelectric  single  crystals,  however,  an 
“optimum  cut”  was  determined  from  a  single  domain  crystal,  whereas  engineered  multidomain 
state  is  being  utilized  in  this  study.  The  term  “optimum  orientation  ’  includes  crystallographic 
direction  not  only  for  poling  but  for  actuator  driving  to  obtain  hysteresis  (strain)  values  as  low 
(high)  as  possible,  respectively. 

Though  it  was  determined  that  <001>  is  superior  to  <11 1>,  multidomain  ihombohedral 
crystals  may  allow  other  crystallographic  directions  as  optimum  for  higher  strain  values  with 
comparably  low  hysteresis.  To  investigate  this  possibility,  piezoelectric  samples  oriented  along 
<001>  +  a  were  made,  where  a  is  the  degree  of  deviation  from  <001>  toward  <1 1 1>,  as  shown 
in  figure  9,  using  PZN-4.5%PT  crystal.  In  figure  10  remnant  polarization  (Pr)  approaches  the  Pr 
<,„>  value  (~43  pC/cm")  as  orientation  (a)  approaches  <1 1 1>  (54.7“),  respectively.  Ec  value  also 
increase  with  increased  a.  Figure  11  presents  unipolar  strain  behavior  as  a  function  of  a  for 
<001>  poled  PZN-4.5%PT  crystals.  At  low  field,  strain  values  do  not  vary  significantly  as  a 
function  of  a,  attributed  to  strain  induced  by  domain  reorientation  when  a^O.  For  example,  the 
strain  value  of  a  sample  (a=10®)  is  higher  than  that  for  a  <001>  oriented  sample  (a=0®).  This 
increased  domain  motion  ivith  increased  a  is  presented  in  figure  12,  with  hysteresis  vs.  a 
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calculated  from  the  area  of  polarization  vs.  unipolar  E-field  curves  (20kV/cm).  The  amount  of 
hysteresis  can  be  translated  into  degree  of  depoling,  because  more  depoling  consequently 
involves  more  domain  reorientation,  resulting  in  increased  hysteresis.  Therefore,  domain 
becomes  more  unstable  with  the  increased  a.  From  these  observations,  the  closer  to  <001>  is 
riiombohedral  crystal  orientation,  the  more  stable  domain  configuration  and  the  higher  strain 
value  are  achieved.  Therefore  the  optimum  crystallographic  orientation  is  <001> . 

High-field  strain  behavior  is  shown  in  figure  13  for  <001>  oriented  PZN-4.5%PT  crystals. 
Strain  saturation  becomes  more  significant  as  a  increases  (orientation  approaches  <1 1 1>).  At 
high  fields,  strain  level  decreases  with  increased  a,  the  result  of  polarization  saturation  and 
subsequent  strain  saturation.  However  crystals  oriented  close  to  <001>  (a<20®)  exhibit  no 
saturation,  and  furthermore,  an  apparent  jump  in  strain  to  values  over  1%.  The  question  anses, 
does  this  strain  jump  originated  from  domain  reorientation?  This  question  is  addressed  in  the 

following  section. 

4.4  Strain  vs.  E-field  behaviors  for  <00 1>  oriented  rhombohedral  crystals 

Strain  as  a  function  of  electric  field  for  various  <001>  oriented  rhombohedral  crystals  (pure 
PZN,  PZN-4.5%PT,  PZN-8%PT,  and  PMN-24%PT)  are  presented  in  figure  14.  Also,  E-field 
induced  strains  of  various  electromechanical  ceramics  such  as  soft  PZT  (PZT-5H),  hard  PZT 
(PZT-8)  and  electrostrictive  ceramics  (PMN-PT)  are  compared.  Piezoelectric  coefficients  (dss) 
directly  calculated  from  the  slope  of  strain  vs.  electric  field  curves  confirmed  the  piezoelectric 
coefficients  determined  by  the  low-field  resonance  method.  Strains  as  high  as  0.6%  were 
observed  with  low  hysteresis  for  these  crystals,  significantly  larger  than  that  for  polycrystallme 
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ceramics.  The  limitation  of  achievable  strain  for  polycrystalline  ceramics  is  the  result  of  low 
piezoelectric  coefficient,  polarization  saturation  and  subsequent  saturation  on  strain  and 
breakdown  strength. 

Figure  15  schematically  presents  engineered  domain  state  and  their  piezoelectric  response 
imder  bias  for  rhombohedral  crystals  oriented  and  poled  along  <00 1>.  When  actuation  is  induced 
by  an  E-field  along  <00 1>,  the  polar  direction  is  expected  to  incline  close  to  E-field  direction  in 
each  domain  (step  A  in  figure  1 5),  possibly  resulting  in  increased  rhombohedral  lattice  distortion. 
Domain  reorientation  is  not  necessary  during  this  step  because  neighboring  domains  must 
involve  equal  amount  of  induced  distortion,  nullifying  all  strains  in  individual  domains,  the 
reason  of  stable  domain  configuration  and,  therefore,  the  origin  of  hysteresis-free  or  low- 
hysteresis  strain  behavior  for  <00 1>  oriented  rhombohedral  crystals. 

Strain  saturation  for  <00 1>  oriented  rhombohedral  crystals  was  investigated  until  dielectric 
breakdown,  as  presented  in  figure  16.  Far  from  saturation,  the  strain  abruptly  increased  with 
strain  levels  as  high  as  0.8%  being  achieved  for  all  crystals  tested.  Strain  levels  >1.2%  for  the 
P2^-4.5%PT  crystal  were  the  result  of  higher  breakdown  voltage  (compare  with  figure  13).  The 
observed  strain  behavior  is  believed  to  be  associated  with  an  E-field  induced  rhombohedral- 
tetragonal  phase  transition.  Polarization  inclination  towards  <00 1>  finally  result  in  collapse  of  all 
polarization's  into  <001>  direction,  step  B  in  figure  17.  This  E-field  induced  phase  transition  is 
more  apparent  in  figure  17,  presenting  the  E-field  induced  strain  behavior  of  the  <00 1>  oriented 
PZN-8%PT  crystal.  At  ~120kV/cm,  the  induced  strain  along  with  high  dielectric  breakdown 
strength  resulted  in  a  strain  levels  as  high  as  1 .7%.  The  piezoelectric  coefficient  (das)  ~  480  pC/N 
calculated  directly  fixim  the  slope  of  strain  vs.  E-field  in  the  high  field  region  (between  40  and 
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120  kV/cm  in  figure  17)  corresponded  to  values  determined  for  tetragonal  phase  with  a  dss  value 
(-500  pC/N,  in  figure  5).  Individual  steps  of  piezoelectricity  (step  A)  and  phase  transition  (step 
B)  shown  in  figure  15  are  given  in  figure  17.  Although  phase  transition  is  a  likely  explanation  for 
the  ultrahigh  strain  level,  direct  observation  of  phase  transition  using  in-situ  XRD  is  required. 

5.  Conclusion 

Pseudocubic  <00 1>  oriented  relaxor  based  rhombohedral  crystals  such  as  (l-x)PZN-xPT 
(x<0.09)  and  (1-x)  PMN-  xPT  (x<0.35)  exhibited  actuation  levels  not  available  with  current 
piezoelectric  ceramics.  Ultrahigh  piezoelectric  coefficients  (d33)  >  ~2500  pC/N  and  strain  levels 
up  to  0.6%  with  low  hysteresis  were  observed.  Optimum  crystallographic  orientation  was  found 
to  be  pseudocubic  <001>  for  rhombohedral  crystals,  exhibiting  high-strain  and  low-hysteresis 
behavior.  Domain  instability  could  explain  the  composition-independently  inferior  piezoelectric 
properties  for  <1 1 1>  oriented  rhombohedral  crystals.  Ultrahigh  strain  levels  up  to  1.7%  could  be 
achieved  for  <001>  oriented  rhombohedral  single  crystals  as  a  result  of  E-field  induced 
rhombohedral-tetragonal  phase  transition.  Other  relaxor  based  rhombohedral  crystals  are 
expected  to  exhibit  similar  strain  vs.  E-field  behavior.  <001>  textured  polycrystalline  ceramics 
are  expected  to  exhibit  high  strain  behaviors.  Though  clearly  promising  candidates  for  high 
performance  actuators,  further  investigations  in  crystal  growth  and  pre-stress  testing  are  required 
for  single  crystal  piezoelectrics  to  become  the  next  generation  material  of  actuators. 
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Figure  Caption 


Figure  1  Ternary  diagram  depicting  MPBs  in  PZT  and  Relaxor-PT  systems  for  piezoelectric 
ceramic. 

Figure  2  Piezoelectric  coefficient  (das)  as  a  fimction  of  transition  temperature  (Tc)  for 
piezoelectric  ceramics,  including  PZT.  modified  PZTs,  and  Relaxor-PT  systems.  Data  were 
compiled  from  references,  commercial  brochures  and  internal  investigations. 

Figure  3  Strain  vs.  E-field  behavior  for  various  electromechamcal  ceramics. 

Figure  4  An  as  grown  crystal,  an  aligned  sample  using  the  Laue  camera,  a  cut  and  polished 
sample,  an  electroded  sample  for  direct  strain  observation  as  well  as  for  low  field  planar  mode 
resonant  sample,  and  a  low  field  longitudinal  mode  resonant  sample,  for  PZN-9.5%PT, 
respectively  (from  left  to  right). 

Figure  5  das  as  a  fimction  of  crystal  composition  and  orientation. 

Figure  6  Polarization  and  strain  vs.  E-field  (bipolar)  curves  for  PZN  crystals  oriented  along 
<1 1 1>  (a,b)  and  <00 1>  (c,d). 

Figure  7  Strain  vs.  E-field  (unipolar)  curves  for  PZN  crystals  oriented  along  <in>  (a)  and 
<001>  (b). 

Figure  8  Schematics  of  engineered  domain  configuration  for  <00 1>  poled  rhombohedral  crystals. 
Figure  9  Schematic  diagram  of  sample  preparation  for  investigating  optimum  crystallographic 

orientation. 

Figure  10  Ferroelectric  hysteresis  for  PZN-4.5%PT  crystals  oriented  along  <001>  +  a,  where  a  is 
the  degree  of  deviation  from  <001>  toward  <1 1 1>,  (a)  a=0®,  (b)  a=10®,  (c)  a=20®,  (d)  a-30  , 
(e)  a=40®. 
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Figure  1 1  Strain  vs.  E-field  (unipolar)  PZN-4.5%PT  crystals  oriented  along  <001>  +  a,  where  a 
is  the  degree  of  deviation  from  <001>  toward  <11 1>,  (a)  a=0®,  (b)  a=10®,  (c)  a=20®,  (d) 
a=30®,  (e)  a=40°. 

Figure  12  Hysteresis  as  a  function  of  a,  where  a  is  the  degree  of  deviation  from  <00 1>  toward 
<11 1>.  Hysteresis  values  are  calculated  from  the  area  of  umpolar  polarization  curve  when 
maximum  E-field  is  20  kV/cm. 

Figure  13  Strain  vs.  E-field  (unipolar)  PZN-4.5%PT  crystals  oriented  along  <001>  +  a,  where  a 
is  the  degree  of  deviation  from  <001>  toward  <11 1>.  Maximum  field  is  limited  either  by 

dielectric  breakdown  or  by  voltage  limit  of  the  apparatus. 

Figure  14  Strain  vs.  E-field  behavior  for  <00 1>  oriented  rhombohedral  crystals  of  PZN-PT  and 
PMN-PT  and  for  various  electromechanical  ceramics. 

Figure  15  Schematic  diagram  of  domain  configurations  in  <00 1>  oriented  rhombohedral  crystals 
under  bias  (step  A  -  piezoelectricity,  step  B  -  induced  phase  transition). 

Figure  16  Strain  vs.  E-field  behavior  for  <001>  oriented  rhombohedral  crystals  of  PZN-PT  and 
PMN-PT  and  for  various  electromechanical  ceramics.  Maximum  field  was  limited  either  by 

dielectric  breakdown  or  by  voltage  limit  of  the  apparatus. 

Figure  17  Strain  vs.  E-field  behavior  for  <001>  oriented  PZN-8%PT  crystal.  Maximum  field  is 

limited  by  voltage  limit  of  the  apparatus. 
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Table  1  Dielecuic  and  piezoelecttic  properties  as  a  fl^on 

ihombohedral  PKZn.oNbwPs  (PZN)  and  PKZnmNbMPr  -  8 /.PbTiOs  (PZN  T). - 


Crystal 

Orientation 


Compo-  mode 
sition 


coup-  S33^(10*  Dielectric  Loss  d33(pC/N) 

ling  ‘V/N)  Constant _ 


*  Samples  were  poled  by  field  cooliiiB  from  200»C  with  applying  20kV/cm. 


Relaxor 


Pb(B,, 82)03 


PbZr03  MPB-I  PbTi03 


S.-E.  Park  and  T.  R.  Shrout 
Figure  1.  Ternary  diagram  depicting  MPBs  in  PZT 
and  Relaxor-PT  systems  for  piezoelectric  ceramic  [20]. 
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Figure  3  Strain  vs.  E-field  behaviors  for  various  electromechanical  ceramics 
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I'igurc  4  An  as  grown  crystal,  an  aligned  sample  using  the  Laue  camera,  a  cut  and 
polished  sample,  an  electroded  sample  for  direct  strain  obser\ation  as  well  as  tor  low- 
field  planar  mode  resonant  sample,  and  a  low  field  longitudinal  mode  resonant  sample, 
for  PZN-9.5“oP  r.  respecti\  ely  ( from  left  to  right). 
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Figure  5  d,,  as  a  function  of  crystal  composition  and  orientation 


1 1 1  Orientation 


001  Orientation 


Electric  Field  (kV/cm)  Electric  Field  (kV/cm) 


S,-E.  Park  and  T.  R.  Shrout 

Rgure  6  Polarization  and  strain  vs.  E-field  curves 

for  PZN  crystals  oriented  along  <11 1>  (a,b)  and  <001  >  (c,d). 
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Figure  7  Polarization  and  strain  vs.  E-field  curves  (unipolar) 

for  PZN  crystals  oriented  and  poled  along  <1 1 1>  (a)  and  <001  >  (b). 
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Figure  8  Schematics  of  engineered  domain  configuration 
for<001>  poled  rhombohedral  crystals 
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Figure  9  Schematic  diagram  of  sample  preparation 

for  investigating  optimum  crystallographic  orientation 
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Figure  1 1  Strain  vs.  E-fieid  (unipolar)  PZN-4.5%PT  crystals  oriented  along  <001  > 
where  a  is  the  degree  of  deviation  from  <001  >  toward  <1 1 1  >. 

(a)  a=0*,  (b)  a=10*,  (C)  a=20",  (d)  a=30“,  (e)  a=40*. 
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Figure  12  Hysteresis  as  a  function  of  a,  where  a  is  the 
degree  of  deviation  from  <001  >  toward  <1 1 1>.  Hysteresis  values 
are  calculated  from  the  area  of  unipolar  polarization  curve 
when  maximum  E-field  Is  20  kV/cm. 
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Figure  13  Strain  vs.  E-field  (unipolar)  PZN-4.5%PT  crystals  oriented  along  <001  a. 

where  a  is  the  degree  of  deviation  from  <001>  toward  <1 1 1>. 

Maximum  field  Is  limited  either  by  dielectric  breakdown  or 

by  voltage  limit  of  the  apparatus. 
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Figure  14  Strain  vs.  E-field  behaviors  for  crystals  of  PZN-PT  and  PMN-PT, 
and  for  various  electromechanical  ceramics. 
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Figure  15  Schematic  diagram  of  domain  configurations 
rhombohedral  crystals 
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Figure  16  Strain  vs.  E-field  behaviors  for  crystals  of  PZN-PT  and  PMN-PT, 
and  for  various  electromechanical  ceramics.  Maximum  field  was  limited 
either  by  dielectric  breakdown  or  by  voltage  limit  of  the  apparatus. 
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Figure  17  Strain  vs.  E-field  behavior  for  <001  >  oriented  PZN-8%PT  crystal. 
Maximum  field  is  limited  by  voltage  limit  of  the  apparatus. 
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Multilayer  ceramic  actuators 

Kenji  Uchino*  and  Sadayuki  Takahashit 


Figure  1 


Examples  of  multilayer,  multimorph  and  multilayer-moonie  structures. 


Piezoelectric  and  electrostrictive  ceramic  materials  have  now 
become  key  components  in  smart  actuator/sensor  systems 
used  as  precision  positioners,  miniature  ultrasonic  motors 
and  adaptive  mechanical  dampers,  in  particular,  multilayer 
structures  have  been  intensively  investigated  in  order  to 
improve  their  reliability  and  to  expand  their  applications. 

There  have  been  recent  developments  in  the  USA,  Japan  and 
Europe. 
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Introduction 

Piezoelectric  actuators  arc  forming  a  new  field  midway 
between  electronic  and  structural  ceramics  (l**,2-4]. 
Application  fields  can  be  classified  into  three  cate¬ 
gories:  positioners,  motors,  and  vibration  suppressors.  The 
manufacturing  precision  of  optical  instruments  such  as 
lasers  and  cameras,  and  the  positioning  accuracy  for 
fabricating  semiconductor  chips,  which  arc  adjusted  using 
solid-state  actuators,  is  of  the  order  of  0,1  mm.  Tiny  motors 
smaller  than  1  cm^  arc  often  required  in  office  or  factory 
automation  equipment,  and  conventional  electromagnetic 
motors  which  have  sufficient  energy  efficiency  are  rather 
difficult  to  produce.  Ultrasonic  motors  whose  efficiency 
is  insensitive  to  size  arc  superior  in  the  mini-motor  area. 
Vibration  suppression  in  space  structures  and  military 
vehicles  using  piezoelectric  actuators  is  also  a  promising 
technology. 

Multilayer  structures  arc  mainly  used  for  practical  appli¬ 
cations  because  of  their  low  drive  voltage,  high  energy 
density,  quick  response  and  long  lifetime  [5,6*],  Figure  I 
illustrates  multilayer,  multimorph  and  multilayer-moonie 
structures,  which  will  be  covered  in  this  review.  Recent 
investigations  have  focused  on  the  improvement  of 
reliability  and  durability  in  multilayer  actuators. 

This  article  reviews  the  investigations  of  device  structures, 
reliability  issues,  and  recent  applications  of  multilayer 
actuators,  comparing  developments  in  the  USA,  Japan  and 
Europe. 

Multilayer  structures 

Two  preparation  processes  arc  possible  for  multilayer 
ceramic  devices:  a  cut-and-bond  method  or  a  tape-casting 
method.  The  tape-casting  method  requires  expensive 


fabrication  facilities  and  sophisticated  techniques,  but  is 
suitable  for  the  mass  production  of  thousands  of  devices 
per  month.  Tape-casting  also  provides  thin  dielectric 
layers,  leading  to  low  drive  voltages  of  40-100  V  [7*,8I. 

A  multilayer  actuator  with  interdigital  internal  electrodes 
has  been  developed  by  Tokin  [9,10].  In  contrast  to  the 
conventional  electrode  configuration  in  Figure  I,  line 
electrodes  arc  printed  on  piezoelectric  ceramic  green 
sheets,  and  arc  stacked  so  that  alternating  electrode  lines 
arc  displaced  by  one-half  pitch  (sec  Fig.  2).  This  actuator 
generates  motions  at  right  angles  to  the  stacking  direction 
using  the  longitudinal  piezoelectric  effect.  Long  ceramic 
actuators  up  to  74  mm  in  length  arc  manufactured. 

A  three-dimensional  positioning  actuator  with  a  stacked 
structure  has  been  proposed  by  PI  Ceramic  (Fig.  3), 
in  which  shear  strain  is  utilized  to  generate  x  and  y 
displacements  [llj. 

Composite  actuator  structures  called  ‘moonics’  and  ‘cym¬ 
bals*  have  been  developed  at  The  Pennsylvania  State 
University  to  provide  characteristics  intermediate  between 
the  multilayer  and  bimorph  actuators.  These  transducers 
exhibit  an  order  of  magnitude  larger  displacement  than 
the  multilayer,  and  much  a  larger  generative  force  with  a 
quicker  response  than  the  bimorph  [12,13*].  The  device 
consists  of  a  thin  multilayer  piezoelectric  element  and  two 
metal  plates  with  narrow  moon-shaped  cavities  bonded 
together  as  shown  in  Figure  1.  A  moonic  5x5x2.Smm-» 
in  size  can  generate  a  20  mm  displacement  under  60 
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Figure  2 


eight  times  as  large  as  the  displacement  of  a  multilayer 
of  the  same  size  (141.  This  new  compact  actuator  has 
been  used  to  make  a  miniaturized  laser  beam  scanner  [14]. 
Moonic/cymbal  characteristics  have  been  investigated  for 
various  constituent  materials  and  sizes  [15,16]. 

Reliability  of  multilayer  actuators 

As  the  application  fields  expand,  the  reliability  and  dura¬ 
bility  issues  for  multilayer  actuators  become  increasingly 
important.  The  reliability  of  ceramic  actuators  depends 
on  a  number  of  complex  factors,  which  can  be  divided 
into  three  major  categories;  reliability  of  the  ceramic  itself, 
reliability  of  the  device  design,  and  drive  technique. 

Compositional  changes  of  actuator  ceramics  and  the  effect 
of  doping  arc  the  primary  issues  used  in  stabilizing  the 
temperature  and  stress  dependence  of  the  induced  strains. 
A  multilayer  piezoactuator  for  use  at  high  temperatures 
(150*0  has  been  developed  by  Hitachi  Metal,  using 
Sb303  doped  (Pb,Sr)(Zr,T0O3  ceramics  [17].  Systematic 
data  on  the  uniaxial  stress  dependence  of  piezoelectric 
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characteristics  have  been  collected  on  various  Navy  PZT 
(Pb(Zr,Ti)03  materials  (181.  The  grain  size  and  porosity 
control  of  the  ceramics  are  also  important  in  controlling 
the  reproducibility  of  actuators  [19].  Aging  phenomena, 
especially  the  degradation  of  the  strain  response,  are,  in 
general,  strongly  dependent  on  the  applied  electric  field 
as  well  as  on  temperature,  humidity  and  mechanical  bias 
stress  [20]. 

The  design  of  the  device  strongly  affects  its  durability  and 
lifetime.  Silver  electrode  metal  tends  to  migrate  into  the 
piezoceramic  under  a  high  electric  field  in  high  humidity. 
Silvenpalladium  alloys  suppress  this  behavior  effectively. 
Resistive  coatings  of  the  device  should  also  be  taken 
into  account  [21].  To  overcome  electrode  delamination, 
improved  adhesion  can  be  achieved  by  using  a  mesh-type 
electrode  or  an  electrode  material  with  mixed  metal  and 
ceramic  (the  matrix  composition!)  powders.  Pure  ceramic 
electrode  materials  have  also  been  developed  using 
semiconductive  perovskite  oxides  (barium  titanate-based 
PTCR  [positive  temperature  coefficient  of  resistivity] 
ceramics)  [22].  The  lifetime  characteristics  of  a  multilayer 
actuator  with  applied  DC  or  unipolar  AC  voltage  at  various 
temperatures  [23, 24**]  and  at  various  humidities  [25] 
have  been  investigated  by  Nagata  and  Kinoshitao.  The 
relationship  between  the  logarithm  of  the  lifetime  and 
the  reciprocal  of  the  absolute  temperature  showed  linear 
characteristics  similar  to  the  Arrhenius  type.  Nevertheless, 
the  degradation  mechanism  remains  a  critical  problem. 

In  multilayer  actuators,  reduction  of  the  tensile  stress 
concentration  around  the  internal  electrode  edge  of 
the  conventional  interdigital  configuration  is  the  central 
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Figure  4 _ _ _ ^ 

Various  internal  electrode  configurations 
in  multilayer  actuators,  (a)  Interdigital,  (b) 
Plate-through,  (c)  Slit-insert,  and  (d)  Float 
electrode. 
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Interdigital 


(b) 

Plate-through 


(C) 

Interdigitai 

w/slit 


(d) 

Interdigitai 
w/float  electrode 


problem.  With  regard  to  the  destruction  mechanism  of 
multilayer  ceramic  actuators,  systematic  data  collection 
and  analysis  have  led  to  considerable  progress  [26-32,33*I. 
Two  typical  crack  patterns  are  generated  in  a  conventional 
interdigitai  multilayer  device;  one  is  a  Y-shaped  crack 
located  on  the  edge  of  an  internal  electrode,  and  the 
other  is  a  vertical  crack  located  in  a  layer  adjacent  to  the 
top  or  bottom  inactive  layer,  connecting  a  pair  of  internal 
electrodes. 

To  overcome  this  crack  problem,  three  electrode  config¬ 
urations  have  been  proposed  (see  Fig.  4):  plate-through, 
interdigitai  and  slit,  and  interdigitai  and  float  electrode 
types.  The  ‘float  electrode’  type  is  an  especially  promising 
design  which  can  be  fabricated  using  almost  the  same 
process  as  the  conventional  multilayer  actuator,  and  can 
lead  to  much  longer  lifetimes  (34*1.  An  empirical  rule 
‘the  thinner  the  layer,  the  tougher  the  device’  (27,28) 
is  also  very  intriguing,  and  will  receive  more  theoretical 
investigation  in  the  near  future. 

Failure  detection  or  lifetime  prediction  methods  are 
expected  to  remarkably  increase  the  reliability  of  multi¬ 
layer  actuators.  Acoustic  emission  and  surface  potential 
monitoring  are  promising  methods  (35).  The  Pennsylvania 
State  University  has  developed  a  modified  multilayer 
actuator  containing  a  strain  gauge  as  an  internal  electrode 
(36).  This  internal  strain  gauge  electrode  can  detect  the 
crack  initiation  sensitively  and  monitor  the  field-induced 
strain  [37**1. 

When  considering  drive  techniques  for  ceramic  actuators, 
the  pulse  drive  and  AC  drive  require  special  attention. 
The  vibration  overshoot  associated  with  a  sharp-rise 
step/pulse  voltage  causes  a  large  tensile  force,  leading 
to  delamination  of  the  multistacked  structure,  while 
the  long-term  application  of  an  AC  voltage  generates 
considerable  heat.  A  special  pulse  drive  technique  using 
a  mechanical  bias  stress  is  required  in  the  first  case, 
and  heat  generation  can  be  suppressed  by  changing 
the  device  design.  An  analytical  approach  to  the  heat 
generation  mechanism  in  multilayer  actuators  has  been 
reported,  indicating  the  importance  of  a  larger  surface 
area  [38].  Heat  generation  in  piezoelectric  ceramics  is 


mainly  attributed  to  the  P-E  (polarization-electric  field) 
hysteresis  loss  under  large  electric  field  drive  [39*, 40). 
For  ultrasonic  motors,  antiresonance  drive  is  preferable  to 
resonance  drive  because  of  higher  efficiency  and  lower 
heat  generation  for  the  same  vibration  level  (41). 

Applications  of  multilayer  actuators 

Table  1  compares  a  variety  of  ceramic  actuator  develop¬ 
ments  in  the  USA,  Japan  and  Europe.  Additional  details 
will  be  described  below. 

USA 

The  principal  target  is  military-oriented  applications  such 
as  vibration  suppression  in  space  structures  and  military 
vehicles.  Substantial  up-sizing  of  the  actuators  is  required 
for  these  purposes,  A  typical  example  is  found  in  the 
aircraft  wing  proposed  by  NASA  (42).  A  piezoelectric 
actuator  was  installed  near  the  support  of  the  wing, 
allowing  immediate  suppression  of  unwanted  mechanical 
vibrations.  Several  papers  have  been  reported  on  damper 
and  noise  cancellation  applications  (43,44). 

Passive  dampers  constitute  another  important  application 
of  piezoelectrics,  where  mechanical  noise  vibration  is 
radically  suppressed  by  the  converted  electric  energy 
dissipation  through  Joule  heat  when  a  suitable  resistance, 
equal  to  an  impedance  of  the  piezoelectric  element  l/(oC, 
is  connected  to  the  piezoelement  (45). 

A  widely-publicized  application  occurred  with  the  repair  of 
the  Hubble  telescope  launched  by  the  Space  Shuttle.  Mul¬ 
tilayer  PMN  (Pb(Mg,/3Nb2/3)03)  electrostrictive  actuators 
corrected  the  image  by  adjusting  the  phase  of  the  incident 
light  wave  (Fig.  5)  (46).  PMN  electrostrictors  provided 
superior  adjustment  of  the  telescope  image  because  of 
negligible  strain  hysteresis. 

Japan 

Japanese  industries  seek  to  develop  mass-consumer  prod¬ 
ucts,  especially  mini-motors  and  micropositionets,  aiming 
at  applications  such  as  office  equipment  and  cameras/video 
cameras.  Tiny  actuators  smaller  than  1  cm^  are  the  main 
focus  in  these  products. 
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Table  1 


Carimlc  actuator  developments  in  the  USA,  Japan  and  Europe. 

USA 

Japan 

Europe 

Target 

Military-oriented  products 

Mass*consumer  products 

Laboratory  equipment 

Category 

Vibration  suppressor 

Mini-motor  positioner 

Mini-motor  positioner 

Vibration  suppressor 

Application  field 

Space  structure 

Military  vehicle 

Office  equipment 

Cameras 

Precision  machines 

Automobiles 

Lab  stage-stepper 

Airplanes 

Automobiles 

Hydraulic  systems 

Actuator  size 

Up*8izing 
(30  cm) 

Down-sizing 
(1  cm) 

Intermediate  size 
(10cm) 

Major  manufacturer 

AVX/Kyocera 

Morgan  Matroc 
ttek  Opt.  Systems 

Burleigh 

AlliedSignal 

Tokin  Corp. 

NEC 

Hitachi  Metal 

Mitsui-Sekka 

Canon 

Seiko  Instruments 

Philips 

Siemens 

Hoechst  Ceram  Tec. 

Ferroperm 

Physik  Instrumente 

Figure  S _ _ _ 

Articulating  fold  mirror  using 
Pb{Mgi/3Nb2/3)03  multilayer  actuators. 


A  dot  matrix  printer  was  the  first  widely-commercialized 
product  using  multilayer  ceramic  actuators  [47].  Each 
character  formed  by  such  a  printer  was  composed  of  a 
24  x  24  dot  matrix  in  which  a  printing  ribbon  was  impacted 
by  a  multiwire  array.  The  printing  element  was  composed 
of  a  multilayer  piezoelectric  device  with  a  sophisticated 
hinge  lever  magnification  mechanism.  The  magnification 
by  a  factor  of  30  resulted  in  an  amplified  displacement 
of  0.5  mm  and  an  energy  transfer  efficiency  greater  than 
50%.  A  modified  impact  printer  head  has  been  developed 
by  Tokin,  using  new  interdigital  internal  electrode  type 
actuators,  which  have  lowered  the  production  cost  (48). 
A  color  ink-jet  printer  has  also  been  commercialized  by 
Seiko  Epson,  using  multilayer  piezoactuators  [49]. 


Automotive  applications  by  Toyota  Motor  have  recently 
been  accelerated.  Multilayer  actuators  have  been  intro¬ 
duced  to  an  electronically  controlled  suspension  [50]  and 
a  pilot  injection  system  for  diesel  engines  [51]. 

Efforts  have  been  made  to  develop  high-power  ultrasonic 
vibrators  as  replacements  for  conventional  electromagnetic 
motors.  The  ultrasonic  motor  is  characterized  by  low 
speed  and  high  torque,’  which  is  contrasted  with  the  ‘high 
speed  and  low  torque’  of  the  electromagnetic  motors. 
Since  the  invention  of  Ji-shaped  linear  motors  [52],  various 
modifications  have  been  reported  [53,54**].  The  Mitsui 
Petrochemical  model  is  of  particular  interest  because  the 
motor  body  is  composed  of  only  one  component  prepared 
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Figure  6 _ 

Monolithic  multilayer  piezoelectric  linear 
motor. 


by  a  cofiring  method  as  illustrated  in  Figure  6  [54**).  A 
maximum  speed  of  200mm/s  and  a  maximum  thrust  of 
60  gf  have  been  reported  for  this  motor. 

Camera  motors  utilize  a  traveling  elastic  wave  induced 
by  a  thin  piezoelectric  ring.  A  ring-type  slider  in  contact 
with  the  "rippled*  surface  of  the  elastic  body  bonded 
onto  the  piezoelectric  can  be  driven  in  both  directions  by 
exchanging  the  sine  and  cosine  voltage  inputs.  Another 
advantage  is  its  thin  design,  making  it  suitable  for 
installation  in  cameras  as  an  automatic  focusing  device. 
Nearly  80%  of  the  exchange  lenses  in  Canons  ‘EOS* 
camera  series  have  been  replaced  by  ultrasonic  motor 
mechanisms  (55*]. 

Intriguing  research  programs  arc  undenvay  in  Japan  on 
the  vibration  damping  of  earthquakes,  using  piezoelec¬ 
tric  actuators  [56,57).  Active  damping  of  a  multilayer 
piczoactuator  has  been  tested  using  an  actual  size  H-type 
steel  girder,  and  has  been  verified  to  be  effective  during 
earthquakes. 


Europe 

Ceramic  actuator  development  began  relatively  recently 
in  Europe  with  a  wide  range  of  research  topics.  Several 
major  manufacturers  arc  currently  focusing  on  laboratory 
equipment  products  such  as  laboratory  stages  and  steppers 
with  rather  sophisticated  structures. 

Figure  7  shows  a  walking  piezomotor  with  four  multilayer 
actuators  made  by  Philips  [58].  Two  short  actuators 
function  as  clampers  while  the  longer  two  provide  a 
proceeding  distance  in  an  inchworm  mechanism.  Physik 
instrumente  has  developed  more  complicated  two-leg 
type  walkers  [59]. 

Conclusions 

The  intensive  development  of  piezoelectric  actuators 
began  in  Japan  18  years  ago  and  then  spread  worldwide. 
\t  present,  attention  is  being  focused  on  practical  device 
applications.  This  article  has  reviewed  several  reliability 
issues  for  the  multilayer  ceramic  actuators  as  well  as  new 


Figure  7 _ 

An  inchworm  using  multilayer 
piezoelectric  actuators. 
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accuacor  structures,  and  compared  the  developments  of 
recent  applications  in  the  USA,  Japan  and  Europe. 

The  markets  in  USA  arc  chiefly  limited  to  military  and 
defense  applications,  and  it  is  difficult  to  estimate  the 
amount  of  sales.  The  current  Navy  needs  include  smart 
submarine  skins,  hydrophone  actuators,  and  prop  noise 
cancellation.  Smart  aircraft  skins  arc  an  Air  Force  objec¬ 
tive,  while  the  Army  requires  helicopter  rotor  twisting, 
acroscrvoclastic  control  and  cabin  noisc/scat  vibration 
cancellation. 

Meanwhile  in  Japan,  piezoelectric  shutters  (Minolta  Cam¬ 
era)  and  automatic  focusing  mechanisms  in  cameras 
(Canon),  dot-matrix  printers  (NEC)  and  part-feeders 
(Sanki)  arc  mass-produced  by  tens  of  thousands  of  pieces 
per  month.  During  the  commercialization,  new  designs 
and  drive-control  techniques  of  the  ceramic  actuators  have 
been  developed  mainly  over  the  past  few  years.  A  number 
of  patent  disclosures  have  been  generated  by  NEC, 
TOTO  Corporation,  Matsushita  Electric,  Brother  Industry, 
Toyota  Motors,  Tokin,  Hitachi  Metal,  and  Toshiba. 

If  wc  estimate  the  annual  sales  in  2001  (neglecting  the 
current  economic  recession  in  Japan),  ceramic  actuator 
units,  camera-related  devices  and  ultrasonic  motors  are 
expected  to  reach  $500  million,  $300  million  and  $150 
million,  respectively.  A  total  of  $10  billion  is  a  realistic 
target  for  the  overall  sales  actuator-related  products  [55*]. 

Future  research  trends  can  be  divided  in  two  ways: 
up-sizing  in  space  structures  and  down-sizing  in  office 
equipment.  Further  down-sizing  will  also  be  required  in 
medical  diagnostic  applications  such  as  blood  test  kits  and 
surgical  catheters. 

Key  words  for  the  future  of  multilayer  ceramic  actuators 
arc  ‘miniaturization’  and  ‘hybridization .  Layers  thinner 
than  10  ^im,  corresponding  to  current  multilayer  capacitor 
technology,  will  also  be  introduced  in  actuator  devices 
replacing  the  present  100  pm  sheets.  Piezoelectric  thin 
films  compatible  with  silicon  technology  arc  a  focus  in 
micro-electromechanical  systems.  Ultrasonic  rotary  motors 
as  small  as  2  mm  in  diameter  [60]  and  a  two-dimensional 
micro-optica  I -scanner  (611,  both  of  which  were  fabricated 
on  a  silicon  membrane  arc  good  examples. 


Nonuniform  configurations  or  heterostructures  of  different 
materials,  layer  thickness,  or  electrode  patterns  will 
be  adopted  for  practical  devices.  Functionally  gradient 
piezoelectric  actuators  now  being  prototvT)cd  indicate  a 
new  trend  [62], 
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Multilayer  piezoelectric  actuators  when  driven  under  high 
frequency,  generate  significant  heat,  which  influences  the 
reliability  and  other  piezoelectric  properties.  In  this  paper, 
heat  generation  in  various  types  of  multilayer  PZT-based 
actuators  was  studied.  Experimental  results  showed  that 
heat  generation  is  mainly  caused  by  ferroelectric  hysteresis 
loss  in  the  stress-free  state.  A  simplified  analytic  method 
was  established  to  evaluate  the  temperature  rise,  which  is 
useful  for  the  design  of  multilayer  and  other  high-power 
actuators. 

1.  Introduction 

Multilayer  piezoelectric  actuators  have  been  under  inten¬ 
sive  development  during  the  past  decade’-*  and  have  led 
to  some  applications  such  as  impact  printer  hcads,^  *  ultrasonic 
linear  motors,’  and  jr-y  stages.’®  Today,  low  voltages  are  the 
trend  for  electronic  systems,  cars,  and  medical  devices."  This 
trend  has  made  the  multilayer  structure  with  thin  layers  neces¬ 
sary.  Efforts  have  been  made  and  thin  layers  down  to  20  p.m 
have  been  achieved.* 

In  order  to  obtain  maximum  displacement  of  the  actuator,  an 
electric  field  of  1-3  kV/mm  is  usually  necess^.  For  some 
applications,  like  the  commercialized  high-speed  impact  printer 
head  with  a  printing  speed  of  3  kHz,*’  a  high-frequency  voltage 
is  applied  to  the  actuator.  Under  such  continuous  high  fre¬ 
quency,  and  driving  with  a  high  electric  field,  fairly  large  heat 
generation  has  been  observed, Heat  generation  influences  the 
reliability  and  other  propenies  and  may  also  limit  the  applica¬ 
tions  of  the  actuators.  Therefore,  a  comprehensive  investigation 
of  heat  generation  becomes  very  important. 

This  paper  has  studied  heat  generation  in  various  types  of 
multilayer  PZT-based  actuators.  A  simplified  analytic  method 
was  established  to  evaluate  the  temperature  rise,  which  is  very  use¬ 
ful  for  the  design  of  multilayer  and  other  high-power  actuators. 

II.  Experimental  Procedure 

Multilayer  PZT-based  ceramic  actuators  of  various  types 
(see  Table  I)  were  prepared  by  a  tape  casting  technique.  Tht 
thickness  of  the  active  layers  was  20  or  40  jxm.  After  cutting 
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and  the  application  of  the  end  termination,  the  samples  were 
poled  by  applying  a  high  voltage  at  an  elevated  temperature. 

The  heat  generation  experiment  was  carried  out  by  monitor¬ 
ing  the  temperature  rise  of  the  actuators  under  continuous  volt¬ 
age  driving.  The  actuators  were  floated  in  free  air  supported  by 
two  rigid  paper  legs.  Under  this  condition,  heat  conduction 
through  the  holder  was  minimized.  When  the  actuator  was  put 
directly  on  a  certain  surface,  the  heat  conduction  through  the 
sample  surface  was  increased.  The  former  sample  holder  was 
used  in  this  study. 

The  actuators  were  driven  by  continuous  triangular  unipolar 
voltages.  The  voltages  were  generated  by  an  NF  1940  Function 
Synthesizer,  amplified  through  an  NF  4010  High  Speed  Power 
Amplifier  (the  voltage  being  in  the  range  of  300  V,  current 
range  1.2  A),  and  then  applied  to  the  actuators. 

The  temperature  rise  of  the  actuators  was  measured  by 
attaching  a  very  fine  (0  =  0.001  in.)  chromel— alumel  thermo¬ 
couple  (CHAL-001 ,  OMEGA  Engineering  Inc.)  to  the  center  of 
the  sample  surface.  Silicone  grease  was  used  as  a  connective 
material.  The  temperature  was  read  through  a  thermometer 
(OMEGA  HHl  1 )  with  an  accuracy  of  O.rC. 

III.  Results 

(1 )  Temperature  Rise  Phenomena 

Heat  generation  for  various  size  actuators,  while  driving  under 
3  kV/mm  at  300  Hz,  is  shown  in  Fig.  1.  It  is  seen  from  the 
figure  that  the  actuator  temperature  initially  increases  exponen¬ 
tially  with  an  increase  in  the  driving  time  (0»  and  then  it  reaches 
equilibrium.  Temperature  rise  IT  is  defined  as 

A7'=7'“ro 

where  T  and  To  indicate  the  maximum  and  initial  actuator  teni- 
peratures.  The  relationship  between  AT  and  v./A,  where  is 
the  effective  volume  where  electric  field  is  applied,  and  A  is  the 
total  surface  area  of  the  sample,  is  shown  in  Fig.  2.  Temperature 
rise  AT  is  approximately  proportional  to  the  value  vJA,  because 
the  heat  generation  is  assumed  to  be  proportional  to  the  “Ve” 
defined  above,  and  the  heat  dissipation  is  proportional  to  the 
surface  area  A. 

Heat  generation  depends  on  driving  frequency  (/).  The  fre¬ 
quency  dependence  of  the  temperature  rise  AT  is  given  in  Fig.  3 
for  various  applied  electric  fields.  At  the  low  electric  field  level, 
AT  increases  almost  linearly  with  increasing  driving  frequency. 
As  the  electric  field  level  becomes  higher  (>1  kV/mm),  the 
AT  ^/curves  gradually  deviate  from  the  linear  trend. 

Heat  generation  is  also  electric  field  (£)  dependent.  The 
temperature  as  a  function  of  time  for  the  actuator  with  dimen¬ 
sions  of  7  mm  X  7  mm  X  2  mm  is  shown  in  Fig.  4.  Figure  5 
gives  the  temperature  rise  AT  as  a  function  of  applied  electric 
field  at  various  frequencies.  The  AT  is  initially  proportional  to 
£*.  But  the  curves  then  deviate  from  this  square  relationship 
under  higher  electric  field  levels. 
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Table  L  Various  Sizes  of  Multilayer  Actuators 


Actuator  type 

Dimensions 

(mm) 

Layer  thickness 
(pjn) 

Surface  area 

A  (mm*) 

Effective  volume' 

V,  (mm^) 

yJA 

(mm) 

^33 

4.5  X  3.5  X  2 

20 

63J 

18.8 

0.293 

7X7X2 

40 

154 

72.3 

0.470 

5  X  5  X  20 

40 

450 

308 

0.882 

17  X  3.5  X  1 

20 

160 

40.2 

0.251 

Multimorph 

31  X  9X0.3 

20 

582 

34.7 

0.078 

'Effective  volume  is  a  volume  that  actually  generates  heat 


Fig.  L  Heal  generation  for  various  actuators  (driven  at  3  kV/mm. 
300  Hz). 


(2)  P-E  Hysteresis  Loss 

Heat  generation  in  piezoelectric  actuators  is  considered  to  be 
caused  by  losses,  such  as  mechanical  loss  and  dielectric  loss.*^"'^ 
Since  the  major  contribution  to  dielectric  loss  is  from  the  ferro¬ 
electric  hysteresis  loss  under  high  electric  fields  (higher  than 
the  threshold  field  at  which  domain  wall  motions  start),‘^  *®  in 
this  paper,  the  ferroelectric  loss  was  studied.  With  the  well- 
known  Sawyer-Tower  circuit,  the  polarization  (P)-electric 
field  (E)  hysteresis  loops  were  observed.  The  P— E  hysteresis 
loops  of  different  size  actuators  are  almost  the  same,  because 
all  of  the  actuators  are  made  of  the  same  composition.  However, 


Fig.  3.  Temperanire  rise  vs  frequency  under  various  applied  electric 
fields  (data  from  the  actuator  with  dimensions  of  7  nun  X  7  nun  X 
7  nun). 


the  hysteresis  loops  depend  on  temperature,  frequency,  and  the 
electric  field.  At  different  sample  temperatures,  hysteresis  loops 
are  slightly  different  Figure  (5  shows  the  P-E  hysteresis  loops 
at  25®  and  98®C,  while  driving  under  3  kV/mm,  300  Hz.  The 
hysteresis  areas  represent  the  ferroelectric  loss  per  driving 
cycle,  tt.  The  ferroelectric  loss  m  as  a  function  of  sample  temper¬ 
ature  is  shown  in  Fig.  7.  It  is  seen  that  u  decreases  slightly  with 
an  increase  of  the  sample  temperature.  The  P— E  hysteresis 


Fig.  4.  Heat  generation  while  driving  under  various  electric  fields 
400  Hz  (data  from  the  actuator  with  dimensions  of  7  nun  X  7  mm 
2  nun). 
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Electric  Field  (kV/mm)  Sample  Temperature  CO 

Fig.  5.  Temperature  rise  as  a  function  of  applied  electric  fields  (data  Fig,  7.  P-E  hysteresis  loss  as  a  function  of  sample  temperatures 
from  the  actuator  with  dimensions  of  7  mm  X  7  mm  x  2  mm).  (3  kV/mm,  300  Hz  i. 


loss  also  decreases  slightly  with  increasing  frequency  (Fig.  8). 
However,  it  increases  noticeably  when  the  applied  electric  field 
is  increased  (Fig.  9).  The  hysteresis  loss  is  approximately  pro¬ 
portional  to  £■"  (Fig.  10). 

IV.  Theoretical  Analysis  and  Discussion 

According  to  the  law  of  energy  conservation,  the  rate  of 
heat  storage  in  the  actuator  resulting  from  heat  generation  and 
dissipation  effects  can  be  expressed  as  (assuming  uniform  tem¬ 
perature  distribution  in  the  sample) 

?o  - 

where  v,  p.  c  are  total  volume,  density,  and  specific  heat, 
respectively. 

As  staled  above,  heat  generation  is  considered  to  be  caused 
by  losses.  Thus,  the  rate  of  heat  generation  {qo)  in  the  actuator 
can  be  expressed  as 

<7o  =  h/v.  (2) 

where  u  is  the  loss  of  the  sample  per  driving  cycle  per  unit 
volume.  /  is  the  driving  frequency,  v,  is  the  effective  volume 
which  actually  generates  heat  (as  defined  above). 

Since  the  conduction  heat  transfer  rate  through  the  lead  wire 
in  this  experiment  is  negligibly  small  (i.e..  roughly  1/25  of  the 


convection),  the  rate  of  heat  dissipation  ( q^)  from  the  actuator 
is  the  sum  of  the  rates  of  heat  flow  by  radiation  {q,)  and 
convection 

4ou.  = 

-  aeA(  r  -  rj)  +  KAiJ  -  7,.)  (3) 

where  cr  is  the  Siefan-Boltzmann  constant  (5.67  _X  10"* 
W/(m*‘K'*)),  e  is  the  emissivity  of  the  sample,  and  is  the 
average  convective  heat  transfer  coefficient. 

Thus,  by  substituting  q^  and  from  Eqs.  (2)  and  (3), 
Eq.  (1)  becomes 

«/v.  -  .4  [aeir  -  rj)  +  KiT  -  To)]  =  (4) 

The  above  equation  can  be  written  in  the  form 

ufv,  “  Ak(T)(T  -  To)  =  vpc(dT/dt)  (5) 

where 

k{T)  =  aE(  T-  +  TiXr  +  T„)  +  K  (6) 

is  defined  as  the  overall  heat  transfer  coefficient.  The  first  term 
of  k(T).  aetr-  -  Ti)(T  +  T.,).  is  5.58  W/(m-K)  at  25'’C  and 
9.08  W/(m--K)  at  125'’C.  if  we  take  ^of  the  sample  as  0.93 
(similar  to  rough  fused  quartz),  while  h.  has  the  magnitude  of 


Electric  Field  (kV/mm)  Electric  Field  (kV/mm) 

Fig.  6.  P-£ hysteresis  loops  at  different  sample  temperanires:  (a)  25°.  (b)  98‘’C  (3  kV/nun.  300  Hz). 
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Fig.  8.  -P-*E  hysteresis  loss  as  a  function  of  frequency  (f  =  25T, 
£  =  2kV/inm). 


Fig.  10.  P-E  hvsteresis  loss  as  a  function  of  electric  fields  {7*  = 
25°C,/=  300HzV 


6-30  W/(m'‘K)  for  free  convection  in  air*’  (our  results  in 
Figs.  12  and  13  suggest  that  h,  is  close  to  20  W/(m’-K)  at  free 
convection  sjaie),  and  it  increases  when  increasing  the  forced 
convection  is  about  10-500  W/(m**K)  for  forced  convection 
in  air’*).  Therefore,  we  can  take  kiT)  as  a  constant  relatively 
insensitive  to  temperature  change.  Then,  the  solution  to  Eq.  (5) 
for  the  actuator  temperature  as  a  function  of  time  ( f )  is 


r-io- 


k(T)A 


^ -(4t<rM/vpr>»  j 


The  time  constant  in  Eq.  (4)  is 


(7) 


vpc 

'^~k{f)A 


(8) 


As  r  the  maximum  temperature  rise  in  the  actuator 

becomes 


Ar  = 


w/v. 

k(T)A 


(9) 


At  f  0.  where  there  is  no  heat  dissipation  effect,  the  initial 
rate  of  temperature  rise  is 


(10) 

\dt  pc*v  7 

'  '»-*0 

Curve  fit  to  the  experimental  data  with  the  solution  form 

r=ro  + Aid  (ID 

is  provided  in  Fig.  1 1  as  an  example.  Fitting  cur\'es  well  match 
the  experimental  data,  which  indicates  that  the  analytic  result 


adequately  describes  the  heat  generation  phenomena.  Conse¬ 
quently,  we  can  calculate  the  total  loss  u  of  the  actuator  through 
Eq.  (10).  The  calculated  results  are  shown  in  Table  II.  The 
experimental  data  of  P-E  hysteresis  losses  are  also  listed  for 
comparison.  It  is  seen  that  the  P-E  hysteresis  losses  are  nearly 
equal  to  the  total  losses.  Since  the  experiments  were  carried 


Fig.  11.  Curve  fit  with  the  equation  T  =  +  ST(\  -  e  (data 

from  the  actuator  with  dimensions  of  7  mm  X  7  mm  X  2  mm,  driven 
under  3  kV/mm,  300  Hz). 


Electric  Field  (kV/mm)  Electric  Field  (kV/mm>  Electric  Field  (kV/mm) 


Fig.  9,  P-E  hysteresis  loops  of  various  applied  electric  fields:  (a)  I,  (b)  2.  (c)  3  kV/mm  (7*  -  25 'C./ -  300  Hz). 
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Table  IL  Loss  and  Overall  Heat  Transfer  Coefficient  (£  =  3  kV/mm,/  =  300  Hz) 


Actuator 

4.5  X  3.5  X  2  mm 

7  mm  X  7  mm  X  2  mm 

1 7  mm  X  3.5  mm  x  I  mm 

Total  loss  ( X 10^  J/m^) 
u  =  (pcv/Ae)(dr/d/),_,o 

19.2 

19.9 

19.7 

£-£  hysteresis  loss  (X 10^  J/m^) 

18.5 

17.8 

17.4 

it(r)(W/(m-K)) 

38.4 

39.2 

34.1 

out  without  external  stress  applying  to  the  samples,  it  can  be 
concluded  that  under  the  stress-free  condition,  heat  generation 
in  piezoelectric  actuators  is  mainly  due  to  ferroelectric  hystere¬ 
sis  loss,  and  mechanical  loss  is  negligibly  small. 

Regarding  the  heat  generation  from  piezoelectric  vibrators, 
Hirose  et  aL  reported  intriguing  results.*^  When  a  piezoelectric 
vibrator  is  driven  at  the  resonance  point,  the  main  contribution 
to  the  total  loss  is  from  the  mechanical  loss  under  a  small 
vibration  velocity  level,  while  the  contribution  from  the  dielec¬ 
tric  loss  increases  drastically  above  a  cenain  critical  vibration 
velocity  such  as  0.2-0.3  m/s,  where  the  actual  heat  generation 
starts  to  be  observed.  In  our  experiment,  the  driving  frequencies 
are  far  below  the  resonances.  (For  example,  the  resonance  fre¬ 
quency  of  sample  7  mm  X  7  mm  X  2  mm  is  236  kHz,  while  the 
test  conditions  are  below  2  kHz.) 

By  determining  the  time  constant  t,  the  overall  heat  transfer 
coefficient  k{T)  can  be  calculated  through  Eq.  (8).  The  values 
of  k{T)  for  various  actuators  while  driving  under  3  kV/mm. 
3(X)  Hz  are  provided  in  Table  II.  The  k{T)  as  functions  of 
frequency  and  electric  field  are  shown  in  Figs.  12  and  13.  It  is 
seen  from  Fig.  12  that /:(!)  is  relatively  insensitive  to  frequency 
under  the  low  electric  field  level,  while  it  increases  gradually 
with  increasing  frequency  when  the  applied  electric  field  level 
is  higher  (>1  kV/mm).  From  Fig.  13.  it  is  seen  that  k{T) 
increases  with  increasing  the  applied  electric  field.  The  increase 
in  k(T)  may  be  caused  by  the  convection  increase  due  to  the 
vibration  level  (i.e..  vibration  velocity)  change  with  E  and /.  If 
we  summarize  the  frequency  dependence  of  u  and  /.(D,  we  can 
see  that  u  is  insensitive  to  frequency,  k{T)  is  also  relatively 
insensitive  to  /  at  the  low  electric  field  level,  but  gradually 
increases  with  increasing  frequency  at  higher  electric  field  lev¬ 
els  (>1  kV/mm),  Thus,  according  to  Eq.  (9),  temperature  rise 
is  almost  linearly  proportional  to /under  the  low  electric  field 
level,  while  at  higher  electric  field  levels  (>l  kV/mm),  the 
relationship  between  AT  and  /  deviates  from  the  linear  trend 
due  to  the  increase  of  k{T)  at  higher/.  This  analysis  coincides 
with  the  experimental  results  shown  in  Fig.  3.  For  the  electric 
field  dependence,  ii  is  proportional  to  £%  while  k{T)  increases 
significantly  with  increasing  £.  That  is  why  AT  is  nearly  pro¬ 
portional  to  £“  under  low  electric  field  as  shown  in  Fig,  5, 
while  deviating  from  this  square  relationship  at  a  higher  electric 
field  level. 


Fig.  12.  Overall  heat  transfer  coefficient  k(T)  as  a  function  of  fre¬ 
quency  (data  from  the  actuator  with  dimensions  of  7  mm  X  7  mm  X 
2  mm). 


Fig.  13.  k{T)  as  a  function  of  applied  electric  fields  (400  Hz.  data 
from  the  actuator  with  dimensions  of  7  mm  X  7  mm  X  2  mm). 


Equation  (9)  can  be  used  to  evaluate  the  temperature  rise  of 
piezoelectric  actuators  at  the  undamped  state  and  high  field 
(higher  than  the  threshold  field  at  which  domain  wall  motions 
occur)  applications  providing  that  the  equilibrium  temperature 
is  lower  than  the  Curie  point.  By  simply  measuring  the  P-E 
hysteresis  loss  of  the  actuator,  the  temperature  rise  can  be 
estimated  by  ufvJkiDA,  where  k(T)  can  be  obtained  from 
Figs.  12  and  13.  if  the  sample  has  a  size  similar  to  the  actuator 
in  the  figures.  If  the  actuator  has  a  significantly  different  size 
from  the  one  in  Figs.  12  and  13,  the  at  £  =  0,  /  =  0 
can  be  estimated  through  the  empirical  correlation  for  free 
convection  described  in  Ref.  20.  But  notice  that  is  frequency 
and  electric  field  dependent.  If  the  actuator  is  used  under  a 
clamped  state  (as  in  most  applications),  the  heat  conduction 
through  the  sample  surface  is  increased,  and  the  mechanical 
loss  may  need  to  be  considered. 

V.  Conclusion 

Heat  generation  in  multilayer  piezoelectric  actuators  was 
studied.  Under  the  stress  free  condition,  heat  generation  is 
mainly  due  to  ferroelectric  hysteresis  loss,  not  to  mechanical 
loss.  The  temperature  rise  of  the  actuator  can  be  estimated  by 
the  expression  AT  =  ttf\Jk(T)A,  The  proposed  analytic 
method  is  useful  for  the  design  of  multilayer  and  other  high- 
power  actuators  from  the  viewpoint  of  heat  generation. 
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ABSTRACT: 

The  moonie  actuator  is  a  versatile  performer  that  fills  the  gap  between  the  multilayer 
actuator  and  the  bimorph  actuator.  This  papa  describes  the  effect  of  geometrical  change  in  the 
endcap  on  the  moonie's  actuator  characteristics.  The  effect  of  cavi^  itoensions,  which  are  the 
key  parameters  for  transformation  and  amplification  of  radial  motion  into  axial  motion,  on  the 
characteristics  of  the  moonie  actuator  were  investigated.  Actuator  properties  of  the  newly 
designed  cymbd  actuator  were  investigated  with  different  ceramic  driving  element 
Characteristics  of  the  different  composite  ceramic  actuator  designs  were  compared  based  on  their 
displacement  response  time  and  generative  forces. 

1.  INTRODUCTION 

Piezoelectric  ceramics  are  widely  used  as  a  sensor  many  years.  In  recent  years,  there  are 
lots  of  attention  for  the  actuator  application  of  the  piezoelectric  ceramics.  Piezoelectric  ceramic 
exhibits  Poisson’s  effect;  expansion  in  longitudinal  direction  and  shrinkage  in  transverse 
direction,  when  an  electric  field  is  applied  pa^el  to  the  polarization  direction.  There  are  two 
actuator  designs  which  use  this  phenomenon  to  enhance  the  displacement  for  practical 
applications:  multilayer  ceramic  actuators  with  internal  electrodes  and  bimoph  actuators.  The 
advantages  of  the  multilayer  actuator  are  its  generative  force  and  quick  response  sp^. 
However,  small  displacement  and  high  capacitance  of  multilayer  actuators  make  them  impractical 
for  certain  applications.  On  the  other  hand,  bimorph  actuators  show  larger  displacement 
However,  low  generative  forces  and  small  displacement  are  their  disadvantages.  There  is  a  need 
for  an  actuator  with  intermediate  level  displacement  and  generative  force.  The  moonie  and 
cymbal  type  of  flextensional  and  roto-flextensional  actuators  were  designed  to  fill  the  gap 
between  multilayer  and  bimorph  actuators  (1*3].  The  ceramic-metal  composite  moonie  and 
cymbal  actuators  consist  of  ceramic  driving  element  epoxy  or  metal  alloy  bonding  layer  and 
metal,  polymer  or  glass  endcaps.  Piezoelectric,  electrostrictive,  artd  aiitiferroelectric-feTioiElectric 
type  of  ceramic  materials  can  be  used  as  a  driving  element  The  basic  configuration  of  the 
m^nie  and  cymbal  actuators  are  shown  in  Hgure  1.  The  metal  endcaps  serve  as  mechanical 
transformer  for  converting  and  amplifying  the  lateral  motion  of  the  coamic  into  a  large  axid 
displaoeroent  normal  to  the  endcaps.  Bodi  the  d3i  (sd32}  d33  coefficients  of  the  piezoelectric 

ceramic  contribute  to  the  axial  displacement  of  the  composite. 
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The  objective  of  this  study  is  to  evaluate  the  actuator  characteristics  of  the  moonie  and 
cymbal  transducers.  Cavity  beneath  the  endcaps  plays  an  important  role  for  transforation  and 
amplification  of  lateral  motion  into  axial  motion.  The  effect  of  cavity  dimensions  on  the 
chmcteristics  of  the  moonie  were  studied  systematically  altering  the  cavity  diameter  and  cavity 
depth  while  keeping  the  endcap  thickness  constant 


2.  EXPERIMENTAL  PROCEDURE 


The  composite  actuators  were  made  from  electroded  PZT*5A  (Niobium  doped  Lead 
Zirconate  Titanate).  0.9PMN-0.1PT  (Lead  Magnesium  Niobate-Lead  Titanate)  and  PNZST 
(Ptio.» Nbo.Ol[(ZtO.«Sno.4)o.»lTlo.l)6lo?803) (12-7  mm  in I  mm  m 
metal  endcaps  (12.7  mm  in  diameter  with  thickness  ranging  from  0.3  to  3.0  mm).  The  endcaps 
for  moonie  were  machined  from  brass  with  composition  30%  Zn-70%  Cu.  A  punch  die 
designed  was  used  for  the  fabrication  of  the  cymbal  endcaps.  The  ceramic  disk  and  the  endcaps 
were  bonded  together  around  the  circumference  with  Emerson  &  Cuming  epoxy  (Eccobond  45), 
taking  care  not  to  fill  the  cavity  or  short  circuit  the  electrodes.  The  displacement  of  a  composite 
actuator  at  0  01  Hz  was  measured  with  a  Linear  Voltage  Differential  Transducer  (LVDT)  having 
a  resolution  of  approximately  0.05  \m.  Fastest  response  speed  of  the  actuators  was  calcula^ 
from  the  resonant  frequency  data  which  was  obtained  with  a  Hewlett-Packard  Irapedance/Gain- 
Phase  Analyzer  HP-4194A. 


cza  Brass  Endcaps  CZJ  Piezoelectric  Ceramic  —  Bonding  Layer 

Dimensions:  (all  in  mm)  Endcap  diameter,  din=  12.7  mm,  PZT  diameter;  dn=  12.7  mm 
Cavity  diameter;  dc=  variable;  3.0, 5.0, 7.0, 9.0  /  Cavity  depth,  h=  variable;  0.3, 0.5, 0.7, 0.9 
Metal  cap  thickness;  tni=l.0,  PZT  thickness;  tp=  1.0  Bonding  layer  thickness;  tb=  0.020 

Figure  1  The  cross  sectional  view  of  the  moonie  (on  the  left)  and  cymbal  (on  the  right) 
actuator.  Arrows  show  the  displacement  direction 


3.  EVALUATION  OF  THE  EFFECT  OF  GEOMETRY  ON  THE  CHARACTERISTICS  OF 
THE  MOONIE  TRANSDUCER 

By  altering  the  geometric  parameters,  it  is  posable  to  change  the  effective  compliance  of  the 
composite  moonie  structure  and  tailor  the  desired  actuator  properties.  The  effect  of  the  endcap 
thickness  on  the  performance  of  the  moonie  actuator  has  been  studied  expenracn^y^ 
correlation  with  FEA  methods.  The  results  has  been  reported  in  earlier  publicattons  {3-7].  The 
displacement  shows  an  inverse  relation  with  endcap  thickncM.  The  di^la^ment  of  me  inoome 
actuator  increases  with  decreasing  thickness,  because  of  the  increased  flexibibty  of  the  endt^ps. 
A  moonie  actuator  with  0.3  mm  thick  brass  endcap  and  with  the  fixed  cavity  diameter  ratio  at 
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80%,  which  is  the  ratio  of  the  inner  diameter  of  the  endcap  to  the  outer  diameter,  exhibits  a  22 
um  displacement  under  an  applied  electric  field  of  1.0  kV/mm.  Generative  force  of  the  moonie 
actuator  decreases  drastically  with  decreasing  endcap  thickness.  Although  the  flexural  motion 
and  displacement  are  largest  at  the  center  of  the  endcaps,  the  generative  force  is  the  lowest  at  that 
point.  Displacement  decreases  dramatically  when  moving  from  the  center  to  the  edge. 
Conversely,  generative  force  increases  when  moving  from  the  center  to  the  edge,  where  it 
approaches  that  of  ?ZT  ceramic.  Therefore,  it  becomes  possible  to  tailor  the  desired  actuator 
properties  by  changing  the  contact  area  over  the  brass  endcaps. 

Brass  endcaps  with  3,  5,  7,  and  9  mm  cavity  diameters  Md  0.3, 0.5, 0.7,  and  0.9  mm 
crescent  shaped  cavity  depths  were  fabricated  by  matrix  combination  while  keeping  the  thickness 
constant  at  1.0  mm,  and  diameter  fixed  at  12.7  mm.  Figure  2  shows  the  relation  between 
displacement  and  cavity  diameter.  The  displacement  amplification  is  exponentially  related  to  the 
cavity  diameter  of  the  endcaps  (at  constant  cavity  depth)  and  linearly  related  to  the  cavity  depth 
(at  constant  cavity  diameter).  Fastest  response  time-cavity  diameter  relations  at  various  constat 
cavity  depths  of  the  moonie  actuators  is  shown  in  Figure  3.  Response  lime  increases  rapidly 
with  increasing  cavity  diameter  and  decreases  slightly  with  increasing  cavity  depth. 


2345*78>I» 
Cavity  Diameter  (mm) 


234  56789  10 

Cavity  Diameter  (mm) 


Beuie2.  The  displacement-cavity  diameter  relation  of  the  moonie  acmator  (on  tte  left) 

Rgure  3.  Fastest  response  time-cavity  diameter  relation  of  the  moonie  actuator  (on  the  right) 

The  effective  coupling  factor,  keff.  can  be  derived  from  the  series  (resonant)  and  p^el 
(antiresonant)  resonant  frequencies  of  the  transducer  [9].  The  effective  coupling  factor  for  the 
first  resonance  mode  of  moonie  actuators  with  varying  cavity  depth  ^d  cavity  diameters  were 
calculated  from  their  admittance  spectra.  Figure  4  shows  the  calculated  effecuve  coupling 
keff  of  the  moonie  actuator  as  a  function  of  the  cavity  diameter  at  vanous  cavity  depths,  ^e 
planar  coupling  factor  of  PZT-5A  disc  with  12.7  mm  in  diameter  and  1^  mm  in  thicks  u 
around  0.65.  ’Ae  coupUng  factor  decreases  to  0.54  because  of  the  load  effect  when  the  IOT-5A 
ceramic  disc  is  sandwiched  between  brass  endcaps.  The  moome  actuator  with  a  cavity  di^erer 
of  3  0  mm,  and  a  cavity  depth  of  0.9  mm  has  the  highest  effecuve  coupling  factor  of  0.42  As 
Uie  cavity  diameter  inercas^,  the  effective  coupling  f^tor  decrcpid  linearly  down  *0  0.1 1  for  a 
cavity  diameter  of  9.0  mm  with  a  cavity  depth  of  0.3  ram.  Effecuve  coupling  factor  depends 
heavily  on  the  cavity  diameter  and  shows  low  sensitivity  to  the  cavity  depth. 
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Rgure  4  Effect  of  cavity  diameter  on  the  effective  coupling  factor,  keff,  of  the  raoonie 
tiMsducer  at  various  cavity  depths,  (on  the  left). 

Figure  5  Displacement  hysteresis  of  different  endcap  designs,  (on  the  right). 


4.  CYMBAL  ACTUATOR 

Finite  element  analysis  in  correlation  with  experimental  results  showed  that  there  is  a  high 
compressive  stress  on  the  endcaps  just  above  the  bonding  layer.  Mechanical  losses  due  to  this 
stress  concentration  reduces  the  effective  force  transfer  and  from  the  PZT  to  the  endcap.  A  new 
actuator  design  with  truncated  conical  endcaps  was  developed,  named  the  “Cymbal”  (3-4).  The 
cymbd  actuator  exlubits  higher  displacement  with  larger  generative  force  along  with  cost 
effective  manufacturing.  Displacement  hysteresis  of  the  different  endcap  designs  are  shown  in 
Figure  5.  With  the  cymbal  endcaps  displacement  values  enhanced  twice  as  much  as  moonie’ 
value.  Displacement  characteristics  of  the  cymbal  actuators  with  different  ceramic  materials  as 
the  driving  unit  are  shown  in  Figure  6.  Cymbal  actuator  with  PZT-5A  exhibits  linear 
displacement  The  cymbal  actuator  with  the  P\fl4-PT  type  ceramic  driving  element  shows  larger 
displacement  with  lower  losses.  However,  it  exhibits  non-linear -displacement  characteristics 
according  with  the  relaxor  characteristics  of  the  PMN-PT  ceramic.  PNZST  is  an 
antifenoelectric-ferroelectric  type  phase  transition  materiaL  Unique  properties  of  this  material  is 
that  the  volumetric  expansion  under  applied  electric  field.  With  cymbal  endcap  design  this 
volumetric  expansion  is  convened  to  the  negative  axial  displacement 

5  SUMMARY 

Displacement  of  the  moonie  actuator  is  approximately  inversely  proportional  to  endcap 
thickness,  exponentially  proportional  to  the  cavity  diameter,  and  linearly  proportional  to  the 
cavity  depth.  Response  time  of  the  moonie  actuator  increases  with  increasing  cavity  diameter. 
Moonie  shows  position  dependent  behavior  with  maximum  displacement  at  the  center  and 
maximum  force  at  the  perimeter.  Desired  actuator  properties  can  te  easily  t^ored  by  changing 
the  contact  area  over  the  endcaps.  The  moonie  actuator  can  give  displacements  of  20  to  greater 
than  100  itm,  inaximum  generative  force  of  3  N  to  several  N's,  and  response  time  of  20  to  100 
)isec.  By  altering  some  of  the  dimensions  on  the  moonie  design,  it  is  easy  to  tailor  an  actuator 
with  desired  properties.  Cymbal  actuator  shows  twice  as  big  as  displacement  than  moonie  for 
the  studied  size  and  much  higher  generative  force  due  to  increased  contact  surface. 
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Figure  6.  Displacement  hysteresis  various  ceramic  discs  and  cymbal  actuators  with  this  ceramiM 
as  driving  element  a)  PZT-SA,  b)  Cymbal  with  PZT-5A.  c)  PMN-PT.  d)  Cymbal  with  PMN-PT 
e)  PNZST,  0  Cymbal  with  PNZCT 
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Several  features  of  the  various  solid  state  actuator  designs  are  listed  in  Table  I.  Similar 
dimensions  for  each  actuator  were  selected  to  make  a  fair  comparison.  Advantages  of  the 
moonie  and  particularly  the  cymbal  actuators  are  the  easy  tailoring  the  desired  actuator  properties 
by  altering  the  cavity  size  and  endcap  dimensions  with  cost  effective  manufacturing. 


Table  1  Comparison  of  the  Solid  State  Actuator  Designs 


mmmm 

BintiiTU 

Cymbal 

m^nrnrmi 

References 

[10] 

[101 

[11] 

[3^1 

[5-8] 

Dimensions 

5x5x12.7 

(UWxT) 

12.7x10x1 

(LxWxT) 

d>  12.7  mm 

T=  0.5  mm 

d>  12.7  mm 

T=  1.7  ram 

12.7  mm 

T=  1.7  mm 

100 

100 

450 

100 

100 

Displaceraent 

(pm) 

10 

35 

20 

40 

20 

Generative 
Force  (N) 

900 

0.5-1 

1-3 

15 

3 

Fastest 

Response 

Time(psec) 

1-5 

100 

100 

5-50 

5-50 

Cost 

mm 

med. 

med. 

low 

low 

ACKNOWLEDGMENTS 

The  authors  would  like  to  express  their  gratitude  for  the  support  to  the  following  agencies 
and  organizations:  Office  of  Naval  Research  Contract  no.  N00014-92  J  1510,  National  Science 
Foundation  Grant  no.  DMR-9223847,  Spanish  Science  Ministry  (CICYT  MAT94-807  and 
IX3ICYT  PR94-028),  Turkish  Science  and  Technology  CouncU  (TUBITAK)  and  Middle  East 
Technical  University,  Ankara-TURKEY. 

REFERENCES 

[1]  R.E.Newnham.  Q.C.Xu  and  S.  Yoshikawa,  U.S  Patent  #  4,999,619  (March  12,1992). 

[2]  Y.Sugawara,  ICOnitsuka,  S.Yoshikawa,  Q.Xu,  R.E.Newnham,  and  K.Uchino, 
J.Am.Ceram.Soc.,  voL  75,  no.4,  pp.  996-999,  (1992). 

[3]  A.  Dogan  Ph.D.  thesis  The  Pennsylvania  State  University  (PSU),  1994 

[4]  A.  Dogan  and  R.E.  Newnham  “  Flextensional  Cymbal  Transducei^’  PSU  Invention 
disclosure  No.  94-1395,1994 

[5]  A.Dogan,  Q.C.  Xu,  K.Onitsuka,  S.Yoshikawa,  K.Uchino  and  R.E.Newnham, 
Ferroelectrics,  vol.  156,  pp.  1-6,  (1994) 

[6]  R.E.Newnham,  A.Dogan,  Q.C.  Xu,  K.Onitsuka,  S.Yoshikawa,  K.Uchino  IEEE  - 
UFFC ,  Baltimore  1993 

[7]  K.Onitsuka,  A.Dogan,  J.Tressler,  Q.C.Xu,  S.Yoshikawa,  and  R.E.Newnham, 
Ferroelectrics,  vol.  156,  pp.  37-42,  (1994) 

[81  A.Dogan,  S.Yoshikawa,  K.Uchino  and  R.E.Newnham,  lEEE-UFFC  Proceedings 
Vol  U  pp  935-939,  France  1994 

[9]  Mason,  W.P.,  Physical  Acoustics  Vol  1  Part  A,  Academic  Press,  New  Yoii:,  (1964) 

[10]  Information  was  extracted  from  Tokin  Co.  product  catalog  1993 

[11]  Haertling  G.  Bull  of  Am.  Ceramic  Soc.  Vol  73  No  1,  pp  93-96, 1994  (1994) 


